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SUMMARY 


The aerodynamics of three-dimensional flow over swept wings 
with sharp leading edges is discussed, and an attempt is made to 
establish the character and effects of the spiral-conical vortex 
formed at the sharp leading edge of a swept wing 

Potential flow equations for a flat plate, with the addition of a 
point vortex near the leading edge, are used to determine the 
streamwise flow conditions at various span stations. For a given 
vortex location, the vortex strength and airfoil circulation are 
determined by satisfying the Kutta conditions at the leading 
and trailing edges. The lift due to circulation is modified by 
the presence of the leading-edge vortex 

The spanwise variation of the lift coefficient is obtained by 
combining, at a determined location, the point vortex with the 
potential flow over a flat plate. Finally, the calculation of span 
load distribution, allowing for the presence of the rolled-up free 
vortex at the leading edge, is discussed. 

The theory presented herein is nonlinear, and the expression 
for lift due to circulation is shown to be nonlinear 


List OF SYMBOLS 


Cartesian coordinates, physical plane 

Cartesian coordinates, circle plane 

n + 12&, circle plane 

n — 1&, circle plane 

radial distance, circle plane 

angle between 7 and ¢ in circle plane 

angle of attack 

x + ty, physical plane 

free vortex intensity 

vorticity required to satisfy the Kutta condition at 
the trailing edge 

bound vorticity 

velocity component in direction of n, circle plane 

velocity component in direction of £, circle plane 

velocity component in direction of x, physical plane 

velocity component in direction of y, physical plane 

force components 
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air density 

drag and lift components 

section lift coefficient 

velocity along free vortex core 

angle between axis of free vortex core and its projec- 
tion in plane of wing 

sweep angle of projection of axis of free vortex core in 


plane of wing 


Subscripts 


= core of the free vortex 
= reattachment point of the stagnation streamline over 


the vortex 


INTRODUCTION 


—_ AERODYNAMICS OF THE FLOW over a wing with 
a round leading-edge airfoil section can be con 
sidered to be completely nonviscous at low angles of 
attack because the regions of flow separation are limited 
and the potential flow boundaries are essentially the 
surface of the wing. With a sharp leading-edge airfoil 
section, flow separation takes place right at the leading 
edge, and the potential flow boundaries are consider 
ably modified over a large portion of the airfoil surface. 
For swept wings with a sharp leading-edge airfoil sec- 
tion, the flow is further complicated by the presence 
of a spiral-conical vortex at the leading edge of the 
wing where separation takes place. With the intro- 
duction of the conical vortex into the potential flow 
picture, the resulting equations are nonlinear, and it 
will be shown later that the modified lift due to circu- 
lation is nonlinear. 

The conical vortex at the wing leading edge was shown 
by Winter! in 1936. Somewhat later,” tests on a large- 
scale delta wing with sharp and round leading edges 
showed the importance of the effects of the conical 
vortex on wing lift. Studies of a 47.5 degree swept- 
back wing*® show that both the chordwise and spanwise 
load distributions for a sharp leading-edge wing cannot 
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be determined by conventional methods, such as 


Multhopp, and Falkner 
The downwash characteristics for 


W eissinger, even for very 
low angle of attack. 
the same wing were reported,’ and it was shown that 
two distinct trailing vortex patterns exist, clearly not 
a conventional result. Studies® of the conical vortex 
formed at the sharp leading edge of a delta wing were 
also made by the NACA. 


ration region and consequently the size of the conical 


The magnitude of the sepa 


leading-edge vortex were shown to be so large that fur 
ther neglect of this phenomenon, when calculating span 
load distribution or downwash, is inconceivable. 

To determine the shape, vorticity distribution, and 
location of the leading-edge conical vortex with its 
trailing vortex system is a formidable mathematical 
problem. Kaden* and Westwater’ have studied the 
rolling-up of the trailing vortex sheet as it leaves the 
trailing edge of a straight wing. This problem is com- 
paratively simple compared to the rolling-up of the vor 
tex sheet at the leading edge, and yet no real mathe 
matical solution exists for it. Recently, Stern*® has 
obtained a mathematical solution for the problem of 


the rolling-up of a trailing vortex sheet, and it may be 





l 


possible to extend his work to inciude the rolled-up 
vortex at the leading edge of a sweptback wing at a 
later date. 

Legendre? analyzes the leading-edge vortex of a highly 
sweptback wing by taking two-dimensional transverse 
cuts across the wing. The problem then reduces to 
flow across a flat plate in a manner similar to that 
worked out by Riabouskinsky'’ (see example 3, p. 254 
In this problem, stationary vortices form behind the 
flat plate at specified points in such a way that the ve 
locity is finite at the edges of the plate. Actual flow 
studies'' have been made on a delta wing on the upper 
surface of the wing at an angle of attack of 17°, as 
shown in Fig. 11 of the reference. The vortex fila 
ments at different points along the span are shown, and 
it can be seen that the filaments are not straight lines 
but helicoidal. It is also clearly shown that vorticity 
is constantly changing spanwise, which offers experi 
mental proof that the free vortex intensity varies along 
the leading edge of the wing. 

In this paper an attempt is made to establish th 
character of the spiral-conical vortex formed at the 


sharp leading edge of a swept wing and to replac 
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the viscous separated flow by a potential flow in two 
dimensional streamwise planes. The spanwise varia 
tion of the lift coefficient is obtained by combining, 
at a determined location, the point vortex with the po- 
tential flow over a flat plate. At some spanwise sta 
tion the stagnation streamline over the point vortex 
will fail to reattach on the airfoil chord. At this sta 
tion, and for all stations outboard of this station, the 
flow is completely separated. The calculation of span- 
wise variation of the lift coefficient will be carried out 
only to the station where complete separation occurs. 
In order to complete the spanwise variation of lift co 
efficient, recourse must be had to experimental data be- 


yond the stall point. 


PHYSICAL INTERPRETATION OF SEPARATION VORTEX 


The general nature of the leading-edge vortex is indi- 
‘ated in Fig. 1. On the right is a plan view of a swept 
wing illustrating the formation of a leading-edge vortex 
at the apex and its growth in size as it progresses out 
board and rearward. <A view from the rear of the wing, 
looking forward, shows that two distinct vortex pat 
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terns exist immediately bebind the wing, and that both 
On the left 


1 is the same type of flow for a very low 


are removed from the plane of the wing 
side of Fig. 
aspect ratio rectangular wing A comparison is 
drawn, inasmuch as for large sweeps the leading edge 
of the swept wing tends toward similarity with the tips 
of the rectangular wing. The vortex for the swept wing 
might be called a trailing vortex formed at the leading 
edge, whereas the one at the left 1s the trailing vortex 
normally referred to. A study of the flow pattern be 
hind a swept wing was made? and the results are re 
produced in Fig. 2. 

A somewhat more detailed insight into the formation 
of the leading-edge vortex can be obtained from Fig. 3 
In the center is a plan view showing one streamline of 
the fluid flow. At the right side is a two-dimensional 
section, and the same streamline is illustrated (adjacent 
to the streamline stagnating on the lower surface) going 
under the airfoil, reversing its direction, and flowing 
up and over the leading edge. A boundary layer forms 
from the lower surface stagnation point to the sharp 
leading edge where it separates. The pressure gradient 


from the stagnation point to the leading edge is nega 
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tive so that the boundary layer at the point of separa- 
tion is very thin and probably can be treated as a vor- 
tex sheet. The vortex sheet is curled up by the main- 
stream, as well as by the induced velocities of the ad- 
jacent vortices situated on the free vortex sheet, and 
then travels spanwise. All the vorticity in the rolled- 
up vortex at the leading edge is initiated in the small re- 
gion indicated in Fig. 3. The individual vortex ele- 
ments making up the vortex sheet can be considered to 
be part of a vortex filament starting at the sharp lead- 
ing edge as indicated. In this manner, it is seen that 


a bundle of vortex filaments can be considered to 
emanate from the wing leading edge (left wing panel of 


) the vortex 


rig. 3). For small angles of attack (2-3 
filaments can be drawn to the wing tip where they are 
then washed aft by the main stream flow. Aft of the 
bundle of vortex filaments, the flow has reattached to 
the surface. For higher angles of attack, complete 
separation takes place in the outboard regions of the 
wing. The main potential flow does not appear to re- 
attach behind the rolled-up vortex sheet, and the vortex 
filaments must be considered to turn aft before reaching 
the wing tips. This does not actually happen as 
sharply as is indicated in Fig. 3; the separation vortex 
will proceed along smooth curves. 


REPLACEMENT OF SEPARATION VORTEX BY POTENTIAL 
VORTEX 

At various spanwise stations the vorticity is con- 
sidered to be concentrated at the center of gravity of 
the vorticity, and a line connecting these points is as- 
sumed to be a straight line which is swept back from 
the leading edge and bends aft near the wing tips. A 
cut normal to the axis of the vortex will yield a single 
Sections normal to the leading edge are 
The stagnation streamline in the 


point vortex. 
illustrated in Fig. 4. 
potential model separates the flow into two potential 
regions, and the outer region is similar to the potential 
portion in the real flow. The point vortex is stationary 
in this plane. The resultant velocity at the point vor- 
tex at any span station is considered to be along the 
vorticity c.g. line. The velocity along this line is illus- 
trated in Fig. 5, and the three components along the 
coordinate axis are shown. For the purpose of this 
paper, the spanwise velocity component, v,,, is not con- 
sidered in determining the flow in two-dimensional 
streamwise planes. The component of velocity normal 
to the wing plane, v,,, will be considered in conjunction 
with the streamwise component, ?;,. 

The sketch on the right side of Fig. 5 illustrates a two- 
dimensional streamwise cut showing a point vortex 
near the wing surface. In this view it must be remem- 
bered that the vorticity is only a component of the total 
vorticity. The resultant vorticity vector is assumed 
to be along the vorticity-center of gravity line, and the 
vector direction would be toward the apex of the wing. 

Having assumed the character of the conical leading- 
edge vortex to be as just described, flow only in two- 
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dimensional streamwise planes will be considered. |, 
is assumed that the location of the center of gravity of 
the leading-edge vorticity across the span of the wing 
is tantamount to determining the greater part of the 
mutual influence effects between neighboring span sta 
tions where a leading-edge vortex exists. 

Fig. 6 illustrates the notation used in determining the 
two-dimensional streamwise 


flow over a 


transforming from the circle plane to the physical 


section by 
plane. The free vortex, A, which is a component oj 
the resultant vorticity of the leading-edge vortex, re 
quires two image vortices, one at the center of the circle 
and one inside the circle. The total vorticity at the 
center of the circle and the vorticity, A, of the fre 
vortex is adjusted to satisfy the Kutta conditions at 
the leading and trailing edges. Considering that the 
Kutta conditions are satisfied in streamwise planes 
rather than in planes normal to the leading edge is a 
further assumption that should be investigated more 
fully. 
probably justified, as will be illustrated later in the 


For sweep angles up to 45°, this assumption is 
paper. 


SPANWISE Lirt DISTRIBUTION 


For a given angle of attack, the two-dimensional 
streamwise flow is fully determined once the position 
of the vortex is known and Kutta conditions at the 
leading and trailing edges are assumed. This can be 
seen from Eqs. (5), (6), (7), (8), (9), and (10) of the 
next section where A, I’, %»,, vz,, ¢z,, and v,, are all shown 
to be functions of r, and 6, and hence x, and y,. 

The variation of C, from Eq. (11) of the next section 
is shown in Fig. 7. The sweep of the leading-edge 
vortex, A,, which yields the value of x, for any par 
ticular wing has been shown*: * to be primarily a func 
tion of the sweep of the leading edge and the angle oi 
attack, and not to vary with planform. 

The value of A, has been determined for the wing 
investigated by the NACA® by examining the pressure 
distribution diagrams and, in accordance with further 
NACA data,°® the location of the vortex is given by 
the points of minimum pressure. For an angle of attack 
of 10.2°, the sweep of the leading-edge vortex, A,, is 
54.3° 
vortex remains close to the leading edge for the inboard 
Outboard, the vortex bends 


The pressure diagrams also indicate that the 


50 per cent of the span. 
aft and stalled-type pressure distributions are obtained 
Accordingly, calculations will be made only out to 50 
per cent of the wing semispan to illustrate the validity 
of the theory presented herein. 

The wing investigated in this paper corresponds to 
that shown in Fig. 9 where the sweep of the quarter 
chord line is 45° and the airfoil section is a circular 
arc. In making the calculations, it was necessary to 
take the sweep of the leading-edge vortex in planform, 
A,, from experimental data. The value of A, for am 
angle of attack of 10.2° is 54.3°. Referring to Fig. 4, 
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itcan be seen that 


V, = w,sind 


Ur, = W, cos 6 sin A, 


Taking the ratio of the two quantities, one obtains 


vy,/Ur, = tan 6/sin A, 


y 


It will be noted that in the latter equation it is not 
necessary to know the resultant velocity of the free 
vortex. This equation constitutes the starting point. 
As stated above, it is necessary to resort to experimental 
lata for the value of A,. Knowing A,, it is then possible 
to determine x,. A value of 6, is then assumed which 
fixes y,. Knowing x, and y,, 7, and @, can be obtained 
from the transformation 


s = [+ (1/f) 
x, = [r, + (1/r,)] cos 8, 
y, = [r, — (1/r,)] sin 6, 


With these data, A and [ can be computed from Eqs. 
Knowing A and I, u,, and vg, can be ob- 
Finally, u;, and v,, 
The ratio of 


d) and (6). 
tamed from Eqs. (7) and (8). 
are determined, using Eqs. (9) and (10). 

to u;, computed in this manner should be equal to 
the value determined previously by tan 6/sin A,. If 
not, it is necessary to try another value of 6, until the 
The to the correct 


process converges. convergence 


value is quite rapid. 


FLOW IN THIS REGION 
FLOWS SPANWISE 


VORTEX FILAMENT 





POTENTIAL FLOW 
STEAMLINE 











SECTION A-A 





It may be of value to define stall at any given station 
as that condition where the reattachment of the stag 
nation streamline over the vortex fails to exist on the 
airfoil. 
stalled section C; should be used. 


When this condition occurs, corresponding 
To determine these 
values, recourse must be had to experimental data. 
Fig. S shows the calculated reattachment points as a 
The 


This chart was calculated 


function of angle of attack and the value of x,. 
variation with y, is small. 
on the following basis: Eq. (1), when differentiated 
with respect to ¢, yields a quartic in ¢, the four roots 
of which correspond to the four stagnation points in 
the circle plane. The Kutta conditions at the leading 
and trailing edges of the wing in the circle plane yield 


the following two conditions: 


me 
= ( 
dé J, Ly 


dF(¢) 
= 0) 
dé F 1 


These conditions produce finite velocities at the leading 
and trailing edges of the wing in the physical plane. 
Having introduced the Kutta conditions at these two 
points in the circle plane, the quartic reduces to (f — 1) 
(¢ + 1) (€¢2 + Aé + B) = 0, the two unknown roots of 
which correspond to (1) the conventional stagnation 
point on the lower surface of the airfoil, and (2) the up- 
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per surface stagnation point of the stagnation stream- 
line over the free vortex. 

For an angle of attack of 10.2°, the value of x, for a 
value of xz of 100 per cent (stall condition) is approxi 
mately 44 per cent of the chord. This corresponds to 
aspan station of about 70°,. This value is high, and at 
this time can only be used as a rough indication of sec- 
tion stall and the point at which the free leading-edge 
vortex leaves the wing and is washed aft by the main 
stream flow. When this happens, the longitudinal and 
directional stabilities can be materially affected because 
of the closer proximity of the rolled-up free vortex to 
the tail surfaces. This phenomenon undoubtedly 
accounts for the major portion of the sharp nonlinear 
effects observed in the Cye, and Cy, curves at the 
higher angles of attack. 

The results of the calculations for section lift co 
efficient are shown in Fig. 9 and are compared to a 
Weissinger calculation and the experimental data for a 
wing’ at an angle of attack of 10.2°. It will be noted 
that the lift coefficient steadily increases spanwise up 
to the point where the stagnation streamline over the 
vortex fails to reattach to the airfoil chord, and at a 
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rate than that predicted by Weissinger’s 


For our particular example, separated flow 


greater 
method. 
occurs experimentally at 50 per cent of the wing semi 
span. Outboard of 50 per cent of the wing semispan, 
comparisons are no longer valid because the leading 
edge vortex is bending aft and the sections are exhibit 
ing stall characteristics. It would appear that the 
model used in this paper accounts for a good portion 
of the three-dimensional effects existing for the inboard 


a sweptback wing through the concept of 


portion of 
the three-dimensional free vortex core-axis 


DEVELOPMENT OF EQUATIONS* 


Expressions for Free Vortex Intensity and Bound Vorticity 


Fig. 5 illustrates the position of the free vortex rela- 
tive to the leading edge of the wing and wing chord 


Fig. 6 shows the position of the free vortex 


plane. 
relative to the unit circle in the ¢-plane. 


potential function in the circle plane for an isolated 


The complex 


* The authors wish to express their appreciation of the assist- 
ance rendered by Messrs. M. Stern and M. Vinokur of the Aero- 
dynamics Division, Republic Aviation Corporation 
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r rpe * . . 9) | te | = =. 
and K’ where K’ is the vorticity required to satisfy the a Sale 2(1 + ¢>) (1 + §,) sin a 
‘i = —_— @s i = 
Kutta condition at the trailing edge. t+ ¢, 
Finally, the expression for the free vortex intensity becomes 
é 211 — ¢,2)(1 — ¢,2) sin a 1+ 7,* — 27,2 cos 26, . . 
K = ae = = 2 sin a (0 
Cl = fs) Ge 4 Fe) r,(1 — 7,2) cos 8, 
and for the bound vorticity, 
: —2 sin a (1 + ¢,¢,) 1l+~r,* . 
rT = = —— = — sin @ (6 
Oe + Gs r, COS 6, 
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The velocity of the point vortex in the circle plane is given by 
Un, — Ws. = +=. [F(¢ — 1K,1n (¢ — c,)] 
(dc f 
But 
F(¢) — 1K, In (¢ — § -fe'* — — ¢€ iK, In (¢ + wing 
« ¢ 
Hence, 
ee rs tae ee all 
6 Be = 5 8 iN, — + 21 ( 
{ ¢ T,— I ¢ 
ta 
5 i é ane « l 
u, — WwW: = e” + — 1h + aI 
¢ i l ( 
The component velocities wu, and v; are: 
l K,r ri. 
up. = — cosa Tr cos (a + 26 = sin @ \4 
° r,- — | yr 
, 1 . | ke r 
v% =+snat+ sin (a + 26,) + > - cos 6 (3) 
fe r= r, — 1 r 
To find the velocity in the physical plane, z, we use the transformation : 
— ee oe es 
Ur, — Wy = (uu, — ;,) —+ 1K lim In 
dz t>¢,d2 2-2 
This expression reduces to 
¢ — ¢ 
Ur — Wy, = (u,, — 2 - (kK — 
¢ l (t-— i}? 
The component velocities, “,, and v,,, are: 
Uy, (%o" — COS 20,) — vz, sin 20 r,> sin 30, — (27,* + 7,) sin 6 
Ue = = = K = (9) 
ry,” — 2 cos 26, + 1/r,’ (7,4 — 2r,* cos 20, + 1 
“ > ne DD . o con Sh ‘ ' ‘ 
Uz,(%r — COS 20,) + Uu,, sin 26 _r, cos 36, + (r, — 2r,°) cos @, 
vy = +A (10 
ro” — 2 cos 20, + 1/r,” (r,* — 2r,? cos 20, + 1)? 


Development of Lift Equation 


To determine the lift around the airfoil section, we 
apply Blasius’ theorem: 


1. {( (dF\? 1. ((/dF\*d¢é je 
*, we \dz - * o\d¢é —" 


where the symbols se and 0 under the integral signs 
denote the physical and circle plane, respectively. 


XY —-1Y= 


Let us consider first the integral taken around a large 


contour, C, including the isolated vortex point, ¢,, out- 


from the integranc 


iA 


We evaluate the residue at ¢ = ¢ 


side the circle. For che integral we obtain 


[ (T)F 
dé = 
7 dé dz 


We next proceed to evaluate the integral around the 


trTe'* 


contour, 0, excluding the isolated vortex point outside 
the cirele. 

To do this, we subtract the integral around ¢, from 
the integral around C. 


l 
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The residue is: 


The integral around ® is then 


: e la 
/ te e+ K e” + - 
J 4) cs." 


Applying Blasius’ theorem, we obtain 
Dhl [ 

1pe 
2 JD 


2pm ‘ wT + | -1 + 7,~%e ** 


- 4h 


D—-iL= 


The expression for lift can then be written as follows: 


AERONAUT 


iAr,e 
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ih 
( | 


v7 — J 


sa led Oo Lai ae K tre ** ) 
+ a'r, te + 
Pe l (vt. 


z . resr-A)-—-T, ‘|r,* — r, cos 20, . 
L = —2rp, Tl —|1—-—~7, 2cos2 (a + @4,) — - sin (a + 6,) Sill @ 
\ tAv," —~ F) J(1 — 7,7) cos 6, 
as re-(r — A) -T r, sin 26 ' 
t, * sin 2 (a + 6,) — : cos (a + 9@,) = Sin @ 
r,(%,- —- 1) (1 -- 7,7) cos 8 
[7,4 (1 -- 2 sin? 26,) —- 27,2 cos 20,+ 1] . a 2 sin 26,(r,7 cos 26, — 1) 
——* sin (@ — 6,) sin? a — a r,cos (a — 6,)sin?a@} (11) 
ry,(1 —- 7,7)? cos? 6, (1 — r,7)* cos? 6 
Note: 47.5° Having Circular Arc Airfoil Sections and Equipped With 
Drooped-Nose and Flaps, NACA RM No. L9G15, September 8, 
; e L 1949 
ies 2 ir 4 Lange, Roy H., and Fink, Marvin P., Studies of the Floi 


(1/2)p-12-4 2p 


It is possible to determine the section lift once we 
know the position of the free vortex, 7, and 6,. 
expressions for K and I are also functions of 7, and 0,, see 
Eqs. (5) and (6). 

It is interesting to note that a similar expression for 
the drag can be written which will yield a higher value 
than the one in which no free vortex exists at the lead- 
ing edge of a wing. At angle of attack and for highly 
swept wings where the free vortex is quite intense and 
large, rather large drag penalties will be incurred. 
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Determination of Loads in the Presence of 


Thermal Stresses 


* 


SAMUEL LEVY?t 


National Bureau of Standards 


SUMMARY 


A method is presented for determining loads on an aircraft 


structure when thermal stresses are also present. Equations 
are given for determining the axial load, shear load, and bending 
moment from the output of gages located at specified positions 


on the structure 


(1) INTRODUCTION 


7 DETERMINATION OF AIR LOADS by strain gage 
instrumentation on airplanes in flight has been an 
important method of securing reliable load information. 
With the advent of high speed flight and its associated 
problems of aerodynamic heating, the reliability of this 
method of load determination has been subject to 
question because superposed on the stresses due to load 
will be the thermal stresses from a variable and unknown 
temperature distribution. 

This paper describes a method for determining loads 
when thermal stresses are also present. The method 
takes account of the facts: (1) that the more common 
types of wire resistance strain gages give an indication 
of strain with a change in temperature even when the 
external load is zero; and (2) that the variable temper- 
ature in a highly redundant structure, such as a wing, 
causes internal stresses which are superposed on the 
stresses due to externally applied loads. 

The basis of the formulas developed in this report 
is the same as that used by Gauss for his formulas for 
numerical integration.! In this method when u gages 
are to be used, it is possible to so place them that the 
results will be rigorously true when the stress distribu- 
tion can be described by a polynomial of degree no 
greater than (2x — 1). In those cases where the stress 
in a structure can be closely represented by a poly- 
nomial, the integration of the stress in this manner can 
be expected to give adequate accuracy. Some of the 
formulas given in this paper have the same numerical 


constants as those Gauss obtained. 


(2) TEMPERATURE COMPENSATION OF GAGES 


The output of strain-sensitive wire strain gages at- 
tached to an unloaded member the temperature of 


Presented at the Structures Session, Twenty-Second Annual 
Meeting, IAS, New York, January 25-29, 1954 

* This work was conducted at the National Bureau of Stand- 
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which is changed from one value to another depends 
on: 

(1) The temperature coefficients of expansion of the 
wire and of the member. 

(2) The temperature coefficient of resistance of the 
wire. 

(3) The strain coefficient of resistance of the wire. 

The interaction of these effects was thoroughly dis- 
cussed in 1940 in a report by de Forest and Leader- 
man.” They show that, by suitably choosing a pair of 
wires to be connected either in series or parallel in a 
gage, it is possible to obtain a gage which will be tem 
perature compensated when attached to a member 
made of a specific material. Such a gage will have an 
output proportional to the strain in the member to 
De Forest and Leaderman show, 
Nichrome 


which it is attached. 
for example, that a gage having 0.0001-in. 
and 0.0015-in. Copel wires in series will be temperature 
compensated for use on aluminum alloy if the ratio of 
lengths of Nichrome wire to Copel wire is about 0.73. 

Substantial progress has been made since 1940 in the 
art of temperature compensation by using combinations 
of wires. As a result, such gages are now commer- 
cially available for use at temperatures as high as 
300°F. 

An alternative method of temperature compensation 
mentioned in reference 2 makes use of an auxiliary gage 
attached at right angles to the expected direction of the 
applied load. This gage and the gage placed in the 
expected direction of the load are connected in adja- 
cent arms of a Wheatstone bridge. Changes in re- 
sistance due to temperature are thus balanced out 
in the Wheatstone bridge. An advantage of this 
method is that, by using identical gages, there is no 
output as a result of changes in temperature. Actu- 
ally, small differences in nominally identical gages will 
result in some output. A disadvantage of this method 
of temperature compensation is that it cannot be used 


where bi-axial stresses may be present. 


(3) DETERMINATION OF LOAD 


The temperature-compensated gages discussed in the 
previous section may be used to measure the strains in a 
structure at the positions where they are attached. 
These strain measurements can be converted into 
corresponding values of stress by multiplying them 
by the Young’s modulus of the structural material. 


The stresses thus obtained are ordinarily present in 
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the structure as the result of both the externally applied 
load and the internal stresses resulting from temper- 
ature gradients in the structure. Since the stresses 
resulting from the latter cause are self-equilibrating, 
they produce no net axial load or moment in the struc- 
ture when their effect is integrated over the cross- 
section. Therefore, when the effect of the total stress 
is integrated over the cross-section the result will be 
the same as though only the stresses due to the exter- 
nally applied load had been integrated. By properly 
combining the outputs of several suitably located 
gages on a structural member, it is possible to carry out 
this integration and thus obtain the value of the ex- 
ternal load, irrespective of the stresses resulting from 
temperature gradients. 


(3.1) Axial Load, Constant Thickness 


When an axial load P is applied to a plate of uniform 
thickness, Fig. la, with temperature gradients across 
its width, the stress consists of a uniform stress dis- 
tribution, Fig. 1b, and a self-equilibrating thermal 
stress distribution, Fig. lc, which combine to give the 
stress distribution in the plate, Fig. 1d. 

The load carried by the plate is given by 


pat{ : o dx 
—h/2 


where o is the stress shown in Fig. 1d and x is the dis- 
tance from the mid-width of the plate. It is shown 
in the Appendix (5.1) that if the curve in Fig. 1d can be 
adequately approximated by a polynomial of third 
degree, the load can be determined with two gages 
located at A and B as shown in Fig. le. The load is 
given by 


P = ht (1/2(04 + op)] (1) 


where o, and og are the stresses at A and B, respec- 
tively. It can be shown by methods similar to those 
given in the appendix, that if the curve in Fig. 1d can 
be adequately approximated by a polynomial of fifth 
degree, the load can be determined with three gages 
located at C, D, and E as shown in Fig. lf. The load 
is given by 


5 | 
P= i] (oc + or) + >| (2) 
IS 9 


When the curve in Fig. ld is so complicated as to re- 
quire a polynomial of seventh degree to adequately 
approximate it, the load can be determined with four 
gages located at F, G, H, and J as shown in Fig. lg. 
The load is given by 


P = ht [0.174(or + o7) + 0.326(0, + o4)] (3) 


The coefficients in Eqs. (1), (2), and (3) are the same as 
those formulas for numerical integra- 


tion.! 


for Gauss’s 
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3.2) Bending Load, Constant Thickness 


When a bending moment ./ is applied to a plate of 
uniform thickness, Fig. 2a, with temperature gradients 
across its width, the stress consists of a linearly varying 
stress distribution, Fig. 2b, and a self-equilibrating 
thermal stress distribution, Fig. 2c, which combine to 
give the stress distribution in the plate, Fig. 2d. The 
moment carried by the plate is given by 


M=1 { ox dx 


where o is the stress shown in Fig. 2d. It is shown in 


the Appendix (5.2) that if the curve in Fig. 2d can be 
adequately approximated by a polynomial of third 
degree the moment can be determined with two gages 
located at J and A as shown in Fig. 2e. The moment 


is given by 
M = 0.1076 ht (oy —_ OK) (4) 


For determining the bending moment, three gages 
have no advantage over two because the third gage 
falls at the neutral axis and gives no additional informa- 
tion. Where the curve in Fig. 2d is so complicated 
as to require a polynomial of seventh degree to ade- 
quately approximate it, the moment can be determined 
with four gages, located at L, 7, N, and O as shown 


in Fig. 2f. The moment is given by 


M = h?t [0.05365(0, — oo) + 0.06445 (oy — on)] (5 


(3.3) Axial Load, Variable Thickness 


When an axial load P is applied to a member of vari- 
able thickness ¢, Fig. 3a, with temperature gradients 
across its depth, the stress consists of a uniform stress 
distribution, Fig. 3b, and a self-equilibrating thermal 
stress distribution, Fig. 3c, which combine to give the 
stress distribution in the member, Fig. 3d. The load 
carried by the member is given by 

wh /2 
P= | ot dx 


where o is the stress, Fig. 3d, and f¢ is the thickness of 


the member, Fig. 3a. It is shown in the appendix 
that if the curve in Fig. 3d can be adequately approxi- 
mated by a polynomial of third degree, the load can be 
determined with two gages located at P and Q or R 
and S as shown in Figs. 3e and 3f, respectively, for 
symmetric and nonsymmetric cross-sections. For sym- 
metric cross-sections, Fig. 3e, the gages are placed a 
distance a from the top and bottom of the cross-section, 


where 


c -| kx? ar| 


x = distance from mid-depth 
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Fic. 1. Axial load measurement in plate of constant thick- 


ness: (la) member being loaded; (1b) stress due to external load; 
(lc) internal stress due to temperature gradients; (1d) total 
stress on member; (le), (1f), and (1g) gage locations with 2, 3, 
and 4 gages, respectively 
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Fic. 3. Axial load measurement in a member of variable 
thickness: (3a) member being loaded; (3b) stress due to ex- 
ternal load; (3c) stress due to temperature gradients; (3d) total 
stress on member; (3e) location of two gages for ¢ symmetrical 
about center; (3f) location of two gages for ¢ not symmetrical 
[sec Eqs. (6)—-(12 
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Fic. 2. Bending moment measurement in plate of constant 


(2a) member being loaded; (2b) stress due to ex- 
internal stress due to temperature gra- 
gage locations 


thickness: 
ternal moment; (2c) 
dients; (2d) total stress on member; (2e) and (2f 
with 2 and 4 gages, respectively 
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Fic. 4. Bending moment .4 measurement in member of 


variable thickness (symmetric about mid-width (4a) member 

being loaded; (4b) stress due to external moment; (4c) internal 

stress due to temperature gradients; (4d) total stress on member; 
te) location of two gages [see Eqs. (13)-(15)] 
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The load is given by For nonsymmetric cross-sections, Fig. 3f, the gages 
P = (A/2) (op + ao) (R are placed a distance 
b h ' AD — BC — V A*D? — 3B°C? — 6ABCD + 4B*D + 4AC' 
= : (9 
2 (B> — AC) 
from the top and a distance 
V A*D? — 3B°C? — 6ABCD + 4B*D + 4AC 


h we Re = 
4 + 


») 


C= 


from the bottom of the cross-section, where A and C are given by Eq. 


i. a | ii ax | 


2(B? 


10 


— AC) 


(7), and 


(11 


wh /2 
D= | . Ix dv} 


The load for nonsymmetric cross-sections is given by 


A*D — 3ABC + 2B° + A V A*?D? — 3B°C? — 6ABCD + 4B°D + 4AC 
2 V A*?D? — 3B°C? — 6ABCD + 4B°D + 4A4C 


P= 


oR 1 


—~A2D + 3ABC — 2B* + A VA2D? — 3B°C? — 6ABCD + 4B°D + 4AC: 


2 V A*D? — 3B°C? — 6ABCD + 


(3.4) Bending Load, Variable Symmetric Thickness 


When a bending moment J/ is applied to a member 
of variable symmetric thickness, ¢, Fig. 4a, with tem- 
perature gradients across its depth the stress consists 
of a linearly varying stress distribution, Fig. 4b, and 
a self-equilibrating thermal stress distribution, Fig. 
4c, which combine to give the stress distribution in the 
member, Fig. 4d. The moment carried by the member 
is given by 


» 


M -| : otx dx 
h/2 


where o is the stress, Fig. 4d, ¢ is the thickness of the 
member, Fig. 4a, and x is the distance from the mid- 
depth. It can be shown that if the curve in Fig. 4d 
can be adequately approximated by a polynomial of 
third degree the moment can be determined with two 
gages located at 7 and U as shown in Fig. 4e, where 


d = (h/2) —- VE/C (13) 
C is given by Eq. (7), and 
eh /2 
E = | tx! dx (14) 
h/2 
The moment is given by 
M = (C/2) V(C/E) (or — ov) (15) 


where o7 and oy are the stresses at the upper and lower 
gage locations, respectively. 


1B°D + 4AC 


(3.5) Shear Load, Constant Thickness 


When a shear force S is carried by a plate of uniform 
thickness, Fig. 5a, with temperature gradients across 
its width, the stress consists of stresses resulting from 
the applied loads, Fig. 5b, and a self-equilibrating ther- 
mal stress distribution, Fig. 5c, which combine to give 
the stress distribution in the plate, Fig. 5d. The shear 
force carried by the plate is given by 

eh /2 
S=t | 7 dx 


where 7 is the stress shown in Fig. 5d. It can be shown 
that if the stress distribution in Fig. 5d can be ade- 
quately approximated by a polynomial of third degree, 
the shear force can be determined with two gages |’ 
and IV located as shown in Fig. 5e. The shear force is 


given by 
S = (1 2)th (Ty T Tw) 16 


When the stress distribution requires a polynomial ol 
fifth degree to approximate it, the shear force can be 
determined with three gages, Y, Y, and Z located as 


shown in Fig. 5f. The shear force is given by 


S = (th 


1S) (OTx — NTz — Sty) li 


(3.6) Shear Load, Gradually Variable Thickness 


When a shear force S is applied to a plate of gradu- 


ally varying thickness with temperature gradients 


across its width, the gage locations can be the same as 
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those given in Fig. 5, with a modification of Eqs. (16) 
and (17) to take account of the fact that shear flow, i.e., 
thickness times shear stress, is being integrated. This 
gives for two gages located at V’ and JW as shown in 
Fig. 5e, 


ba a (1 2)h (tyr, T twTw) (1S 


and for three gages located at XY, }, and Z, as shown in 


Fig. ot, 


S = [(1/18)h] (S5txrx + 5tz rz + Sty ry (19 


(4) GENERAL 


Formulas similar to those presented, but based on 
fitting curves of higher degree, can be derived. How- 
ever, it is believed that in most cases where additional 
points are considered warranted, it would be better to 
first divide the structure into several parts and then 
use the formulas given for each part. 

When the range in temperature is sufficient to cause 
appreciable differences in Young’s modulus at the vari- 
ous positions in the structure, the present method 
should yield correct results if the correct value for 
Young's modulus at each gage location is used; however, 
if an average value of Young’s modulus is used, some 


error may result. 


(5) APPENDIX 


The equations presented in the text were developed 


as follows: 


5.1) Axial Load, Constant Thickness 


Consider the stress in Fig. 1d to be expressed as the 


polynomial 


to] a+ bx + cx? + dx Al 


where x 1s the distance from the neutral axis, Fig. 6, and 
hc, and d are constants. Let s; and so be the values 
of x at the locations of two strain gages, and o; and on», 
respectively, be the stresses at these locations. De 
note by A, and A, two constants to be determined. We 


will now find the values of A,, A», s;, and s» for which 
K,o,t+ Keo = t o dx A2 
is identically true. We know that 
o=—-ar bs, tT CS, - ds, (A3 


and 


a. = a + bs, + c5:? + ds» (A4) 


$y carrying out the indicated integration, we find that 


i o dx = ah 


e N/es 


ch‘ 12 


Substituting Eqs. (A3), (A4), and (A5) into Eq. 


gives 
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Fic. 5. Shear load measurement in plate of constant thick 
ness 5a) plate carrying shear load; (5b) stress due to external 
load; (5c) internal stress due to temperature gradients; (5d) 
total stress on plate (He ind (df gage locations with 2 and 3 
gages, respec tive ly 

A + A h + D(A ‘ -F Kos T 

c(Ays,? + Kose? h?/12 
d(Kys;3 + Kos 0 (A6) 
Eq. (AG) will be identically satisfied, regardless of the 
values of a, 5, c, and d, if 
A, +A h 0: | 
As; T K .) ): 7 
(A7) 
Ays;*7 + Kos h*/12 " , 
Kis; T K .) 0 
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Fic. 6. 
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Coordinates used in developing formulas. 
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Solving Eqs. (A7) simultaneously gives S$; = 0.3872h; so = —0.3872h; K, = 0.1076h?; 
: nae ; K. = —0.1076h2 (A13) 
S$: = 0.2887h, ss = —0.2887h, Ky = Ko = 0.5h (AS8) 
These results are the basis for Fig. 2e and Eq. (4). 
These results are shown in Fig. le and Eq. (1). 
5.3) Axial Load, Variable Thickness 


(5.2) Bending Load, Constant Thickness , : ; ee 
Consider the stress in Fig. 3d to be expressed by the 


Consider the stress in Fig. 2d to be expressed by the polynomial in Eq. (A1), let s; and s» be the values of x 
polynomial in Eq. (Al). Let s; and s, be the values of at the locations of two strain gages and let o; and 0» 
x at the locations of two strain gages and let o; and o2» be the corresponding stresses. Denote by A, and A, 
be the corresponding stresses. Denote by A; and AK»  twoconstants to be determined such that 
two constants to be determined such that °h/2 
4/2 K, 0, + Ke o2 = | o tdx (A14) 
Kyoit + Kooot = if oxdx (AQ) silt 

—h 


h/2 Using Eq. (Al), 


Now, n/2 *h/2 h/2 
[ atdx = a| tdx + b [ xtdx + 
. bh? dh*® J —k/2 Lye 


h/2 2 h/2 « h/2 
f oxax = = a (A10) h/2 h/2 
—h/2 12 SO c [ xtdx +d x°tdx (A15) 
Substituting Eqs. (A3), (A4), and (A10) into Eq. (A9) 


gives Using the notation in Eqs. (7) and (11), and substi- 
tuting Eqs. (A3), (A4), and (A15) into Eq. (A14) gives, 
a(K,; + Ky) + b(Ays; + Koso — h3/12) + 


(Kit K,—A Ni. + Kn ~ D4 
c(Kys,2 + Kose?) + AA 2 ) + DO( Aisi + Kose re 


c(K,s;? + Koso” aa C) + 


d(Kisy + R252 i h, 80) = 0 (All) d(K,s;3 fe Koss —_ D) = O (A 16) 
Eq. (A11) is identically satisfied, regardless of the Eq. (A16) is identically satisfied, regardless of the 
values of a, b, c, and d, if values of a, b, c, or d, if 
K, + Kz = 0; K, + Ke = A; 
Kis, + Kos» — h?/12 = 0; Kis; + Kos» = B: 
Kis," + Kes.? = 0; (A12) Kis; + Kos2? = C; (A17) 
and | and 
K,s,3 oa Koso* = h® SO = () K,s,3 + Koss? -_ D 
Solving Eqs. (Al2) simultaneously gives The solution of Eqs. (A17) gives 
—AD + BC + V A2D? — 3B°C? — GABCD + 4B*D + 4AC3 cain 
5; = d 
2(B? — AC) 
~AD + BC — V A2D? — 3B°C? — 6ABCD + 4B*D + 4AC3 (A19) 
Ss = AL 
2(B? — AC) 
K A?D — 3ABC + 2B°+ A V A2D? — 3B°C? — 6ABCD + 4B°D + 4AC3 (A20) 
: > ] AS 
2 V A*?D? — 3B°C? — 6ABCD + 4B*D + 4AC? 
K —A°D+3ABC — 2B*+ A V A2D? — 3B°C? — 6ABCD + 4B°D + 44C3 (421) 
‘ 2 V A2D? — 3B°C? — 6ABCD + 4B°D + 4AC! i 
These equations are the basis for Eqs. (9), (10), and (12). For the symmetric case, B = D = 0, and Eqs. 


(6) and (8) are obtained. 
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Some Aerodynamic Effects of Streamwise 
Gaps in Low Aspect Ratio Lifting Surtaces 
at Supersonic Speeds’ 


Z. O. BLEVISS? ano R. A. STRUBLE}? 
Douglas Aircraft Company, Inc., Santa Monica Division 


ABSTRACT 


Low aspect ratio control surfaces may be used on supersonic 
uirplanes and guided missiles. Their deflections will generally 
‘reate gaps which are approximately in the flight direction (e.g., 
ill-movable tails and variable-incidence wings). The aerody- 
namic effects of such gaps have nearly always been neglected in 
theoretical investigations These gap effects are properly 
neglected for high aspect ratio control surfaces but, since these 
effects increase with decreasing aspect ratio, they may be im- 
portant for low aspect ratio control surfaces. 

The purpose of this paper is to obtain some understanding of 
the effects of streamwise gaps in low aspect ratio lifting surfaces at 
supersonic speeds by investigating a simplified configuration 
whose gap effects are expected to be typical. 
consists of a flat lifting half-delta wing which is adjacent and 
normal to a vertical wall such that the wing side edge is parallel 
The usual assumptions 


The configuration 


to the wall and forms a streamwise gap. 
of linearized supersonic wing theory are made and the following 
wing characteristics are studied as functions of gap size, Mach 
Number, and leading edge sweep: spanwise lift distribution, lift, 
center of pressure, drag due to lift, and vorticity at the trailing 
edge. In addition, in connection with wing-tail interference, the 
strengths and initial positions of the rolled-up trailing vortices are 
An illustrative example of the gap effects on wing- 
Viscous effects in the gap, 


investigated 
tail interference is also presented. 
which are important for very small gaps, are discussed qualita- 


tively 
NOTATION 
b = wing span 
( = wing chord 
Cp = coefficient of drag due to lift based on wing area 
C; = coefficient of thrust force due to leading edge suc- 
tion based on wing area 
C; = lift coefficient 
C = coefficient of rolling moment about wing side edge 
based on wing area and span 
3 = coefficient of tail pitching moment about apex of 
tail based on wing area and chord (see Illustra 
tive Example, Section 5) 
tan w Bb : 7 
d = = = subsonic leading edge parameter 


tan yu c 


Based on a paper presented at the Aerodynamics Session, 
Twenty-First Annual Meeting, IAS, New York, January 26-29, 
1953. Revised and received March 12, 1954. 

*This paper is an abbreviated version of “Some Effects of 
Streamwise Gaps on the Aerodynamic Characteristics of Low 
Aspect Ratio Lifting Surfaces at Supersonic Speeds,’”’ by Z. O 
Bleviss and R. A. Struble, Douglas Report No. SM-14627, April 
1, 1953, to which reference should be made for detailed discussion 
and derivations. 

+ Aerodynamics Engineer. 


t Aerodynamicist. Now at Illinois Institute of Technology. 


E(k’, v) = incomplete elliptic integral of second kind 
/ =F, 
di 
J 0 \ i ~§ 
E = E(k’) = E(k’, 1 complete elliptic integral of 
second kind 
I ; 
F = supersonic leading edge parameter 
€ 
F(k’,v) = incomplete elliptic integral of first kind = 
dt 
0 Vil — kyl -—- 
k = = gap parameter 
h 7 ) 
k’ = V1 — k? 
k, = effectively-infinite gap 
K = F(k’,1) = complete elliptic integral of first kind 
Kx = spanwise lift distribution 
M = free-stream Mach Number 
V 21 — d?) 
n = = 
Vili¢V1-@ 
q = !/.,)U? = free-stream dynamic pressure 
= gap width 
l = free-stream velocity 
6% 2 = right-hand wing coordinate system (Fig. 2 
Xe p = spanwise coordinate of center of pressure 
Xen = initial position of rolled-up trailing vortex 
Xo = position of zero vorticity at wing trailing edge 
£ = x/b 
a = angle of attack 
8 =VM-!1 
(x = vorticity at wing trailing edge 


/ 

I = strength of rolled-up trailing vortex 
iy = T/Ub 

m = Mach angle 


l ' x+s\? 
y = aj 

kN b+s 
p = free-stream density 


¢ = perturbation velocity potential 
wing apex angle 


INTRODUCTION 


i hee DEFLECTION OF control surfaces on airplanes 
and guided missiles will generally create gaps. The 
“‘leakage’’ of air through these gaps will tend to equalize 
the pressures across the lifting surfaces and will, thereby, 
affect the aerodynamic characteristics of the lifting sur- 
faces and the characteristics of their trailing vortex 


systems. 
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Two possible low aspect ratio wing, body, and control 
surface combinations for the delta plan form. 


Fic. 1. 


All existing theoretical results for deflected control 
surfaces have been obtained for the case of zero ‘‘leak- 
age’’ (the no-gap case). The purpose of this paper is to 
investigate the effects of gaps, these effects being con- 
sidered as corrections which are to be applied to the 
no-gap results. 

A control deflection is usually a rotation about some 
hinge line. In general, this rotation will create gaps 
which are approximately normal to the hinge line (e.g., 
the gaps formed at the tips of a deflected conventional 
subsonic aileron). In this paper, only the case where 
the hinge lines are approximately normal to the free- 
stream direction (gaps approximately streamwise) will 
be studied. This is the case which occurs most fre- 
quently in practice, especially when the gap effects may 
be large. 

In the immediate neighborhood of a gap the correc- 
tions to the pressures and velocities will generally be 
large. These corrections will decrease with distance 
from the gap. The magnitudes of the corrections to the 
integrated quantities such as forces and moments will be 
functions of the geometry of the gap, the geometry of 
the lifting surfaces affected by the gap, and the Mach 
Number. 

For high aspect ratio trailing-edge-type control sur- 
faces the gap effects on the control and the wing char- 
acteristics and on the trailing vortex system will gen- 
erally be small at all speeds. On the other hand, the 
gap effects on the control characteristics and on the 
trailing vortex system may be large for low aspect ratio 
control surfaces. If, in addition, the adjoining wing has 
a low aspect ratio, the gap effects on the wing charac- 
teristics may also be large. 

The conventional subsonic airplane has high aspect 
ratio trailing-edge control surfaces and hence the gap 
effects are small. Furthermore, since the subsonic air- 
plane configurations are relatively standardized the gap 
effects have long since been accounted for by experiment 


and experience. It should be remarked that gap effects 
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will also be small for the high aspect ratio swept wings 
and control surfaces. 

Supersonic airplane and guided missile configurations 
are far from standardized. Low aspect ratio wings, 
tails, and control surfaces may often be used. For ex- 
ample, low aspect ratio wing and control surface com- 
binations may be achieved by using all-movable wings, 
the wings and control surfaces being one and the same. 
For such cases gap effects may be large over most of the 
speed range. Fig. | illustrates two possible low aspect 
ratio wing and control surface combinations for the 
delta plan form. Fig. 1(b) illustrates the all-movable 
wing and the type of gap which is formed between the 
body and the deflected wing. 

Many supersonic guided aircraft will spend some 
fraction of their flight time at subsonic speeds. The low 
aspect ratio wing and control surface combination may 
show significant gap effects at these speeds. This gap 
problem is very difficult even for the simplified con- 
figuration considered in this paper and no attempt has 
yet been made to study it. 

In the analysis of a gap formed by the rotation of a 
planar surface with respect to another surface, serious 
difficulties arise from the fact that the boundary con- 
ditions must be considered nonplanar even when both 
surfaces are planar (e.g., the conventional aileron). If 
one of the surfaces is nonplanar (e.g., a missile with an 
all-movable wing) the problem is still more involved. 
Such nonplanar problems are very troublesome even 
when treated by linear theory. 

The gap problem is further complicated by the fact 
that viscosity is important when the gap is very small. 
For example, in the wing-body gap problem illustrated in 
Fig. 1(b), viscous effects will be important in the entire 
gap when the wing deflection is very small and, for large 
deflections, viscous effects will still be important near 
the hinge line, where the gap is always small. 

It is apparent from the foregoing discussion that the 
nonplanar-gap problem is not soluble analytically even 
without considering viscosity. Numerical solutions 
would be tedious and would not readily divulge the 
basic nature of gap effects. With this in mind, a very 
simplified model configuration which retains the es- 
problem—.e., a 


sential characteristics of the 


streamwise gap in a low-aspect-ratio lifting surface—is 


This configuration is shown in Fig. 2. 


gap 


investigated. 
It consists of a flat half-delta wing of zero thickness 
adjacent to a wall in such a way that the wing side edge 
is parallel to the wall and forms a streamwise gap. 
This configuration was chosen to simulate roughly the 
all-movable-wing control surface. By reflection in the 
wall this configuration may also be interpreted as a 
delta wing with a central streamwise gap (see the dash- 
line plan form in Fig. 2). The important simplification 
which has been achieved is the change in the type of gap 
so that the boundary conditions may be treated as 


planar. It should be noted that such gaps do actually 
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STREAMWISE GAPS IN LOW 


exist when control surfaces are undeflected but they will 
generally be small in this condition. 

The usual assumptions of linearized supersonic wing 
theory are made and the following wing characteristics 
are studied as functions of gap size, Mach Number, 
and leading edge sweep: spanwise lift distribution, lift, 
center of pressure, drag due to lift, and vorticity at the 
trailing edge. In addition, in connection with wing-tail 
interference, the strengths and initial positions of the 
rolled-up trailing vortices are studied assuming these to 
be given by the areas and centroids, respectively, of 
those regions of the vorticity distributions at the trailing 
edge where the vorticity has the same sign. An illustra- 
tive example of the gap effects on wing-tail interference 
is presented. 

For very small gaps these results must be corrected 
for viscous effects in the gap. These corrections are 
discussed qualitatively in a later section. 


THEORY 


The simplified gap configuration shown in Fig. 2 has 
been described in the Introduction. It was pointed out 
there that assumptions are made which reduce this 
problem to a problem in the usual linearized theory of 
planar supersonic wings. Most of the wing solutions 
used in carrying out this study are either given ex- 
plicitly in the literature or may be easily obtained from 
other related solutions without any new basic deriva- 
tions. These solutions will be used as needed with a 
minimum of explanation about their derivations. 

The quantities which are studied were listed in the 
Introduction. These quantities are studied for the com- 
plete ranges of gap sizes and leading edge conditions. 
The gap size is characterized by the gap parameter 
k = s/s + b, where s is the gap dimension and 6 is the 
span of the half-delta wing. This parameter varies 
from zero to unity as the gap dimension varies from zero 
to infinity. The ratio of the tangent of the wing apex 
angle to the tangent of the Mach angle will be denoted 
by d and its reciprocal will be denoted by f. It is con- 
venient to use d as the leading-edge parameter for the 
case of subsonic leading edges and f for the case of super- 
sonic leading edges, then both parameters vary between 
zero and unity. 

For general values of the leading edge parameters 
analytic solutions are not available or easily obtained 
for the complete range of the gap parameter. However, 
for two limiting values of the leading edge parameter 
such solutions are available. These limiting values are 
f = 0 and d = 0. Formally, f = 0 corresponds to 
either J = o (zero Mach angle) and d finite or b = 
and M finite. For this case there can obviously be no 
gap effects on quantities obtained by integration over 
the wing. For this reason this case does not prove to be 
useful for estimating the variations of integrated quan- 
tities with gap size for other values of the leading edge 


parameter. d = 0 corresponds to the well-known 
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Fic. 2. Simplified gap problem 
slender-body case, i.e., either JJ = 1 and 6 finite or 
b = O (for a given k) and 1 M< o. This case is 


particularly useful for estimating the solutions for in 
termediate values of d and f because it is the only case 
exhibiting variations with gap size for which an exact 
analytic solution is known for the complete range of k. 


(1) Lir1 


Slender-Body Solution (d = 0 


It is well known that for lift problems the slender 
body solution may be obtained by solving the two 
The 


potential problem in the cross-flow plane consists here 


dimensional incompressible cross-flow problem. 


of a flow at infinity of magnitude al’ directed normal to 
two flat plates (each of width 7) which are separated by 
a gap of width 2s. The solution to this problem was ob- 
tained long ago in connection with the problem of 
minimum induced drag for a wing with a gap in an in- 
compressible potential flow.’ From the potential 
on the top side of the wing at the trailing edge and the 
relationship between the spanwise lift distribution and 


this potential 
Ux) = g(x,O0+,¢ 
l 
the spanwise lift distribution is 


la l iD 
z E(k’, v) — — F(k’, »| 2a) 
l1—k K 


For zero gap this reduces to 


I(x) = 


= tagb V1— 727 (2b) 


1(x) 
h 0 
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Fic. 3. Spanwise lift distributions for various gap sizes -slender- 


body theory (d = 0). 


and for infinite gap 


i(x)| = 


k=1 


fagob V (1 — &) (2c) 


where & = x/b. 

The spanwise lift distributions for several gap sizes 
are shown in Fig. 3. The large effects of extremely small 
gaps should be noted (compare the k = 0 and k = 10-6 
curves). 

The slope of the lift curve is 


: 27 tan w , E 
Cia = =i"? — 2 = (3a) 
(1 — k)? K 
(C1.)p=0 = 2a tan w (3b) 
(Cr.)pa1 = w tan w (3c) 
To bring out the gap effect, the lifts with and without 
a gap are compared for a fixed angle of attack. This 
comparison is given by the lift ratio 
Ci (1 + k? — 2E/K) 
> = - (4a) 
(Cra)k=0 (1 — k)? 
which, for infinite gap, reduces to 
(Cz) = 1 
; a (4b) 


(Cradk 0 2 
The lift ratio is plotted in Fig. 4 as the solid d = 0 
Here the large effects of small gaps are strik- 
It is seen that the lift ratio decreases 


curve. 
ingly exhibited. 






































2 4 6 8 . 1.0 


i Oo 4 ok 


Fic. 4. Lift ratio versus gap size for various leading-edge condi- 


tions 
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with increasing gap size and attains a minimum value of 
one-half for infinite gap. 

The curve for the case of the leading edge perpendicu- 
lar to the flow, f = 0, is alsoshown. The one-parameter 
family of curves for intermediate values of the leading 
edge parameters is bracketed by these two known 
The only results which are available in the 
intermediate curves are for the 
These results 


curves. 
literature 
limiting values of zero and infinite gaps. 
will be necessary for the construction of the complete 


for these 


curves and they will now be listed. 


Subsonic Leading Edges (0 < d < 1) 


Zero Gap (k = 0O).—This case corresponds to the 
narrow delta wing, for which the results are well known. 


The spanwise lift distribution‘ is 


lagb 
I(x) me l — ~° a 
(Vv 1 — d*) 
x , 
iix)) = aqgb V1 — # db 
d=1 7 
and the lift curve slope’ is 
; 2rd 
B(Cia)e=0 = (6) 
k(V 1 — d?’) 


This case corresponds to the 
The spanwise lift distribution 


Infinite Gap (k = 1). 
narrow half-delta wing. 
was obtained, by integration, from the pressure distribu- 


tion.® 78 The result is 


4V 2agbV #1 — 2) 


I(x) = (7a) 
E(n)\¥1 + V1 - @&? 
8V 2 
l(x)| = agbV #(1 — &) (7b) 
d=1 w 
The lift ratio for infinite gap” ° is 
(5 ee E(V 1 — d?) 
= (3) 


(Cra)e=0 V2E(n) 41 + V1—d 


Supersonic Leading Edges (0 < f < 1) 


Zero Gap (k = O).--The spanwise lift distribution? is 


tagbf f? r 
Ke) = | a + %) cos! - rr 
rV 1 — f f(l1 + x) 
f7— . 
(1 — 2%) cos™! O<Z< f (9a) 
f(l — #  . 
fagbf(1 — &) 
I(x) = i.e 3 (9b 
V1 —- f 


Since the gap affects only the region within the Mach 
cone from the apex, Eq. (9b) is correct for all gap sizes. 
The lift curve slope is 


B(Cra)e-0 = 4 (10) 
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STREAMWISE GAPS IN LOW 


Infinite Gap (k = 1 The spanwise lift distribution 


was obtained from the pressure distribution,® * * with 
the result 
Sagb 
(x) = TVaHf— 2) 4 
aV1l +f ; 
‘ic = 2) Mi +. £) — 2 
cos! : G=2<= Ff (hi) 
wvI-s fi 
The lift ratio for infinite gap*” is 
(Cie); | 
“=! = (12) 


(oer 0 Vi-cef 


Effectively-Infinite Gap k, 

The foregoing infinite-gap results are actually valid 
jor finite gap sizes for all leading edge conditions except 
i= 0. The range of finite gap sizes varies with the 
leading edge condition. These remarks follow from the 
fact that, as the gap size is decreased from infinity, the 
aerodynamic characteristics of the wing are unaffected 
until the wing detects the presence of the wall, i.e., until 
the wing lies in the disturbance field which is reflected 
from the wall (or the disturbance field of the reflected 
wing). The minimum gap size for which the infinite-gap 
results are valid will, for obvious reasons, be referred 
to as the effectively-infinite gap. This gap size is de- 


termined by 


(13) 


k, is plotted in Fig. 4 as the dash-dot line which passes 
obliquely through the family of curves. 

It remains now to complete this family of curves in 
the range 0 < k < k,—the range within which the aero- 
dynamic characteristics vary with gap size. In this 
range the exact solution is known only for d = 0 and 
f = Q. For other leading edge conditions it is, in 
principle, possible to obtain the exact solution by pro- 
gressively extending the solution from k, toward de- 
creasing k. However, such calculations become in- 
creasingly difficult as the number of reflections between 
the wing and the wall increases. The complete exact 
solution is not attempted here. Instead, an estimated 
solution is obtained by using the foregoing results for 
the exact end points at k = 0 and k = k, and for the 
d = 0 solution in conjunction with some existing in- 
formation concerning the applicability of slender-body 
theory to nonslender configurations. 


Estimated Solutions for 0 < k < k, 

A detailed explanation of the reasoning leading to the 
estimated solutions may be found in the aforementioned 
Douglas report. 
/ = 0 curve describes the effects of a slender gap, i.e., 
a gap with zero reduced aspect ratio, in conjunction 


with a slender half-delta wing. 


ASPECT 


Very briefly, the reasoning is: The 


For a general value of 
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the leading edge parameter the gap is effectively infinite 
for k = k, and the gap may be considered slender, inso 
far as the gap effects on spanwise lift distribution and 
related quantities are concerned, if there are more than 
two or three reflections of the Mach waves in the gap. 
It should be possible, by some reasonable and simple 
modification of the d = 0 curve, to account for the 
nonslender wings and obtain an estimate of the exact 
solution for other values of d (and f). Such an estimate 
is obtained here for any d (or f) different from zero by 
Oandk = k,, 
by a curve which is a simple affine transformation of the 

= QOandk I, 
the curve for any 


connecting the two known end points, & 


known d = O curve which connects & 
If g(k) is the curve ford = Oand h(k 
other d, then 


1 — htk,) k 
h(k) = 1 — po z ( ) O<k <b, 
1 — g(1) k, 


The estimated solution for several values of d and f is 
shown in Fig. 4 as dashed lines. 

The accuracy of this estimated solution is difficult to 
ascertain. For very small gaps (not including zero gap) 
the accuracy is probably poor because in this range of k 
the slopes of the curves are very large and, hence, even 
small errors in the shapes of the curves will result in 
large errors in the values of the lift ratio. However, this 
inaccuracy iS unimportant since in the same general 
range of & the exact solution is itself incorrect due to 
viscous effects. Outside of the range of k where viscous 
effects are important it is expected that the estimated 
solution represents a reasonable engineering approxima- 
tion to the gap effects. 

Some indication of the accuracy of this estimate may 
be obtained by extending the exact solution from k, for 
one reflection. One reflection is defined as that range of 
k in which only the original disturbances produced by 
the wing alone affect the wing after reflection from the 
wall. The extension for one reflection is easily carried 
out; the extension for more than one reflection would 
involve a prohibitive amount of numerical calculation. 

The results for one reflection are shown in Fig. 4 as 
the solid curves which extend below &, and are ter- 
minated by the short vertical lines. For supersonic 
leading edges the final results were obtained in closed 
form and are shown for the complete range of the lead- 
ing edge parameter. For subsonic leading edges the final 
integral was evaluated numerically for the one inter- 
mediate value of the leading-edge parameter d = 0.5. 

It is seen that in the range of & for one reflection the 
exact and estimated solutions are identical within the 
This fact together with the 
foregoing about the construction and 
accuracy of the estimated solution indicates that the 


accuracy of the graph. 


discussions 


estimated solution is probably a good approximation to 
the exact linear solution over most of the range of k, the 
region near k = 0, where the slopes of the curves are 
large, being excluded, as previously remarked. 
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Fic. 5. Spanwise coordinate of center of pressure versus gap size 


for various leading-edge conditions 


(2) SPANWISE COORDINATE OF CENTER OF PRESSURE 


Exact results are available or are readily obtained 
when the gap is zero or effectively infinite, k = O or 
k = k,, for all leading edge conditions (for all d and f), 
and also for all gap ratios, k, when the leading edge 
sweep is zero or 90 degrees (f = 0, d = 0). The exact 
solution has been extended for one reflection for the case 
of supersonic leading edges. The extension for one re- 
flection has not been carried out for subsonic leading 
edges. The formulas will not be presented. 

Using these results, the general plot of the spanwise 
coordinate of the center of pressure versus gap size for 
various values of the leading edge parameters may be 


constructed in the same way as for the lift case. This 
plot is shown in Fig. 5. The k = 0 end points are de- 
noted by the circles on the axis of ordinates. The 


uppermost end point is labeled and the other three end 
points correspond in an obvious way to the three non- 
zero values of f. 

As in the lift case, the exact and estimated solutions 
are identical for one reflection within the accuracy of 
the graph. Again, the estimated solution is expected to 
approximate well to the exact solution. 

It is of interest to note the properties of the solution 
for subsonic leading edges. All of the curves pass 
through the same end points atk = Oandk = 1. The 
curve k = k, lies along the line ¢,, = 0.5 and the d = 0 
curve lies below this line. From these facts it follows 
that the d = 0 curve is the lower bound for this family 
of curves and the d = | curve is the upper bound. It is 
seen that these bounding curves do not differ greatly. 

The maximum possible outward movement of the 
center of pressure occurs for subsonic leading edges 
(with k = k,) and is about 7.5 per cent of the span. 


(3) STREAMWISE COORDINATE OF CENTER OF PRESSURE 


The results may be summarized very briefly. When 
k = Oandk = k, the configuration is conical and the 
center of pressure must lie two-thirds of the wing chord 
length aft of the wing apex. As # is decreased from k,, 
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for a given value of d or f, the center of pressure moves 
first rearward then forward again as k approaches zero. 
The maximum rearward shift of the center of pressure 
is relatively small, being less than about 2 per cent of 
the wing chord length. It should be noted that even 
this small shift might be important in some _ hinge 


moment problems. 


(4) DraG DvuE To LiFt 
The coefficient of drag due to lift may be written 
Cp = Cla —™ Cr 14 
where Cr is the coefficient of the thrust force due to 
leading edge suction and is zero for supersonic leading 
edges. 

To bring out the gap effects, it is of interest here to 
compare the drags with and without a gap for the same 
lift. This comparison is obtained by forming the drag 
ratio which, in the absence of leading edge thrust, is 


(15) 


(€5/Cr* (Cras -() 
(Cp 'Cr?), 0 am 


The slender-body (d = 0) results for zero and infinite 
gaps are well known.” * '° These show that the drag 
ratio for infinite gap is the same with or without leading 
edge suction and is equal to 2. The complete slender- 
body results have been calculated for all gap sizes and 
the comparison of the drag ratios with and without 
leading edge suction is shown in Fig. 6. It is seen that, 
at least for the slender-body case, the gap effects are 
only slightly affected by leading edge suction. 

In the general plot of drag ratio versus gap size, some 
assumption about the leading edge thrust must be made 
for subsonic leading edges. In a real fluid, no leading 
edge thrust is to be expected for sharp leading edges. 
Presumably, then, varying amounts of leading-edge 
thrust can be obtained depending upon the cross-sec- 
tional shape of the leading edge. In any case, the as- 
sumption of no leading-edge thrust provides a con- 
servative representation of the gap effects. Fig. 6 is 
based upon this assumption. The family of curves is 
readily obtained from Fig. + using Eq. (15). 
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Fic. 6. 
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STREAMWISE GAPS IN LOW 


(5) Winc-Tatt INTERFERENCE 


Large effects of gaps on wing-tail interference can, in 
turn, result in large effects on the stability and control 
of an airplane or guided missile. Since wing-tail inter- 
ference is produced by the vortex system which trails 
behind a wing, it is important to investigate the gap 
effects on the trailing vortex system. 

In this section the gap effects on the half-delta wing's 
trailing vortex system are investigated. This investiga- 
tion is based on a simplified theory of the trailing vortex 
system suggested by the investigations of Lagerstrom 
ind Graham!! and Spreiter and Sacks.'* This simplified 
theory is based primarily upon the fact that the vortex 
sheet behind a low aspect ratio wing rolls up very 
rapidly if the wing lift is not near zero and the theory 
appears to be satisfactory for engineering purposes. 

The trailing vortex system behind the half-delta 
wing is obtained, according to the simplified theory, as 
follows. First, the vorticity distribution at the trailing 
edge is calculated by linear theory. Then, this vorticity 
distribution is replaced by one or more discrete vortices 
starting at the wing trailing edge, it being assumed that 
each region of constant sign is replaced by a single vor 
tex. The streneths and initial positions of these vortices 
are determined by the areas and centroids, respectively, 
of the regions of the vorticity distribution which the 
vortices replace. Finally, the motion of these vortices 
as they proceed downstream 1s determined by assuming 
that the flow in each cross section is two-dimensional 
and incompressible. 

Assuming a vertical wall between the wing and the 
tail, the motions of two types of vortex configurations 
must be considered in the present problem; two vortices 
for the no-gap case and four vortices for the gap case. 
The motions of such vortex configurations have been in- 
vestigated by Lagerstrom and Graham!! and, based on 
their results, it will be assumed here that the vortices lie 
in the free-stream direction. 

The wing-tail interference is finally obtained by de- 
termining an apparent tail twist from the two-dimen- 
sional cross-flow velocity field of the vortex system at 
the tail and then calculating by linear theory the forces 
and moments on the tail due to this twist. It is clear 
that this part of the wing-tail interference is a function 
of the position and geometry of the tail and is not an 
essential part of the gap effect. 


Vorticity Distribution at the Trailing Edge 


The vorticity distributions at the trailing edge can be 
obtained for those cases whose spanwise lift distributions 
were given in Section | by the following relationship: 


U dl(x) 


—_ j (16) 
2q dx 


v(x) = 


The formulas will not be given. Fig. 7 shows the slen- 


der-body results for several gap sizes. 


ASPECT 
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Vorticity distributions at the trailing edge for various 
gap sizes—slender-body theory (d = 0) 


Vortex Strengths 


From the basic assumptions it is easily shown that 
the circulation is zero for any closed path in the fluid 
which does not intersect the half-delta wing or its trail- 
ing vortex system. This means that the areas under 
the positive and negative portions of the vorticity dis- 
tribution at the trailing edge are equal. Then the two 
vortices which are assumed to start at the trailing 
edge have equal strengths and are of opposite sign. 
The vortex which is nearer to the wall will be referred to 
as the inboard vortex and will be denoted by the sub- 
script 1, the other vortex being outboard and denoted 
by the subscript 2. 

The vortex strength is related to the spanwise lift 


distribution by 


V 
I(x) (17) 
2q 


where x» is the coordinate of the point of zero vorticity, 


r=r.=-—T,= 


y, at the trailing edge. To bring out the gap effect, the 
vortex strengths with and without a gap were com- 
pared for the same angle of attack. This comparison 


was obtained by forming the vortex strength ratio 
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(a) (b) 


Illustrative example of gap effects on wing-tail inter- 
ference 


Fic. 8. 


where [ denotes the nondimensional vortex strength 
r/Ub. 

The formulas will not be given. 
vortex strength ratio versus gap size for various leading 
edge conditions is nearly identical with the lift ratio 
plot (Fig. +). 

It should be noted that, even for zero gap, an inboard 
vortex exists, theoretically, at the wing-wall juncture. 
Of course, as must be the case, it may be neglected since 
its reflected counterpart cancels it. 


The general plot of 


Vortex Positions 


The initial positions of the inboard and outboard 
vortices may be obtained by the following integrations: 


1 "x0 
S25, = r / xy(x) dx = 
1 0 


+ os 


| "i 
/ I(x) dx (18a) 
l(xo) 9 


b 
1 (x ) dx 


(18b) 
l(xo) % 


Kegs = HOF 

The behavior of the initial vortex positions as func- 
tions of gap size and leading-edge condition may be 
summarized as follows: For any given value of d or f, 
both vortices move outward as k increases from zero. 
This outward movement increases as f increases from 
zero to unity and then remains the same for all values of 
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d. The zero-gap position of the inboard vortex is at the 
wing-wall juncture for all leading edge conditions. The 
zero-gap position of the outboard vortex is at the half- 
span when f = 0, moves outward as / increases to unity, 
and then remains fixed at about 0.80 for all values of d. 
The maximum possible outward movement of both 
vortices as the gap is increased occurs for subsonic 
leading edges and is about 10 per cent of the span. 


Illustrative Example 


A numerical example which illustrates how the 
presence of a gap may affect the wing-tail interference 
has been carried out. 
A lifting half-delta wing adjacent to a vertical wall is 
considered with and without a gap. A flat-plate tail 
whose plan form is identical with that of the wing and 


whose geometrical angle of attack is zero is assumed to 


The example is shown in Fig. 8. 


be situated several chord lengths aft of the wing in the 
free-stream direction. This wing-tail configuration 
would correspond to that of a missile whose body and 
tail are at zero angle of attack and whose wing is at in- 
cidence. Since the vortices are assumed to trail in the 
free-stream direction they will lie in the plane of the tail 
for this example. 

Fig. S(a) shows the wing-tail configuration for the no- 
gap case. The vertical wall is assumed to extend from 
Fig. S(b) shows the wing-tail configuration 
Aft of the wing the vertical wall ex- 


wing to tail. 
for the gap case. 
tending between wing and tail is assumed to be in line 
with the wing side edge so that the tail lies directly aft 
of the wing. This was chosen so as to simulate, 
roughly, the wing-tail positions for a wing-body-tail 
configuration with a gap formed by wing incidence. 
The gap and no-gap cases are compared for the same 
This will indicate the effect of 
The gap size was chosen large 


wing angle of attack. 
neglecting the gap. 

enough so that viscous effects would probably be negli- 
gible for usual wing thicknesses. Assuming k = 0.05 and 
sonic leading edges, the following data were obtained 
from the formulas for the strengths and initial positions 


of the vortices: 


k= 0 k = 0.05 
m 4 rs r 
a= — = 1,27 =——=094 
Qa T Qa Qa 
Tv = a 
E.,, = — = 0.79 .,, = 0.09 
, { 
£.5. = 0.88 


From these vortices and their reflected counterparts, 
the apparent tail twist distributions a, and a’, (see Fig. 
S) were determined. The lifts and moments on the tail 
due to these twist distributions were calculated by linear 
theory. 

The lift, pitching moment about the tail apex, and 
rolling moment about the wall-tail juncture were calcu- 
lated for the tail. These results together with the re- 
sults for the wing and for the wing plus tail are summar- 
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STREAMWISE GAPS IN LOW 


ized in the table shown in Fig. 9. The sign conventions 
for the tail are the same as for the wing. The contribu- 
tions of the inboard and outboard vortices are given 
separately for the gap case. 

Taking the zero-gap case as the basis for comparison, 
it is seen that the gap decreases all of the forces and 
moments which are induced on the tail; the induced tail 
lift is decreased 55 per cent, the induced pitching 
moment about the tail apex is decreased 37 per cent, 
and the induced rolling moment about the tail side edge 
is decreased 34 per cent. In addition, the gap decreases 
the wing lift 26 per cent and decreases the rolling moment 
of the wing about its side edge 14 per cent. 

Some indication of the gap effects on the stability and 
control of a complete missile may be obtained by con- 
sidering the wing-plus-tail combination. For the com- 
bination the gap actually increases the lift 45 per cent 
and increases the rolling moment 40 per cent. The gap 
decreases the pitching moment of the combination about 
The 
decrease is due primarily to the changes in the wing and 
tail lifts and hence varies from about 55 per cent to 


any pitch axis lying between the wing and the tail. 


about 26 per cent as the position of the pitch axis varies 
from the wing center-of-pressure to the tail center-of- 
pressure. 

It should be emphasized that the large tail contribu- 
tions to the forces and moments on the wing-plus-tail 
combination are due to the fact that the tail lies 77 the 
trailing vortex system from the wing. Generally, the 
missile will have some angle of attack and the tail will be 
situated some distance from the trailing vortex system. 
Then, the induced forces and moments on the tail will 
be much smaller relative to the forces and moments on 
the wing and the lift and rolling moment on the com- 





















































bination will probably be decreased by the gap. The 
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Fic. 9. Numerical results for illustrative example (comparison of 


gap and no-gap cases) 
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Fic. 10. Qualitative effects of viscosity for small gaps. 


gap effects on the induced tail forces and moments will, 
however, be the same as indicated above. 

There are, of course, many other considerations in- 
volved in the study of stability and control but this 
simple example indicates the possibility of large gap 
effects. 


DISCUSSION OF VISCOUS EFFECTS 


The foregoing theory is subject to the usual assump- 
tion of linearized supersonic wing theory for steady 
flow. One of the basic assumptions is that the viscosity 
is zero. This assumption breaks down when the gap 
size is effectively reduced by viscosity, i.e., by the 
boundary layers on the wing side edge and on the wall. 
This means that the results must be corrected for 
viscous effects for some range of gap sizes near zero 
gap. In the following qualitative discussion these 
viscous effects are considered in more detail. 

Fig. 10(a) illustrates the effects of viscosity on the 


lift ratio. The region near zero gap is shown greatly 


enlarged. The solid curve represents a typical curve of 
Fig. 4. This sketch is based on the following qualitative 


reasoning: Consider the variation of the lift ratio for a 
wing of nonzero thickness and for a fixed leading edge 
condition as the gap size is decreased from some large 
value. The variation will be given by the nonviscous 
solution until a ‘‘gap size’’ is reached for which the 
viscosity effectively reduces the gap size and, hence, 
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effectively increases the lift. This ‘gap size’’ will be of 
the order of magnitude of some pertinent boundary 
layer thickness in the gap and is denoted in the sketch 
by k,,. 
effects will become progressively larger and will eventu- 
ally dominate and produce a Poiseuille-type flow in the 
gap. 
two-dimensional Poiseuille flow indicate that the slope of 
the viscous solution at zero gap will probably be zero. 
The viscous solution will be a function of Reynolds 
number, Mach Number, leading edge sweep, angle of 
attack, and the thickness ratio and profile shape of the 


As the gap size is decreased below k, the viscous 


Some simple considerations using the solution for 


wing side edge. 

The nonviscous solutions for all but one of the other 
aerodynamic quantities will be modified by viscosity in 
a similar manner. The exception is the inboard vortex 
strength. This is readily seen from the fact that the 
viscous and nonviscous solutions do not coincide for 
zero gap. Fig. 10(b) illustrates how viscosity modifies 
the nonviscous solution for the inboard vortex strength 
(the vortex strength ratio is based on the zero-gap end 
point for the nonviscous problem). 

It should be noted that the viscous solutions for the 
wall and for the reflected-wing problems will be dif- 
ferent. k, will be larger for the wall problem than for 
the reflected-wing problem. 

The viscous effects will probably have to be in- 
vestigated experimentally. 


CONCLUDING REMARKS 


Applicability of Present Results 


The present results are applicable to the residual gaps 
in undeflected control surfaces when k > k,. Although 
these results are not directly applicable to the more im- 
portant case of the rotated gap (the gap formed by the 
rotation of a control surface), they indicate qualita- 
tively the gap effects for the rotated gap and they can 
provide some quantitative conservative estimates of 
these effects. 

The leakage through a rotated gap will have the same 
qualitative effects on the spanwise lift distribution as 
the leakage through a translated gap (the gap which is 
studied in the present paper). Then the gap effects will 
be qualitatively similar for the two types of gaps for all 
of the quantities which were investigated with the 
probable exception of the streamwise coordinate of the 
center of pressure. 

The ‘‘size’”’ of a rotated gap is limited by the maxi- 
mum deflection of the control surface. For this type of 
gap there is no effectively-infinite gap size and no range 
of gap sizes in which the gap effects are constant. The 
nonviscous solution will probably have to be corrected 
for viscous effects for all gap sizes. 

A conservative estimate of the gap effects for the 
rotated gap may be obtained from the present results 
by using the ‘‘maximum linear dimension” of a rotated 
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gap as an equivalent translated gap. Such an estimate 
is admittedly crude but it does give some quantitative 
information about the effects of gap size, Mach Number, 
and leading edge sweep for the rotated gap. 


Recommendations for Future Research 


Much more research is necessary in order to provide 
satisfactory quantitative corrections to the no-gap re- 
sults for rotated gaps and for small translated gaps. 
Some recommendations for this research are given be- 
low. 

There is little hope for a de- 
includes 


Theoretical Research. 
theoretical treatment which 
Because of viscous effects, theoretical investi- 


tailed viscous 
effects. 
gations of nonviscous problems for the rotated gap can 
be justified only if the solutions are relatively easily 
The slender-body problem appears to be the 
The solution of this 


obtained. 
only one satisfying this condition. 
problem would be useful for estimating the nonviscous 
solution for the rotated gap for other leading edge con- 
ditions by comparing this solution with the slender-body 
solution for the translated gap. In this way, any ro- 
tated gap would be related to an equivalent translated 
gap and this equivalence would be assumed to apply to 
the estimated solution for the translated gap for other 
leading edge conditions. 

Experimental Research. The viscous solutions for the 
rotated and translated gaps will have to be determined 
by systematic experimental research. This research 
should include a thorough investigation of the gap 
effects on the trailing vortex system. 

When large control surface deflections occur, the 
present results may be significantly modified by non- 
linear and separation effects. These effects should be 
investigated. 
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Lite Expectancy of Aircraft Under Thermal 
Flight Conditions 


GEORGE GERARD? 
New York University 


SUMMAR\ 


This paper is concerned with certain new structural problems 
which result from flight of aircraft under variabie stress and 
elevated temperature conditions arising from aerodynamic heat 
ing. Since creep is a governing factor in determining life ex 
pectancy under elevated temperatures, a review was made of 
test data under cyclic stress conditions at elevated temperatures 
These data indicate that gust loadings under elevated tempera 
ture conditions may not materially affect the life expectancy of 
the aircraft. Furthermore, the total creep under cyclic loading 
conditions appears to result from the net time spent at a particu 
lar stress level. Based on these data, a cumulative creep hy 
pothesis is suggested for use in analysis of life expectancy under 


variable stress and temperature conditions 


INTRODUCTION 


7 ADVENT OF SUPERSONIC FLIGHT has and will con 
tinue to introduce many new structural problems 
associated with the effects of elevated temperatures re- 
sulting from aerodynamic heating. One significant 
feature of this problem is that the design of a structure 
which operates at elevated temperatures entails analysis 
of the time dependent deformations (creep) which occur 
under applied load. Ina survey of such problems, Hoff! 
has indicated that an entirely new philosophy of design 
may have to be adopted in which piloted aircraft may be 
considered as semiexpendable. Depending upon the 
service application of the aircraft, various criteria may 
be used to indicate the time after which deformation 
due to creep may become excessive, thus necessitating 
retirement of major components of the aircraft struc- 
ture. 

In the design of an aircraft for flight at a Mach Num- 
ber greater than two, the structure during its service 
life is exposed to different stress levels and different ele- 
vated temperatures as a result of maneuver and gust 
The problem of life ex- 
determining the 


loads imposed during flight. 
pectanecy therefore, is not one of 
deformations after a period of time for a structure 
which operates at one stress level and at one tempera- 
ture. It is the more complex one of determining the 
deformations as a function of time when the structure 


is subject to a distribution of different stress levels 
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operating at various elevated temperatures during the 
life of the structure. In effect, one may picture the 
ordinary velocity-acceleration diagram for an airplane 
to which a time axis has been added so that at a constant 
velocity the distributions of accelerations are given as 
functions of time. 

The prediction of life expectancy of an aircraft flying 
under elevated temperature conditions entails an inti 
mate knowledge of the applied stress and temperature 
distributions as well as a complete understanding of the 
behavior of materials under the complex loading con- 
ditions associated with the aircraft operation. From a 
physicist’s and metallurgist’s point of view, the phe 
nomenon of creep even under simple loading is rather 
incompletely understood. Thus, the engineer in facing 
the problem of predicting life expectancy of aircraft 
under elevated temperature conditions is forced to 


With this 


approach, hypotheses which may ignore certain features 


adopt a wholly phenomenological viewpoint. 


of creep of significance to the physicist or metallurgist, 
may be acceptable to the engineer provided the results 
obtained are commensurate in accuracy with the initial 
data used. 

In life expectancy analysis, the initial information of 
the problem consists of the basic creep data for the ma- 
terials under consideration and the applied stress and 
temperature levels which the aircraft will experience. 
The latter is completely statistical in nature whereas 
basic creep data also indicate statistical variations from 
a mean value. 

In combining these data to predict life expectancy, it 
is necessary to introduce both a failure criterion and:a 
life expectancy hypothesis. At first, the failure criterion 
may be purely arbitrary in nature and the wisdom of the 
selection ultimately will be proved or disproved by ex- 
perience. Among the failure criteria which may be 
selected are the following: stress rupture, creep buck- 
ling, initiation of third-stage creep, and specified defor- 
mations such as 0.5, 1, or 2 per cent strain. The cri- 
terion selected may vary depending upon how critical is 
the function of the structural component. 

In conjunction with a suitable failure criterion, it is 
necessary to use a life expectancy hypothesis which 
combines the initial data to predict the life expectancy. 
Since the initial data are statistical in nature, the results 
obtained by use of such hypotheses must be ex peri- 
mentally verified by use of statistical methods. 
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AIRCRAFT LOADING CONDITIONS 


As a result of various maneuvers, an airplane flying 
at a Mach Number greater than two will be subject to 
a distribution of stress levels and temperature levels 
during its service life. In the past, statistical data on a 
great variety of aircraft have been accumulated for 
establishing airload criteria. With temperature as an 
additional factor for flight at Mach Numbers greater 
than two, statistical flight load data should include the 
distribution of stress levels as a function of temperature 
level. Until such data are available, it is pertinent to 
consider qualitatively the creep conditions resulting 
from various types of maneuvers. 

In straight and level flight, calculations have indi- 
cated that the stress levels for lg load factor may be 
relatively low from a creep standpoint. However, the 
airplane may be in straight and level flight for a large 
proportion of the total time spent under elevated tem- 
perature operating conditions. Thus, the amount of 
creep resulting from this maneuver may be large. 

If the airplane flies at a relatively constant Mach 
Number, then the temperature conditions will be rela- 
tively constant. However, if the speed of the aircraft 
is increased still maintaining a straight and level flight 
path, then the temperature will increase although the 1g 
stress level may remain unchanged. Consequently, due 
to thermal stresses, the structure will be intermittently 
subjected to levels higher than those corresponding to 
lg. Consequently, the creep will be accelerated relative 
to that for constant velocity flight. 

With aircraft now in operation, it is well known that 
intermittent stresses due to gusts are superposed on the 
lg stress level. Although relatively little is known 
about gust conditions in the upper atmosphere, it is 
possible to obtain some estimate of the magnitude of the 
intermittent stresses by use of theory developed from 
relatively low altitude gust data. Such calculations in- 
dicate, for example, that in straight and level flight 
under the conditions of a 30 ft. per sec. gust, the stresses 
resulting from the incremental load factor may be ap- 
proximately only 50 per cent of the Ig stress level under 
thermal flight conditions 

In almost any other maneuver under elevated tem- 
perature conditions, high intermittent stress levels will 
be applied. These stresses result from accelerations im- 
posed by maneuvers and thermal stresses arising from 
velocity changes which induce variations in thermal con- 
ditions. Although such stresses may be applied for rela- 
tively short periods of time as compared to straight and 
level flight, the high stress and temperature levels may 
contribute significantly to creep of the structural com- 


ponents, 


PROPERTIES OF MATERIALS UNDER ELEVATED 
TEMPERATURE CONDITIONS 


The effects of elevated temperatures upon the 


physical properties of materials generally involve losses 
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in strength and also introduces time as an additional] 
parameter of the problem. It is the purpose here to 
briefly discuss the various effects caused by elevated 
temperature exposure and then consider in some detai 
those factors which are of concern in analysis of life ex. 
pectancy of an airframe. 

An obvious feature of short-time elevated tempera- 
ture exposure is the change in stress-strain characteris- 
tics of the material such as decreased yield and ultimate 
strengths and lower elastic modulus. For materials 
which are subject to age hardening, elevated tempera- 
ture exposure for short times can result in an increase in 
properties for a limited temperature range. Generally, 
however, elevated temperature exposure is associated 
with a deterioration in physical properties of the ma- 
terial. Under no-load conditions for long periods of time 
deterioration may become severe. 

For long-time exposure to elevated temperatures 
under load, creep of the materials becomes of im- 
portance. Because of the long interest in this subject, 
there is a considerable amount of literature on tension 
creep under constant load and constant temperature 
conditions. However, when one considers creep under 
compression loads, the problem of creep buckling of 
columns and plates assumes paramount importance in 
the case of aircraft design. Creep buckling of columns 
has received some attention in the past few years al- 
though extensive test data are still lacking. The statis- 
tical nature of the creep problem should be empha- 
sized since not only the basic creep properties indicate 
statistical variations, but initial imperfections are of 
considerable importance in creep buckling of columns. 
There is apparently no published literature on the im- 
portant problem of creep buckling of plates under com- 
pressive and shear loadings. 

As discussed previously, the airframe during its serv- 
ice life is subject to a distribution of different stress 
levels and various elevated temperatures. Conse- 
quently, in life expectancy analysis of components of the 
airframe in which we are concerned with the long time 
aspects of elevated temperature effects, creep data under 
varying stress and temperature conditions becomes of 
major importance. In the following section a review is 
presented of some of the literature on this subject. 


Alternating Stress, Constant Temperature Creep 


Lazan* has studied the effect of superposition of 
rapidly alternating stress upon direct tensile stress at 
elevated temperatures. The alloys investigated were 
the heat resistant alloys which have been developed for 
gas turbine applications and include: N-155, S-816, 19- 
9DL and Vitallium. Some of the conclusions of this in- 
vestigation are: 

(a) The test data show that at 1,500°F. the super- 
position of up to 67 per cent alternating stress on a given 
tensile mean stress does not generally affect the allow- 
able static mean stress for a given rupture life or dy- 
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LIFE 


namic creep by more than & per cent for the materials 
tested. 

(b) At 1,350°F., the presence of 67 per cent alternat- 
ing stress may reduce the allowable static mean stress 
by as much as 25 per cent. In general, the shorter the 
hours to failure, the greater the reduction caused by the 
superposition of alternating stress. 

(c) The addition of an alternating stress may in 
many cases have actually increased the allowable static 
mean stress at 1,500°F. for a given rupture time or 
creep. Increases of up to 8 per cent for 25 per cent al- 
ternating load and up to 3 per cent increase for 67 per 
cent alternating load were observed. 

A detailed study* of the alloy N-155 was reported in 
which alternating stresses were superposed on mean 
stresses at elevated temperatures in a number of dif- 
ferent types of testing machines. The test data for 
fracture in 50, 150, and 500 hours in a temperature range 
between 1,000° and 1,500°F. support Lazan’s conclu- 
sions given above. 

Some tests on 145-T aluminum alloy by Manjoine* 
indicated that the minimum creep rate or rupture time 
in dynamic creep tests may be higher or lower than that 
for a constant load test of the same mean stress. For 
the 14S-T alloy tested at 400°F., the minimum creep 
rate for an alternating load test of 20 cycles per sec. is 
higher for a mean stress above 25,000 lbs. per sq.in. and 
is lower for lower stresses. The rupture time is poorer 
for a mean stress above 22,000 Ibs. per sq.in. and is 
greater than that for a constant load test below this 
stress. 

The results summarized above have all been obtained 
by superposing a relatively rapidly alternating tension 
stress upon a mean tensile stress at elevated tempera- 
tures. As such, these results would have application to 
an airplane flying under gust conditions at elevated 
temperatures. Calculations of the magnitude of the al- 
ternating stress resulting from gust conditions indicate 
it is generally a small fraction of the mean stress. The 
results on the heat resistant alloys and 145-T aluminum 
alloy indicate that to a first approximation the effects of 
small alternating stresses superposed on a mean value 
may not appreciably effect the creep characteristics 
under the corresponding mean value of stress. There- 
fore, it may be tentatively concluded that the creep 
strength of an airplane under gust conditions may not 
be appreciably altered from the steady loading condi- 


tions. 


Cyclic Stress, Constant Temperature Creep 


Under general maneuvering conditions, the airplane 
may be exposed intermittently to high levels of stress 
under elevated temperature conditions. In such cases, 
it is appropriate to consider high or low stress pulses 
superposed at definite intervals upon a mean stress 
level. This type of loading is designated “cyclic load- 
ing”’ to distinguish it from the rapidly alternating type 
of stress considered in the previous section. 
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Yerkovich® investigated type 347 
stainless steel at 1,200° and 1,500°F., N-155 alloy at 
1,500° and 1,800°F., and H.S. 21 Type alloy at 1,500°F. 
In these tests various 


Guarnieri and 


under cyclic stress conditions. 
overstress to normal stress ratios were used in which the 
overstress was applied for 5, 10, and 25 per cent of the 
total test time. Various overload time patterns such as 
shown in Fig. 3, were used and it was concluded that for 
a given ratio of overload time to normal load time, the 
creep characteristics were not critically affected by 
various timing patterns of overload application. 

The test data are presented in a convenient form 
where normal stress is plotted as the ordinate and time 
as the abscissa. For a stated percentage of creep or 
fracture, the data for 5, 10 and 25 per cent overload fall 
between the 0 and 100 per cent overload steady creep 
curves. It was found that reasonable calculations of 
cyclic creep characteristics from the steady load creep 
data could be made by assuming that the average creep 
rate is equal to the creep rate due to the overstress plus 
the creep rate due to the normal stress weighted accord- 
ing to time spent at each stress level. 

In applying this simple type of analysis to their test 

Yerkovich indicated that dis- 
calculated specified 


Guarnieri and 
crepancies between 

amounts of creep or times for fracture and the experi- 
mental results were of the order of 25 to 50 per cent. On 


a stress basis, however, this amounted to a difference on 


data, 
times for 


the order of 5 to 10 per cent. 

Two cyclic stress test programs are now in progress 
under Air Force sponsorship at the University of Cali- 
fornia (Berkeley)® and Cornell Aeronautical Labora- 
tories.’ Materials under investigation include 24S-T3, 
75S-T6 aluminum alloys, and FS-1H magnesium alloy 
at temperatures of 300°, 450°, and 600°F. The loading 
patterns consist of one hour loaded followed by one hour 
unloaded cycles, and 1.5 hours loaded and 0.5 hour un- 
loaded cycles such as shown in Fig. 2. The stress levels 
were chosen generally so as to produce fracture in 
several hundred hours. 

A survey of test data contained in progress reports 
on these contracts reveals the following: 

(a) FS-1H: tests at 300° and 450°F. indicate that 
the creep strain is roughly in proportion to the net time 
spent under load. 

(b) 75S-T6: tests at 300°, 450°, and 600°F. indicate 
that the net time to strain a fixed amount for cyclic and 
steady loading is approximately the same. 

(c) 24S-T3: tests at 300° and 450°F. indicate that 
roughly equal times under load produce equal creep 
strains. At 600°F., tests conducted under both con- 
tracts indicate that cyclic loading accelerates creep as 
compared to the effect at lower temperatures. In some 
cases, total test time is about equal for equal strains for 
both cyclic and constant loads. 

Except for 24S-T3 at 600°F. for which there ap- 
pear to be some anomalies in the test data, the test 
results for the materials investigated indicate that to a 
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first approximation, the net time to strain a fixed 
amount for cyclic and constant loading is approxi- 
mately the same. 


LIFE EXPECTANCY ANALYSIS 


The available data on alternating stresses under ele- 
vated temperature conditions indicate that when the 
alternating stress level is a fraction of the mean stress, 
the creep curves remain unchanged to a first approxi- 
mation. Since the alternating stresses due to gust loads 
appear to be less than 50 per cent of the mean stress, it 
appears that gust loading should not decrease to any 
appreciable extent the life expectancy based on the 
mean stress levels. 

Under cyclic load conditions in which the loading is 
applied in a repetitive pattern, the available test data 
for both overloading and underloading cycles indicate 
that to a first approximation the total creep can be con- 
sidered to be composed of the creep resulting from the 
net time under each load. Furthermore, the results 
were not critically affected by the different timing pat- 
terns used. 

These data suggest that for a constant temperature, 
the total strain under alternating and cyclic stress con- 
ditions can be obtained from an equation of the form 


«= Velo, kit) (1) 


where the parentheses (o;, k;r) are used to indicate a 
function of stress o;, relative time #; and life expectancy 
T. 

For the purposes of this analysis, which is concerned 
with relatively long-time creep data, the total strain « 
for a given stress can be assumed to be composed of two 
components as shown in Fig. 1. 


ee = Coy + 6; (2) 


The strain, €0;, is the value at the ordinate obtained by 

extrapolation of the steady state creep portion of the 

curve. It represents the slip resulting from instan- 

taneous elastic and plastic deformation under load and 

the transient creep phase. The other strain component, 

5,, represents the creep associated with viscous flow. 
By substitution of Eq. (2) into Eq. (1) 


«= VY e:(o:, Rit) + Y5,(o;, kit) (3) 


SCIENCES—OCTOBER, 1954 

The first term in Eq. (3) representing slip, requires some 
interpretation, since €, is obtained at the ordinate 
where ¢ = 0. The slip term can be appropriately con- 
sidered as a time average value in the following man- 


ner: 
& = > ev:k; (4) 

By use of Eq. (4), Eq. (3) can be written as 
€= Yd eik:(o;) + >5,(0;, k;7) 5 


To demonstrate the application of Eq. (5) to some 
typical problems, the loadings discussed previously will 
be considered. 


a) Alternating Stress 
The alternating loading can be represented by 
o = o,(1 + a@sin ft) (6) 
where a = o,/0m, Om = Mean stress, and o, = alternat- 
ing stress increment. 
For small values of a and since equal times are spent 
above and below o,,, then to a good approximation 
& = Em 
>5,(0;, kit) = bn 
Therefore, from Eq. (5) 
o = Em + bm (7) 


Eq. (7) indicates that for small values of a, the basic 
creep curve for the mean stress remains unaltered. The 
test data of references 2 to 4 indicate this to be a satis- 
factory approximation for values of @ less than ap- 
proximately one-half. 
b) Cyclic Understressing 

As shown in Fig. 2, the cycle consists of o; for relative 
time k, and zero for ko. Under these conditions, Eq. (5) 
reduces to 

= Gy + 6,ki7 (8) 

The creep curve corresponding to Eq. (8) is shown in 
Fig. 2 relative to the basic creep curve. 
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Fic, 2. Creep under cyclic understressing conditions. 





in 


pe 


fre 


; some 
dinate 
y con- 

man- 


(4) 


some 
ly will 


(6) 
-rnat- 


spent 
ion 


(7) 
basic 
The 
satis- 
1 ap- 


vn in 





c Cyclic Overstressing 
For this case, the cycle consists of o; for relative time 
k, and oz for k, as shown in Fig. 3. Here, Eq. (5) be- 


comes 
e = eR} + €yoko + 5,ky7( 04) + bok2t (a2) (9 


The cyclic creep curve corresponding to Eq. (9) is shown 


in Fig. 3. 
HYPOTHESIS 


In the previous section, attention was directed toward 
determining the total strain after a given time under 
certain simple loading cycles. The method of analysis 
used and the cases considered were based on limited ex- 
perimental evidence under constant temperature con- 
ditions. By contrast, to predict life expectancy of air- 
craft under elevated temperature conditions, it is neces- 
sary to determine the time required to satisfy a given 
failure criterion under complex loading conditions which 
may involve both varying stresses and temperatures. 

The initial information for the analysis of life ex- 
pectancy under elevated temperature conditions con- 
sists of basic creep data for the material and the relative 
frequency distribution of applied stress which the air- 
plane is to experience. Such data are shown in a con- 
venient form in Fig. 4 in which the creep spectrum for a 
given temperature and failure criterion is plotted with 
¢ o, aS a function of time, ft. The term, ¢,, is the short 
time stress at ultimate load factor. Similarly, the load- 
ing spectrum is plotted using the same ordinate versus 
the time kr, where & 1s the ratio of the time at a given 
stress ratio level to the life expectancy 7. To determine 
T, it is necessary to introduce an hypothesis which com- 
bines both the loading and physical property data. 

The hypothesis to be considered proceeds from the 
following assumptions which have as their basis the 
limited experimental! evidence on alternating and cyclic 


loading considered previously : 
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(a) Under constant elevated temperature and varying 
stress conditions, the total creep is due to the net time 
spent at each stress level. 

(b) The varying loads are randomly applied. 

As indicated in Fig. 4, the ratio of time of applied 
stress to time for failure at a given stress level, o;, is 
k, r/t;. o; and all other 


stress levels are zero, then k; = 


In the simple case where o 
1.0 and the life expect- 
ancy 7 is simply equal to ¢;. In more complex cases, 
such as for the loading spectrum depicted in Fig. 4, the 
assumption that the total creep results from the net 
time spent at each stress level leads directly to the 


following cumulative creep hypothesis: 


kit 
> (o,;) = 1 (10) 
t 
By rearranging Eq. (10), the life expectancy 
~a 
dt oF (o,) (11) 


To illustrate the application of Eqs. (10) and (11) to 
creep and loading spectra such as shown in Fig. 4, the 
ratio k, 7/t; for each stress level can be presented in the 
form shown in Fig. 5. The distribution of creep damage 
increments for the hypothetical case illustrated in Fig. 5 
indicates that the curve peaks sharply and thus a major 
proportion of the total creep results from a rather 
narrow stress increment. There is some reason to be- 
lieve that such a distribution of creep damage incre- 
ments may be typical of those to be expected in aircraft 
life expectancy analyses. 

In the most general case, aircraft operation under ele- 
vated temperature conditions will involve both variable 
stresses and temperatures. Although no data are avail 
able under cyclic temperature conditions, it appears 
plausible that the cumulative creep hypothesis may be 
extended to this most general case. 

The initial data for this problem are similar to those 
given in Fig. 4, although the stress-relative time curves 
and the basic creep data for a selected failure criterion 
are given for the series of temperatures, 7);, pertinent to 
The creep damage distributions for each 
The life expect- 


the problem. 
temperature will be similar to Fig. 5. 
ancy hypothesis by extension of Eq. (10) becomes 
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The life expectancy 
k; ia 
r=1/> > — («, T)) (133) 
} 1 L; 


By dividing the creep damage distributions into a series 
of stress increments as shown in Fig. 5, computation of 
the life expectancy can be facilitated by use of a tabular 
summation scheme. 


CONCLUSIONS 


Since creep is a governing factor in determining the 
life expectancy of the aircraft under elevated tempera- 
ture conditions, a review was made of test data under 
alternating and cyclic stress conditions. To a first ap- 
proximation, these data indicate the following: 

(a) For an alternating to mean stress ratio of roughly 
less than one-half, the alternating stresses do not sub- 
stantially affect the creep curve obtained under the 
mean stress conditions. 


SCIENCES—OCTOBER, 1954 

(b) Under cyclic stress conditions, it appears that the 
total creep results from the net time spent at each stress 
level. 

(c) Under cyclic loading conditions, the timing pat- 
tern is apparently unimportant as long as the loading is 
applied in a repetitive fashion. 

Calculations indicate that the alternating to mean 
stress ratio resulting from gust loadings will be less than 
one-half. This, in conjunction with the information 
given in (a), indicates that gust loadings may be rela 
tively unimportant as far as their effects upon life ex 
pectancy are concerned. 

A method of analysis for life expectancy under vari- 
able stress and variable temperature conditions was 
given based upon a cumulative creep hypothesis. The 
analysis is rather simple and direct and is based upon 
the rather limited experimental evidence given in (b 
and (¢c). The proposed method of analysis requires ex- 
tensive testing and verification before it can be applied 
with confidence to the prediction of aircraft life ex- 


pectancy. 
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- . ° ° . . == - yas concentration 
lhe tion, temperature, and velocity distributions throughout the 
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upon entire field of flow. The flame surface is considered to be a source Vv sntahies « F 
. " : e : : : = total jet mass flow 
in (t of heat and a source (or sink) of chemical species dividing the . 
) a ? es . . O = oxygen concentration 
‘ flow into two regions containing either fuel or oxidant . : 
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: : : R = transformed independent variable 
so that only the equations of motion need be solved. It is also R grodact concentration 
i - m at = Tor ct conce allo 
found that such solutions can be made to satisfy the usual bound 1 semmperatue 
sq . rT. = e CTU c 
ary conditions at the flame front and the conservation laws. The 
solution of the equation of motion for the case of variable density Subscripts 
Third has been given by Howarth in the two-dimensional case, and an , mn ee 
951 approximation is developed for the axially symmetric case s = reference point 
f Tem- Under these conditions it is found that the diffusion flame tem = atmosphere far from jet 
ASTM perature is the same as the adiabatic flame temperature of a 
stoichiometric fuel-atmosphere mixture. 
ic and rhe length of the diffusion flame is calculated and compared (1) INTRODUCTION 
npera- with some experimental values for the axially symmetric case 
rhe agreement is not univ ersally good, and possible causes of the HE SIZE OF LAMINAR DIFFUSION FLAMES burning 
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is determined once the velocity field is known.* The 
known solutions of the velocity field of a jet are dis- 
cussed, and an approximate solution for the axially 
symmetric heated jet is given. Finally, the flame 
lengths calculated by this method are compared with 
some previously published experimental values. 


(II) THe RELATION BETWEEN VELOCITY 
AND TEMPERATURE OR CONCENTRATION IN A LAMINAR 
JET OF VARIABLE DENSITY 


The equation of motion of a laminar jet of variable 
density assuming no pressure gradients or gravity 
forces is, for the two cases of two-dimensional and 
axially symmetric motion, respectively, 


Ou ‘ Ou 1 Oo ( ~ 
iu v => a) 
Ox Oy p Oy “j oy = 


Ou és Ou 10 ( = . 
1 v = r ) 
Ox or pr Or . or — 


and the corresponding equation of continuity, 
(O/Ox) (pu) + (O/Oy) (pv) = O (2a) 
(O/Ox) (pur) + (0/O0r) (per) = 0 (2b) 


If the jet gas composition is not identical with the 
surrounding atmosphere, there will be a mixing of both 
caused by diffusion. The diffusional flow of any spe- 
cies 7 in a velocity field of the type considered is deter- 
mined by the diffusion equation, 


oc; rm Ww; 410 (v =) 

l v = im 3¢ 
Ox oy p Oy Oy ata 
OC; 4 Oc; i @ (v 4) ' 

t v - im pr a 
Ox Oy pr Or M Or _ 


where c; is the mass fraction and thermal diffusion is 


1 


neglected. This equation holds in the region where no 
sources exist. 

The corresponding energy equation becomes (again 
for the region where no heat sources exist) 


ra) re) 
pu (cpl) + pv (c¢,7) = 
x Oy 


0: 
a) (. or) + (*)’ ie) 
oy Ov “ oy 1“ 


_—__ a 
pu “ae (Cp) + pv oy (Cpl) = 


Ld (r *) in (“) (4t 
) 
r Or f Or ‘ or ; 


* A similar conclusion is drawn by Schwab’ for the case of 
turbulent diffusion flames, under the assumption that the mean 
velocity is determined by the momentum transport theory. The 
author is indebted to G. Markstein and D. Olshevsky, of the 
Cornell Aeronautical Laboratory, Inc., for calling his attention 
to this paper and providing a translation 
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Under the condition that D;,,o = yu it can readily be 
seen that a solution of the diffusion equation, Eq. (3), 
is 
c= B+ Cu 5 


where B and C are arbitrary constants determined by 
the boundary conditions. For common gases, D;;p/u 
varies between 1.3 and 1.6, where D;; is the coefficient 
of self diffusion; hence, this condition is not satisfied 
exactly in any real case. 

A corresponding simple solution of Eq. (4) may be 
found if the Prandtl Number c,u/\ is unity and the 


specific heat c, is constant, 
Cpl + (u?/2) = E+ Fu (6 


[This was first shown by Crocco® for the two dimen- 
sional case, Eq. (4a), and by Pai’ for the axially sym- 
metric case, Eq. (4b).] For common gases the Prandtl 
Number is somewhat less than one, and the specific 
heat is reasonably constant for relatively small changes 
in temperature. 

While it is possible that solutions of Eqs. (3) and (4 
could be obtained in cases where the above restrictive 
conditions are not exactly met, the simplicity of the 
above solutions, which reduces the problem to solving 
Eq. (1) and evaluating the constants of integration in 
Eqs. (5) and (6) from known boundary conditions, 
makes them extremely attractive. As in the case of 
heat transfer in a laminar boundary layer, the devi- 
ation of the Prandtl Number (or the ratio Dp/y) from 
one may not greatly effect the solution. 


(III) THE APPLICATION OF JET SOLUTIONS TO THE 
DIFFUSION FLAME 


Under the simplifications outlined above, once the 
solution of the equation of motion, Eq. (1), is deter- 
mined, the temperature and concentration fields may 
be of the forms of Eqs. (5) and (6), respectively. It is 
not immediately apparent, however, that the particular 
boundary conditions for the diffusion flame may be 
satisfied by solutions of this form, even though, for 
example, such solutions are suitable for the heated jet. 

Burke and Schumann? first proposed that the flame 
surface, which is extremely thin, occurs in that portion 
of the flow where the diffusional flow of the fuel from 
one side and the oxidant from the other are in the 
stoichiometric proportion. Consequently, in the simple 
jet diffusion flame, for example, the gases in the region 
“inside”’ the flame surface are composed of fuel, prod- 
“outside’’ contain 
oxidant, products, and inert gases. The fuel diffuses 
toward the flame front, while the oxidant 


ucts, and inert gases, while those 


outward 
diffuses inward toward it. 
flame front diffuse in both directions away from the 
This picture of the distribution of fuel, 


The products formed at the 


flame front. 
products, and oxidants was confirmed by the measure- 
ments of Hottel and Hawthorne.*® 
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LAMINAR JET 


Since the chemical reaction rate in a stoichiometric 
mixture is much greater than the rate at which the fuel 
and oxidant mix by diffusion, it is the latter process 
that governs the position of the flame front, and the 
thickness of the front (which is the region where both 
fuel and oxidant are present simultaneously) is thus 
much smaller than the other dimensions of the flow.* 
It is the argument of Burke and Schumann? that the 
flame front may be idealized as an infinitesimally thin 
surface in the flow which is a sink for both fuel and 
oxidant and a source of products. 

Assuming that this is the correct picture of a simple 
jet diffusion flame, then the boundary conditions on the 
concentration and temperature distributions may be 
summarized as follows: 

(a) At the flame surface, the concentration of fuel 
and oxidant falls to zero and the gradients are in the 
stoichiometric proportion (since the diffusivities are 
assumed equal for all components). 

(b) At infinity the concentrations of the products 
and any inert diluents added to the fuel is zero, while 
the concentration of the oxidant and inert gases reach 
their atmospheric values. Furthermore, the temper- 
ature also equals the atmospheric value. 

(c) The first space derivatives of the concentrations 
and temperature in a direction transverse to the jet 
axis, evaluated on the axis of symmetry and at infinity, 
are zero like those of the velocity. This follows from 
symmetry and from the fact that there are no sources 
of particles or heat on those boundaries. 

(d) There must be conservation of mass—1.e., the 
mass flow of each atomic species across any plane nor- 
mal to the jet axis must be constant. 

(e) There must be conservation of energy in that the 
enthalpy flux must be similarly conserved. 

It is evident that condition (c) is automatically 
satisfied by solutions of the type of Eqs. (5) and (6) for 
any value of the constants, while (a) and (b) may be 
satisfied by appropriate choice of these constants. 
However, it is not a priori evident that (d) and (e) 
can be satisfied by such a choice of constants. 


(1V) THE DETERMINATION OF CONCENTRATION 


Consider a simple jet diffusion flame of a mixture of 
fuel (F) and diluent (D) exhausting into an atmosphere 
of oxygen (QO) and an inert gas (J). The concentra- 
tions of these four components are chosen to be as fol- 
lows: 


Inside the flame surface: 
[D] = ayu l 


[F = aldo — e€) 


il] = 1 eam’ 





* Measurements of the thickness and the spectroscopic char- 
acteristics of some diffusion flames have been reported by Wolf- 


hard and Parker."! 
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Outside the flame surface: 


[O] = € — dol 
[7] = 1 — € — agu 8) 
[D] = ayu 


where a is the mass ratio of fuel to oxygen in a stoichio 
metric mixture and ¢ is the mass fraction of oxygen 
in the atmosphere. 

This choice satisfies conditions (a), (b), and (c) above 
with the constants a, still undetermined 

Since the sum of all the concentrations must every 
where be unity, the concentrations of products (R) 
is given as follows 

Inside the flame surface: 

[R] = e(1 + a) — (a; + aa (3)U Q) 

Outside the flame surface: 

Q,)u 10) 


[R] = (d2 + a 


It is now possible to write the integral relations for 
the conservation of the various species. Considering 
a plane normal to the jet axis, the flow of fuel particles 


is given by 


Ay j a 
p jaladou — €) + fe(l + a 
Jo | l+a 


| . 
(a; + ad, —a ul udA +4 pX 
. 1 
a 1) an 
(a2 + a3 — a,)u-dA 
l+a 1+, 
where 
MJ = mass flow in the jet nozzle 
8 = massratio of diluent to fuel 
dA = 2nr dr for the axially symmetric flow 


= 2 dy for the two-dimensional flow 


The first integral is taken over the region inside the 
flame surface and the second over the exterior region. 


Upon simplification this equation becomes 


a P VW 
(do + Az; — Q}) } pit dA 
1+ a/ Jo 1+, 


But the integral on the left is just the momentum of 
the jet and, hence, equals the momentum at the mouth 
of the nozzle (P) since there are no axial forces or pres- 


sure gradients. Thus, 
ado +a —-a =(l+a)M all +3)P 11) 
where 
ZW 


+ The diffusional flow is neglected in comparison with the con- 
vective flow, since it is easily shown to be of the order of (Rey- 
nold’s Number with respect to the latter The same argu- 


ment applies to the enthalpy flux 
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Similarly, for the conservation of the diluent, 


1A BM 
P\Qu)u dA = 
J0 1+ 2, 
or 
oM 
aq = (]? 
(1+ B)P 


The flows of oxygen and inert constituents are, re- 


spectively, given by 


7A 
Oo = | E 
J 0 l 7 a 


“ l 
| pl é€ — Qu + (dz + Q3 — an)e | udA 
ow As l + a an 


| p[(1 — €) — a3uju dA 
0 


fe(1 + a) — (a; + ade — a3)ujudA + 


¢; = 
which, upon simplification, becomes 


— Q@Qd2x — aj) 


l 
Qo = 26) + rer (Qs 


(1 — a1) = a3P 


[ pu dA 
/7 0 


The integral Q, in the case of the incompressible 


2; = 


where 


Q= 


jet, increases with distance along the axis* and, hence, 
is not constant--1.e., the jet entrains material at its 
outer edges. However, since this material is drawn 
from the atmosphere, the ratio of the flow of oxygen to 
inert gas must be the same as that in the atmosphere, 
or 

«OQ + (1 (1 + a@)] (as — aa, — a,)P € 


(1 —_ 6)O = a3P l —.-§ 
which simplifies to 


(1 + ea)az — (1 — €) (ade + a;) = O (13) 


Eqs. (11), (12), and (13) may now be solved simul- 
taneously for the constants a» and a; to give 


dy = {{1 + ea(1 + B)]/a(1 + B)} (M/P) (14 


az = (1 — e) (M/P) (15 
which, together with Eq. (12), may now be substituted 
in the expressions (7) through (10) for the concentra- 
tions of all the species. 

The flame surface is located by noting that the fuel 


(and oxygen) concentration on it is zero-—-that is, 


€ ea(l + £) P 


: 16 
as 1+ ea(1 + B) M 


uy = 


Hence, the flame surface is also a surface of constant 


velocity. 
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(V) THE DETERMINATION OF TEMPERATURE 


Since the flame surface is the source of heat, as well 
is products, the temperature will be greatest at this 
surface and will decrease on either side of it. Hence, 
the following is chosen as the temperature distribution 
inside and outside (respectively) the flame surface: 


T = ds — d5u (17) 


T= T. + deu (18) 
The term in Eq. (6) involving u* has been neglected, 
since the flow velocity is extremely low. 

Consider next the enthalpy flux through a plane nor- 
mal to the jet axis at the tip of the flame. For this 
special case Eq. (18) holds throughout the entire plane 
and on the jet axis w has the value of uw, given by Eq. 
(16). The enthalpy flux // is given by 


H = { c,T pu dA 
0 


Cyl oO + AelyP (19) 


where Q and P are the flow and momentum integrals 
previously considered. Since the first term in Eq. (19) 
equals the enthalpy of the atmospheric gases entrained 
in the jet, the second term must equal the enthalpy flux 
Hence, 


ae = (h, Cp) *(M P) 


of the fuel. 


(20)* 


where /y is the heating value per mass of the fuel- 
diluent mixture. 
The flame temperature is now determined from Eq. 
IS) with w = u,, whence 
ms a h, 1 ea(l + 2B) 
T,=T. + : (21) 
& LI + ea{l + 6) 
* The specific enthalpy of the fuel, cp7, is neglected in com- 
parison with the enthalpy content due to chemical reaction 
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which is the same flame temperature as that for a 
stoichiometric premixed flame using as components the 
atmosphere and fuel-diluent systems of the diffusion 
flame. 

To determine da, and ds, it is necessary to write the 
equation of enthalpy flux across a transverse plane 
upstream of the flame tip, 


1; 
[ 1C,(ds — ast) + ahs(1 + B) (aou e)} pudA 4 
7 


Q+ascpP (22) 


+ agu)pudA = c,T 


This equation may be satisfied by letting 
C,fl4 — GE hh, + 8) =o, 


i 


ah + B)de ees Cpls = Cpe 
which simplifies to 


a,j = T.. + ae(1 + 8B) (hy Cy) (23) 


« 


ds = ae(l + B) (hy/c,)-(M/P) (24) 


Thus the temperature distribution inside and outside 
the flame surface becomes, respectively, 


; P h, (: M ) 
- ael + | — 95) 
tT ae P u“ (Zo 


I 


+ (hy/cy)-(M/P)-u (26) 


which expressions are equal at the flame front (u = uy). 


(VI) SOLUTIONS OF THE EQUATION OF MOTION 


Having shown that the concentration and temper- 
ature fields in a diffusion flame may, under certain 
simplifying assumptions, be expressed as linear func- 
tions of the velocity, it only remains to solve the equa- 


tion of motion. 
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lengths with measurements by Wohl, Gazley, and Kapp‘ for 
some city gas- and butane-air mixtures. 


The solution for the laminar incompressible jet (see 
reference 8) is not very useful since, because of the high 
temperatures generated at the flame surface, the den- 
sity may vary by a factor of ten or more throughout the 
field. A solution to the two-dimensional laminar jet 
for variable density was first given by Howarth’ (under 
certain restrictions) but no similar solution for the axi- 
ally symmetric case has been published. 

Howarth’s solution is valid only if the viscosity 
varies with the first power of the temperature—i.e., 


uw = us (T/T;) 


where yp; is the viscosity at a reference temperature 7’. 
Although the elementary kinetic theory gives this 
dependence of viscosity on temperature, experiments 
show that » varies as T to a power less than unity, 


usually 0.7 or 0.8. By choosing the reference tem- 


SCIENCES—OCTOBER, 1954 

perature to be that at the point of maximum transverse 
velocity gradient, the error in the viscous force will 
generally be small compared to its maximum value (at 
the reference point) for those regions where the actual 
viscosity departs appreciably from the assumed value 
Of course, since the temperature at the point of maxi 
mum transverse velocity gradient changes with the 
x coordinate, such a solution can only be expected to 
be a good approximation over a limited region of the 


whole flow field. 


For the case where un = u,(7'/7;), Howarth’s solu- 
tion is* 
P? ‘/s 
“= 0.4543 ( ) sech? £ 
MsPsX 
where & = 0.2752 (Pp,/p,?x?)'*? Y and 
l Y 
y= 4 | Tr dY (27) 
foe SF 
which is not valid near x = 0. 


Pai’ gives a method of successive approximations 
for the two-dimensional case when np = u,(7/7,)”. 
It is possible that a numerical solution of this type could 
be obtained for the particular boundary conditions of 
a diffusion flame. 

There is as yet no exact solution of the Howarth 
type for the axially symmetric compressible laminar jet, 
although a numerical method has been suggested by 
Pai.’ Following the Howarth method, an approximate 
solution for this case is given in the Appendix. 


(VII) FLAME LENGTH 


Provided the flame length is many times the width 
(or diameter) of the mouth of the jet, so that the simple 
jet solutions are valid near the flame tip, the length L 
is easily obtained by noting that, at x = ZL and y (or 
r) = 0, the velocity « = u,. For the two-dimensional 
Howarth solution, 


L = 0.0938 (P?/ psusit;*) 


and for the axially symmetric approximate solution 
(Appendix) : 


3 P [p, p. ' 
L = I 
Sm pwsty \ps py 
Since uy is given 


where J(p,,/p,;) is shown in Fig. 3. 
above in Eq. (16), these become, respectively, 


L 0.0938 pi + ca (1 + 6)” 28) 
,= UdI0d 47 P (Z 
; ex(l +) §$ PsbsP 


3 jl + ea(1 + B)) M p, | (*) , 
Lz < > . I (29) 
Srl ea(1 +B) Sus p, pr 


* In the Howarth solution, all the derivatives of u with respect 
to Y are continuous due to the transformation on y, although 
021 /Oy? is discontinuous at the flame front due to the discontinuity 
in O7/Oy 
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LAMINAR JET 


Since ./ and P are, respectively, the mass flow and 
the momentum at the mouth of the jet, it can be seen 
that in both cases the length increases linearly with the 
average flow velocity at the mouth and with the square 
of the mouth width (diameter). 

(VIII) FLAME LENGTH: COMPARISON WITH 
EXPERIMENT 


The approximate solution of the equation of motion 
for the axially symmetric case (see Appendix) 1s re- 
quired if comparison with experimental measurements 
of flame length is desired, since no measurements have 
been published for two-dimensional flames. Thus, 
the applicability to diffusion flames of the exact Ho- 
warth solution cannot be tested without further experi- 
mental work. 

In Figs. l(a) through 1(e) are shown the theoretical 
flame lengths for some city gas- and butane-air mixtures 
burning in air as compared with the measurements of 
Wohl, Gazley, and Kapp.’ 
the jet gases (viscosity and specific heat) were assumed 
The reference temperature for de- 


The average properties of 


to be those of air. 
termining the viscosity is chosen to be the temperature 
at the point of maximum drag in a transverse plane at 
the flame tip, as found in the Appendix. 

Some of the experimental results of Rembert and 
Haslam! are given in Fig. 2. The measured values for 
large air-gas ratios have not been used, since in such 
cases the inner cone length is a large fraction of the total 
flame length, a type of flow that is not reasonably ap- 
proximated by the solutions previously considered. 

In considering the discrepancies between theory and 
experiment, it is best to recall the restrictions placed 
upon the theoretical results by the following assump- 
tions that were made: 

(1) There are no gravity forces—1.e., natural con- 
vection of the hot gases is neglected. 

(2) The diffusion coefficients of all the species are 
the same and are independent of the mass fraction. 

(3) The diffusivity is equal to the kinematic vis- 
cosity at each point of the jet. 

(4) The diffusivity is proportional to the first power 
of the temperature. 

The first three assumptions are necessary for solu- 
tions of the type where the concentrations and tem- 
perature are linearly related to the velocity. The last 
is necessary for the solution of the equation of motion, 
which in the axially symmetric case is only an approxt- 
mate solution. 

The effect upon the theoretical flame length of the 
failure to meet each of these restricted conditions in 
any actual case can best be determined by more elabo- 
rate analyses. What should also be considered is the 
degree to which the actual jet flames in the experi- 
ments considered approach the flow conditions pre- 
sumed for the laminar jet problem. 
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In two important respects there is a difference be 
tween the idealized and the actual flow. In the first 
case, the theoretical solution assumes no axial velocity 
along the plane x = 0, except at the mouth of the tube. 
In actual experiments the air is permitted to move 
freely in the axial direction past the end of the tube 
under the influence of convective forces. This effect 
would tend to increase the length of the flame because 
of the greater flow velocities. 

An even more important effect is the failure of the 
flow to remain 
Shadowgraph pictures by Wohl, Gazley, and Kapp‘ 
show that in some instances the flow in the outer re- 


laminar throughout the entire jet. 


gions of the jet becomes unstable and eventually be- 
comes turbulent, even though the core flow (in which 
the flame front is located) may remain laminar. The 
instability may be due to convective forces or to condi- 
tions at the lip of the burner tube. In any case, the 
failure of the flow to remain laminar throughout will 
materially affect the velocity and, hence, the flame 


shape. 


APPENDIX 


AN APPROXIMATE SOLUTION OF THE AXIALLY 


SYMMETRIC LAMINAR JET OF NONCONSTANT DENSITY 


In this section an attempt is made to find a solution 
of the equation of motion for the case of uw = yu,(7/7,). 
It is found that the usual type of similarity solution does 
not satisfy the equation of motion exactly, but only 
approximately. 

The equation of continuity (26) is satisfied by de- 
fining the function y¥ such that 
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pur = p(Oyv/or) | (30) 
pyr = —p,(Ow/dx) I 


Before writing the equation of motion in terms of y, 
define a new independent variable R(x, 7), where 


pR=2 [ pr dr (31) 
0 


the integration being performed for x held constant. 
In terms of the independent variables R and x, the 
equation of motion, Eq. (16), becomes 


] 
Vrixr — Vx (Yu oo vn) = 


Ms p ) (£)'{oe - del (39) 
. naaec”™hUCUCrr?hCUT 


A similarity type of solution may be tried. Assume 
WY = (us/p.,)x 2(n) (33) 
where 
a = V p, bs (R/x) (34) 
This makes 
u = (p,/p.)* (2'/nx) (35) 


For the case of a heated jet, the temperature is re- 
lated to the velocity by 


T/T, =1+ Du 


within the assumptions previously made that c,u/A is 
unity, c, is constant, and u* may be neglected. For a 
perfect gas, the density is inversely proportional to the 
temperature, so that the equation relating r and R be- 
comes 
d(r*) = (p,/p)d(R?) 
(1 + Du)d(R?) 


For the similarity type of solution this integrates to 


r\2 2D p, \? g 
= ] : 2 > 
(“.) “+ = (’ ) 7° (36) 


Hence, the equation of motion becomes 


become 


e-“e +7) =~ x 


e-MC'DH o 
g — — ss (94) 
' ” x \p./ n/$ 


Since this equation contains x explicitly, g cannot 
be a function of 7 alone, as was postulated for this simi- 
larity solution, but must be a function of both x and 
n. If it is assumed that the partial derivatives of g 
with respect to x are small compared with those with 
respect to 7, then the former may be neglected in an 
approximate solution and Eq. (37) thus may be inte- 
grated holding x fixed. This approximate solution has 
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at least the virtue that it reduces to the exact solution 
for the constant density case. 

In the two-dimensional solution of Howarth, the 
only effect of a nonuniform density is to stretch the y 
coordinate,* the dependence of the solution on x (when 
y = O) being unchanged from the constant density 
case. In the axially symmetric case this is no longer 
true (since the similarity type of solution fails), yet the 
solution should not radically change its nature, so that 
the principal effect should be a stretching of the ¢ co 
ordinate, with some modification of the dependence 
upon x. 

The approximate equation of motion, Eq. (37), be- 
comes a better approximation near the edge of the jet 
(n very large) or far downstream (x very large). 


Eq. (37) may be rearranged to give 


1 {gg’ d g’ J 
= ( jae ee \(e" ee. \Q if. g a} (28 
dn ' n f dn { n n* f 


where g is a constant to be determined. As in the in- 
compressible jet, this may be integrated once, the con- 
stant of integration disappearing due to the boundary 


conditions at 7» = O (see reference 8) to give 


ge’ tvs £ 
-(:+98)(e-9) 
U} r n 


For the purpose of numerical integration, this may be 


39 


simplified by introducing the variable 

z= pn’ (40 
where / is an undetermined constant, with the result 
that 


—gg’ = 2[1 + gp*(g/z) Jzg” (41 


where the primes refer to differentiation with respect 
to 2. 

From Eq. (31) it can be seen that, as 7 approaches 
zero, the ratio r?/R? becomes 


r?/R? ~ p./ po 


where pp is the density at r = 0. Since the quantity 
in brackets in Eq. (41) is r*/R’, then, by choosing g 
such that g/z is unity when z = 0, g may be found, 
giving 


, 


—gg = 232 + [(p na = l]g}g” (42 


This equation must be integrated numerically except 


when pp = p,, (constant density case), for which 


g = 2/[1 + (2/4)] (43 
The constant » may be determined from the equation 
for the total jet momentum, 


* At a given value of x, the value of y for which the velocity 
is a given fraction of the centerline velocity (y = 0) is larger 
in the heated than in the unheated jet 
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= 4ru, (p;/p,)* p” (g’)* dz (44) 
I (°°) =~ [ (g’)? dz (45) 
Po 1 J0 
3P + ; 
ZOO 
lOmu, \ps pr 


Finally, the velocity « may be determined as 


ar (’ ) , [ p , 
“= . g'(z |! (47) 
SWuUsX \ ps po 


where g’(0) = 1 gives the velocity on the axis. 


Letting 


then 


The reference point s may be chosen to be that at 
which the drag is a maximum-—that is, where ur(Ou 
Or) is greatest. It can readily be shown that 

Ou ; 
ur — « gg 
Or 


so that the maximum value of gg’ defines the point of 


maximum drag and so that the density p, and also yu 
may be determined. 

The approximate solution thus outlined gives the 
velocity profile in a plane normal to the axis when the 
density po on the axis is known. The reference point 
s is somewhat arbitrary and may be chosen to suit the 
particular region of the flow which is of interest. (In 
the case of the diffusion flame, s was chosen to be the 
point of maximum drag in the normal plane at the 
flame tip.) 

In Fig. 3 are shown the values of /(p,/p 9) and 
p.,/ps (where s is the point of maximum drag), calcu- 
lated from the numerical integration of Eq. (42). It 
can be seen that the velocity along the axis of a heated 
jet is less than that for the unheated jet of equal mo- 
mentum at identical axial positions. 


REFERENCES 


!Rembert, E. W., and Haslam, R. T., Factors Influencing 
Length of a Gas Flame Burning in Secondary Air, Industrial and 
Engineering Chemistry, Vol. 17, No. 12, pp. 1286-1238, Decem- 
ber, 1925. 

2 Burke, S. P., and Schumann, T. E. W., Diffusion Flames, 
Industrial and Engineering Chemistry, Vol. 20, No. 10, pp. 998 
1004, October, 1928 

’ Hottel, H. C., and Hawthorne, W. R., Diffusion in Laminar 
Flame Jets, Third Symposium on Combustion, Flame and Ex- 
plosion Phenomena, pp. 254-266, Williams & Wilkins, Baltimore, 


1949 


DIFFUSION 


FLAME 689 


#) 


1 - 








Momentum integral and density ratio at point of maxi- 
mum drag in heated axially symmetric jet 


Fic. 3 


* Wohl, K., Gazley, C., and Kapp, N. M., Diffusion Flames, 
Bulletin No. 1, University of Delaware Engineering Experiment 
Station, Newark, Del., 1948 

5 Schwab, V. A., Relations Between Velocity and Temperature 
Fields of a Gaseous Torch, Journal of Technical Physics, U.S.S.R., 
Vol. II, No 5, pp 431-443, 1941 

6 Crocco, L., Su di un valore massimo del coefficiente di tr 
misstone del calore da una lamina piana a un fluidu scorrente, 
Rendiconti R. Accademia dei Lindei, Vol. 14, Fase. 490-496, 
1931. Also sulla transmissione del calore da. una lamina piana un 


ins- 


fluido scorrente ad alta velocita, L’Aerotecnica, Vol. 12, Fasc. 181 


197, 1932 

? Pai, S. L, 
Fluid, Quarteriy of Applied Mathematics, Vol. 10, No. 2, pp 
141-148, 1952 

8 Goldstein, S., Modern Developments in Fluid Dynamics, pp 
145-148; Oxford, 1938 


Concerning the 


Compressible 


7 , ‘op ref , " 
lxially Symmetric Jet Mixing of a 


* Howarth, L Effect of Compressibility on 
Laminar Boundary Layers and their Separation, Proceedings of 
the Royal Society, A, Vol. 194, pp. 16-42, 1947 
I., Two-Dimensional Jet 
Fluid, Journal of the Aeronautical Sciences, Vol. 16, No. 8, pp 
163-469, August, 1949 

Wolfhard, H. G., and Parker, W. G., “* 
tigation into the Structure of Diffusion Flames,” 
the Physical Society, London, A, Vol. 65, No. 385, pp. 2-19, 1952 


1 Pai, S Vixing of a Compressible 


A Spectroscopu Inves- 


Proceedings of 











Effect of Leading-Edge Separation on the 
Lift of a Delta Wing 


C. E. BROWN* ano W. H. MICHAEL, Jr.t 
Langley Aeronautical Laboratory, NACA 


SUMMARY 


There are a large class of important fluid flows which involve 
the effects of fluid viscosity and yet do not constitute flows having 
restrictive amounts of viscous dissipation. Examples of such 
flows are the flow over slender bodies in which vortex shedding 
occurs, swept wings having separated leading edges and so forth 
The important feature of these flows is the production of free vor- 
tex sheets in which large dissipation ultimately occurs but only at 
great distances downstream of the generating body; hence, these 
flow fields may be approximated by potential flow solutions. 
The present paper presents an application of this technique to 
slender triangular wings having leading edge separation. Curves 
are presented showing the lift versus angle of attack, pressure 
distributions, and integrated span loadings. It is shown that 
leading edge separation leads to nonlinear lift curves with lift 
greater than that obtained from potential flow theory. Compari- 
son of lift curve with experiment is made. 


INTRODUCTION 


rT RECENT YEARS small perturbation theory for com- 
pressible flow has been developed and refined to a 
point where estimates can be made of extremely com- 
plex flow fields and quantitative results can be obtained 
which are often in good agreement with reality. How- 
ever, the theory is limited in use by two important re- 
strictions. The first is the restriction to small perturba- 
tions which allows only a first approximation to the 
effects of compressibility. The second is the restriction 
to flows in which the dissipation due to fluid viscosity is 
negligible. Removal of the limitations requires working 
with nonlinear differential equations with their at- 
tendant great mathematical difficulties. There are, 
however, a large class of important fluid flows which 
involve the effects of fluid viscosity and yet do not con- 
stitute flows having restrictive amounts of viscous dis- 
sipation. Examples of such flows are the flow over 
slender bodies in which vortex shedding occurs, swept 
wings on which separation occurs ahead of the trailing 
edge with the formation of free vortex sheets, etc. The 
important feature of these flows is the production of 
free vortex sheets in which large dissipation ultimately 
occurs but only at great distances downstream so that 
in the vicinity of the body or wing in question potential 
flow theory may be applied with suitable alteration of 
the usual boundary conditions. This technique is not 
new to fluid mechanics, having been used by Kirchhoff! 

Presented at the Aerodynamics Session, Twenty-Second An- 
nual Meeting, IAS, New York, January 25-29, 1954. 

* Head, Supersonic Aerodynamics Section. 

t Project Engineer, Supersonic Aerodynamics Section. 


and von Karman’ to predict the drag of two-dimensional 
shapes; the Karman vortex street represents one suc- 
cessful application. Prandtl, Kaden, and Anton*® in 
vestigated the production of vortex sheets at sharp 
edges; however, their interest was largely in the shape 
and position of the vortex sheet and not in the reaction 
on the generating surface. 

In the present paper, the flow over a slender delta 
wing is considered, and the presence of leading edge 
separation is assumed. The physical flow field over the 
wing is expected to look like the schematic drawing 
shown in Fig. l(a). That is, separation occurs on the 
leading edge and produces the two spiral vortex sheets 
across which the pressure is continuous but the tan- 
gential velocity is discontinuous. For the slender delta 
wing the field is conical since there is no characteristic 
dimension on which to base variations in any quantity 
along conical rays through the origin. The only contri- 
bution of viscosity in these flows is to fix the separation 
point at the leading edge for reasons exactly analogous to 
those justifying the use of the Kutta condition at sub- 
sonic trailing edges. Solution of the problem consider 
ing the spiral vortex sheet was found to be too difficult 
and hence a simplified model, more amenable to calcula- 
tion, was adopted. This model shown in Fig. 1(b) re- 
places the spiral sheet with two concentrated line vor- 
tices above the wing and two feeding vortex sheets con- 
necting the source of vorticity (leading edge) and the 
concentrated line vortices. It is expected that the re- 
sults of the simplified model calculations should give a 
fair estimate of the forces acting on the wing and bring 
out the important features of flows involving separated 
leading edges. 

Since the completion of this work, a paper by 
Legendre® dealing with the same problem has been 
brought to American attention by Adams.?7 Adams? 
pointed out in the Readers’ Forum of the Journal of 
Aeronautical Sciences that Legendre’s solution failed 
to properly account for the forces on the feeding vortex 
sheet and that inclusion of the sheet forces produced a 
result which in effect left an uncancelled finite force in 
the flow field over the wing. In the present paper this 
ambiguity is resolved by use of more appropriate bound- 
ary conditions. Clarification of this point will be made 
in the following sections. 


t Acknowledgment is made to Mr. Mac Adams for his work 
with the senior author on earlier formulation of the problem. 
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SYMBOLS 


, = coordinate along wing in direction of free stream 
= coordinate along wing normal to free stream 
= coordinate normal to wing surface 
= disturbance velocit\ pote ntial 

,w = disturbance velocities in y and ¢ directions, respec 

tively 

VU = Mach Number 

= free stream density 

V = free stream velocity 


= vortex strength (circulation ) 


z = (y + 12) 

J = value of o at vortex point 
= half breadth of wing at any chordwise station, x 
=VY-! 


= mean normal flow velocity over vortex 


= half apex angle of wing 


IW (oe) = complex velocity potential 
0 = vector point coordinate (see Fig. 2 
x = angle of attack 
Ap/q) = (Plocat — Pfree stream)/[(1/2)p V? 
Re = real part of 
” = apparent mass 
l = lift 
ANALYSIS 


The problem to be considered is that of the potential 
flow about a slender delta wing on which there exists 
leading-edge separation. That is, the streamlines of 
the flow which wet the wing do not pass from the 
lower to the upper surface but rather come from both 
surfaces and leave at the leading edge. In reality such 
a condition would produce a conical spiral vortex sheet 
above the wing and the boundary conditions of the 
problem would be that no fluid pass through the wing 
surface and that the pressure across the vortex sheet be 
continuous. Clearly this represents a difficult problem 
since the solution must provide both the shape and 
strength of the sheet. Past experience with vortex 
sheets leads one to hope that the main features of the 
flow can be obtained by replacing the spiral sheet by a 
concentrated vortex near the center of the spiral. For 
conical flow, however, the net vorticity in the spiral is 
linearly increasing in the downstream direction, hence 
the concentrated vortex must be of linearly increasing 
strength also. The increase in strength must be accom- 
plished by a feeding sheet of vorticity in order to 
satisfy Kelvin’s theorem, and thus we fix on the model 
shown in Fig. 1(b) as the most appropriate flow field 
amenable to simple calculation. 

The equation of motion to be satisfied is the well 
worn equation which represents a slightly perturbed 
iain stream of velocity |’ and corresponding Mach 
Number, 7: 


(1 — 3M? )err + Sy + G22 = 0 (1) 


zr Vu 


where g is the disturbance velocity potential, and 
v, y, 2 are Cartesian coordinates fixed to the wing; z 
is measured perpendicular to the flat wing surface, y 
tangential to the wing surface but normal to the free 
stream velocity vector, and x is measured in the flow 
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direction along the plate. If we further restrict the 
flow to highly swept wings, the term ¢g,,(1—./*) may 
be neglected and the equation of motion becomes La- 


place’s equation in y and 2z: 
Puy + ¢::=0 


We therefore use what is commonly known as the 
slender body theory. 

Boundary Conditions.—The conditions at the plate 
are first that the plate is solid and hence the normal 
velocities are zero, and second that the flow separates 
at the plate edges. The conditions in the fe/d are first 
that the disturbances vanish at infinity and second that 
the fluid pressure is continuous. This latter condition 
is however, impossible to satisfy with the assumed 
model and hence must be replaced with one which is 
more compatible. The difficulty lies in the presence of 
the feeding vortex sheet across which there must be a 
pressure discontinuity, but since the assumed vortex 
system represents the true spiral only at a distance it is 
to be expected that in the small regions near the system 
violation of natural conditions might occur. It is, 
therefore, necessary to make the last boundary condi- 
tion less detailed and hence we only require that the 
integral of pressure around the assumed vortex system 
vanish. In simpler terms, we require as a final condi- 
tion that the assumed vortex system (feeding sheet and 
concentrated vortex) have zero net force acting since 
only the wing and not the fluid can sustain forces. Ap- 
plication of this idea to the model then requires that the 
forces on the feeding vortex sheet be cancelled by equal 
but opposite forces on the concentrated vortex; thus 
the concentrated vortex is not force-free as was as- 
sumed in the Legendre solution. The mathematical 
formulation of the preceding boundary condition is as 
follows: at a given station x where the half plate width 
is a, we introduce the quantity o = y + 172, the vector 
distance to a point y, z. The feeding sheet vector force 
per unit length may now be written 


. arf a ; 
ip | (oo — a)T (2) 
dx 


where p is the stream density, oo the location of the con- 
centrated vortex and y is the circulation of the vortex. 
The vector force on the concentrated vortex must be 
produced by a fluid flow normal to the vortex of vector 
velocity of v*. The vector force would be —ipv*y. 
Setting the vector sum of the two forces to zero thus 


. y (2) e-” y (==*} . 
v _ = (5) 
dx Y . a 


T 


vields 


a : 
where € = ( ) and y= x: since the vortex strength 
x 


dx 
is linear in x. Eq. (3) thus requires that the resultant 
flow velocity normal to the concentrated vortex fila- 


+ Note that cosines of small angles are taken equal to unity 
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ment be a function of the vortex position and the wing 
half apex angle «. The complex velocity v* is, however, 
produced by the normal component of the main stream 
plus the normal component of the velocity due to the 
disturbance potential g; or using the complex crossflow 
velocity (v + iw) 


v* = —Ve (ao/a) + (wv + Iw), oa, (4) 
mean value 


and therefore 


(v + 1w)o >, Ve[(200/a) — 1] (5) 


mean value 
The mean value noted in Eqs. (4) and (5) represents 
the average streaming velocity over the vortex position 
and is easily computed by subtracting from the complex 
velocity the velocity of the vortex and then taking the 
limit as o approaches ap. 

The potential solution.—We now look for a solution of 
Eq. (2) which satisfies the aforementioned boundary 
conditions. The unknowns in the problem are first the 
location and strength of the concentrated vortex and 
finally the lift on the plate. We now introduce the com- 
plex velocity potential, W(o, a) the most general solu- 
tion of Eq. (2). W(e, a) is composed of the sum ¢ + 1y 
where y is a function similar to a stream function but 
having no physical significance to the three dimensional 
streamlines of the flow. The dependence of the solu- 
tion on x lies only in the relation between x anda. The 
solution of interest is derived by conformal mapping 
from the flow past two symmetrically placed vortices 
of equal but opposite strength. Thus in the @ plane 
(see Fig. 2) 


1 6 — & 


1 


In - — 1Vaé (6) 
2r 6 + A 


W(é) = -1 


Transforming to the o plane by the substitution @ = 


V o* — a’ yields 


: .Y Vo? — a* — Vao,"-—a 
W(c) = —15— In : _ 
<7 Ver - at Va? -— @& 
tVa Vo? —a?® (7) 


Eq. (7) thus represents the flow normal to an im- 
permeable flat plate having two symmetrically placed 
vortices of equal but opposite strengths located at the 
positions op and —6p. 

The boundary condition that the flow leave the plate 
at the edges yields a relationship between y and Va 


2rVa 


| I 
= + (8) 
¥ V oo” — a’ V G&G”? — a? 


This equation is easily obtained in the @ plane by re- 
quiring the presence of a stagnation point at the 
origin. 

The final boundary condition to be applied [Eq. (8) ] 
requires the calculation of the mean velocity at the 
vortex location and hence it is necessary to consider the 
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effect of the complex function W(c) less the complex 
potential function of the right hand vortex, or 
Wi(c) = W(e) + In (o — ao) Q) 


27 


Differentiating Eq. (9) with respect to o and setting ¢ 
equal to a yields the conjugate of the complex velocity, 
or (v — tw), at the vortex location. Rearranging real 


and imaginary parts of Eq. (5) there is obtained 
(v — iw), = Ve [(25./a) — 1] (10 


and combining Eqs. (7), (8), (9), and (10) we obtain 


1y | oa 
2m L(ao? — a?) + V (a0? — a?) (G2? — a?) 
on oo PM 
\ ‘(ae — a*)(G"? — a’) 6° — a’ 


l a” . Go 
= Ve 2 — | (11 
2 a(ay” — a’) a 


Eq. (11) gives the means of finding the vortex location 
oy, for the real and imaginary parts give two equations 
in the unknown coordinates yo and z. Upon separation 
of Eq. (11) into its real and imaginary parts the simul- 
taneous equations were found to be greatly expanded 
and not amenable to analytic solution. Consequently 
the equations were solved in a numerical manner, 
choosing a value of z) and finding the value yo which 
would give an equal value of the common parameter 
(y/Ve) for the two equations. The results of the 
calculations for vortex position are shown in Fig. 3, 
together with the corresponding value of (@/e) which is 
obtained by use of Eq. (8). It can be seen that as the 
angle of attack increases the vortex center moves in- 
board and up as would be expected from physical con- 
siderations. 

Pressure distribution and lift-—The expression for the 
pressure coefficient is to first order 


‘ ¢y" ’ 
- + gg? — | (12) 
V? z2=0 


on the wing surface. The second term on the right hand 
side does not normally arise when coordinates fixed to 
the free stream are used, however, it is necessary when 
the coordinates are tilted through an angle of attack as 
It does not con- 


Ap 7 | - 9 
q 


= 1/6 


is the case in the present analysis. 
tribute to the lift. The expression for ¢, is obtained as 
(da/dx)(O0¢g/O0a) or in terms of the complex function 
W(a) 


Ap 2€ (2) , 1 Fal 
Ss 2 +a’-— = (13) 
q J Oa /,=0 V*1 66 1], 6 


Calculations were made for a 15-degreet half apex 
angle wing for several values of a/e and the results are 
plotted in Fig. 4. The very low pressure region on the 

+ Other aspect ratios can be easily obtained since for any given 
value of (@/e) both pressure and lift are proportional to & 
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wing upper surface is caused by the presence of the vor- 
tex and the negative pressure peak corresponds approxi- 
mately to the lateral position of the vortex. 

In Fig. 5 are presented the corresponding span load- 
ings which are of interest in that quite a large deviation 
from the elliptical loading can be seen. 

The lift versus angle of attack curve produced by the 
separated wing is shown in Fig. 6; also shown for com- 
parison is the nonseparated case given by Jones.* The 
lift was obtained by integration of the pressure distribu- 
tion and also by computing the apparent mass of the 
system since the lift by Munk’s airship theory® is 


L = aV*m!' 


where !' is the apparent mass of two-dimensional flow 
pattern at the trailing edge. It is interesting to note 
that Jones’ result is obtained by using the irrotational 
flow past a flat plate whereas the present result is ob- 
tained using a special rotational flow pattern; since by 
Kelvin’s theorem the irrotational flow yields the least 
kinetic energy and hence least apparent mass, the lift in 
the present case must always be greater than that com- 
puted by Jones. 


COMPARISON WITH EXPERIMENT 


On Fig. 6 is shown some data obtained at a Mach 
Number of 1.9 in a blow-down jet of the NACA Langley 
Gas Dynamics Laboratory. The wing had a sharp one 
degree wedge airfoil section and a half-apex angle of 
10°. It was a semi-span model mounted on a boundary 
layer scoop-off plate and the forces were recorded on a 
strain gage balance. The data exhibit the same non- 
linearity as predicted but fall somewhat below the 
theoretical values. It is strongly suspected that com- 
pressibility effects at the high Mach Number of the 
tests have modified the flow pattern to a considerable 
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degree. To clarify this latter point, it can be seen from 
Fig. 4, that at Mach Numbers near 2 and moderate 
angles of attack, the pressure on the upper surface 
theoretically passes through absolute zero. Clearly 
then, important modifications of the flow must take 
place before these angles of attack are reached. ‘Tests 
at lower Mach Numbers would have allowed a better 
comparison; however, no such data were found at the 
time of this writing. 

Drag due to lift results are not presented for these 
flat wings as the drag equals the lift times the angle of 
attack, since there is no leading-edge suction. The 
lack of 
directly from the application of the Kutta condition on 


leading-edge suction in the theory results 


the leading edges; hence the edge velocities and pres 


sures are required to be finite. The center of pressure is 


located at the */; root chord position, as would be ex- 
pected for any conical flow. 


CONCLUDING REMARKS 


It has been shown analytically that leading-edge 
separation on slender delta wings leads to nonlinear 
lift curves with lift greater than that predicted by po- 
tential flow theory. For the flow model used in the 
analysis a concentrated vortex and feeding lattice was 
substituted for the more realistic spiral vortex sheet, 
thereby simplifying the problem to a great extent. The 
results should yield results in reasonable agreement with 
the real flows at least at moderate Mach Numbers and 
for slender enough wings. Comparison with experi- 
ments at a Mach Number of 1.9 indicated a qualitative 
check for the theory using a delta wing with SO degrees 
of sweepback. It is hoped that subsequently the prob 
lem of the spiral vortex sheet can be solved to afford a 
check on the approximations used in the analysis. 

No considerations of the stability of the vortex system 
have been given since it would appear necessary to drop 
the assumption of conical flow to test the reaction of the 
system to small disturbances from equilibrium. It is 
certainly to be expected that the real flow would de- 
generate to one producing a Karman vortex street at 
sufficiently high angles of attack. 

In the presented analysis the problem was simplified 
by a knowledge of the separation point on the wing and 
the conical nature o. the flow. In applying the same 
technique to other problems such as the separated flow 
over bodies of revolution a knowledge of the point of 
separation can only come from a boundary-layer analy- 
sis, and to further complicate the matter a step-by-step 
analysis must be made progressing from the nose rear- 
ward. Such a computation for a circular cone might 
prove to be possible with a reasonable amount of work 
since the flow field may still be considered conical if 
laminar flow exists. The solution would certainly be of 
basic interest; however, a solution for the general case 
would appear to require a formidable amount of labor. 


(Continued on page 706) 
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Simplified Analysis of 


Turbulent Boundary- 


Laver Dev deen Along Cylindrical 


Surf faces 


s with Variable Frec. Stream 


Mach Number’ 


HANS U. ECKERT? 
dir Research and Development Command, Wright dir Development Center 


(1) INTRODUCTION 


I ANALYSES OF BOUNDARY-LAYER FLOW along bodies 
of revolution placed with their axis parallel to the 
the thickness of the boundary layer is gener- 
body 


stream, 
ally considered as small in comparison to the 
diameter at any cross section, and terms proportional 
to this ratio are disregarded in the equations. This 
procedure seems well justified at the forward portion 
of the body where the diameter is increasing and the 
boundary layer is spread thin by flow divergence. It 
may become less accurate if the body has an extended 
cylindrical section, and it is probably no longer appli 
cable in case of a boat tail where flow convergence pro- 
duces an excessive thickening of the boundary layer. 
In the following analysis, which is based upon von 
Karman’s integral form of the boundary-layer momen- 
terms originating from the lateral 


tum equation, the 
Since the 


curvature of the surface 
equations become fairly complicated for bodies with 
variable cross section, the treatment will be restricted 
to a constant radius of curvature, that is, flow along 
This case has been treated in 


will be retained. 


the surface of a cylinder. 
simplified manner for a compressible turbulent bound- 
ary layer with uniform free-stream Mach Number in 
an earlier paper by the author.' In the present paper 
the Mach Number is considered as variable. For 
slender bodies in free flight there is, of course, little vari- 
ation in the potential flow Mach Number. A consider- 
able Mach Number gradient can be produced, however, 
along a cylindrical rod used as test object in a wind 
tunnel by variation of the test chamber cross-section 
area; as it has been done, for instance, by Chapman and 
Kester.?? 

Aside from this external flow problem the generalized 
analysis of the paper opens access to treatment of the 
internal flow problem, namely turbulent boundary- 
layer development in a cylindrical channel, which is 
inherently connected with a variation of the Mach 


Number of the potential flow core, due to the displace- 
Received November 27, 1953. Revised and received Febru- 
ary 16, 1954. 
The help of Mrs. A. Matonis and Mr. Kk. Cramer in prepa- 
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ment effect of the boundary layer. The combined 
case of flow in a channel with annular cross sections can 
also be approached. Two-dimensional flow is included 
as the special case of vanishing ratio of boundary-layer 
thickness to surface radius. 

While for this case 
boundary layers, approximate theories have been de- 
veloped by Walz? Tucker,*®: ® Ferrari,’ and Wilson,® 
only one comparatively early treatment is known for 
cylindrical pipes, which is due to Young and 
Since no allowance is made for density 


with compressible turbulent 


flow in 
Winterbottom.® 
variations in the boundary layer and for the influence 
of Mach Number upon surface friction, applicability 
of that method is restricted to low Mach Numbers. 

No method is presently known for treating external 
boundary-layer flow along a slender cylinder with vari- 
able free-stream conditions. f 

The analysis makes use of empirical relations, which 
are known to be approximately valid between time aver- 
aged quantities of turbulent boundary layers. No at- 
tempt is made to go into details of turbulent flow mech- 
anism or to consider effects of surface roughness, free- 
stream turbulence, heat transfer, or Prandtl Number. 
The surface temperature is taken as equal to the total 
Regions where flow separation 


temperature of the gas. 
strong shock 


is imminent, as in the neighborhood 


waves, are excluded from the treatment. 


(II) NOTATION 


= velocity inside boundary layer parallel to the sur- 


u = 
face 

U = velocity of the free stream just outside boundary 
laver 

Ma = Mach Number of the free stream 

Y = ratio of specific heats for air = 1.400 

P = static pressure 

p = density 

m = dynamic viscosity 

v u/p = kinematic viscosity 

T = absolute temperature 

x = distance from leading edge along surface 

y = normal distance from surface 

6 = boundary-layer thickness 


t For laminar incompressible flow this case has been treated 


recently in a paper by Cooper and Tulin 
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1/n = average exponent of velocity distribution power 
law 
r = radius of lateral surface curvature 
6* = displacement thickness of boundary layer in two 


dimensional flow 
6 momentum thickness of boundary laver in two- 


dimensional flow 


T shear stress at the surface 

A 6*/6 

Q) 6/6 

rr) function defined by Eq. (19) 

u* Vv ‘i p = shear stress velocity 

Cc = proportionality factor in ‘‘Law of the Wall,’ Eq. 
(21) 

1/n’ = value of velocity profile exponent as determined 
by wall friction 

a C~l22'/'+D) = proportionality factor in gener 
alized Blasius Law, Eq. (24 

Ris U6/v = Reynolds Number 

B (0 6)? na! —n)I/[n(n’ + 1] = velocity gradient coeffi 


cient in general shear stress relation, Eq. (80) 


ti 27/(p, U*) = coeflicient of local shear stress 


i 
Cr ( dx 
x / 
J 0 


factor defined by Eq. (81) 


= coeflicient of mean shear stress 


€ = 
f = function defined by Eq. (40 
g = function defined by Eq. (41) 
? = function defined by Eq. (44 
r 87/(pu*?) = local coefficient of friction in pipes 
d = inside diameter of pipe 
u = average velocity across pipe section 
Re (tid)/vy = Reynolds Number of pipe flow 
Subscripts 


i denotes incompressible flow conditions 

U’ denotes conditions at edge of boundary layer 
n denotes particular velocity profile 

0 denotes stagnation conditions 


(III) ANALYSIS 


(a) Momentum Equation 
For steady boundary-layer flow along a cylindrical 


surface von Karman'’s momentum theorum reads 


d 


dx « 


d 


dx « 


r pur(r + y)dy —U pu(r + y)dy + 
y y 0 


V o 
dp 
(7+ ydy+rr=0 (1) 
Jy 0 dx 

In case of internal flow (flow in a cylindrical pipe) 7 
has a negative value. 

If the variation of static pressure across the boundary 
layer is disregarded, as it is common in treatment of 
nonseparating boundary layers, the third term reduces 


to 
dpf[. 6° 
dx (va ” 5) 


and Eq. (1) can be written after rearrangement 
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dU [°° a ie 
pu dy — pu(U’ — uj)dy + 
J 0 dx JO 


dx 
dp  1{dU [ j 
6 + uydy — 
dx r Ldx Jo oe 9 
d [ U F 6° dp 0 : 
1 ( —u)yday + + = (2) 
dx Jo ‘ ie 2 dx id 


From the first two integrals of Eq. (2) follow the con- 
ventional definitions of momentum thickness and dis- 
placement thickness for a two-dimensional boundary 
layer. If we assume for the velocity profile a power 
law of the form 


u/U = (y/6)'/" (3) 
as it appears to be realistic according to available ex- 
perimental data,'?~'* ® then we can in a simple manner 
express also the third and fourth integral by two-dimen- 
sional terms for momentum thickness and displacement 
thickness. The only difference is that these terms have 
to be evaluated from a velocity profile with half the ex- 
ponent of the basic profile. Indicating the appropriate 
profile parameter by a suffix the following relations 


can be established 


pu(U — u)dy = p, U6, (4) 


26 
| pudy = p,l'(6 — 6,*) (5) 


76 
{ pu(U — u)ydy = — py, Ube, (6) 


6 
l : 
i} puydy = 5 py V6 — b2n*) (7) 
0 


Though relations similar to Eqs. (6) and (7) have 
been derived and applied in two earlier papers, '* their 
derivation is repeated for convenience in Appen- 
dix 1. 

By introducing Eqs. (4) to (7) into Eq. (2) and fol- 
lowing the common procedure to eliminate the pressure 


gradient by Euler’s equation 


dp _ yy tl (8) 
dx ial dx 
we obtain 
1U 1 ; 
Pl US i + (pyl 0.) : i 
dx dx 
1 1U 1 
ee U: 552,* + (py U750en) = 7 (9) 
2r dx dx 


After carrying out the differentiations, rearrangement of 
terms and division by p,U” Eq. (9) is changed as fol- 


lows: 
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16, 1dU 1 dp 
eo (26, + 6,*) + "he 
dx U dx py, dx 
6 | dbs, A, a6 1 dU 
alae : (26>, + bon*) —_ 
2r 1. dx 6 dx U’ dx 
ld T 
~Pe A, | = : (10) 
pu dx py U? 


We can express variations of velocity and density in 
the free stream by variation of the free-stream Mach 


Number according to the relations 


] dl 5 l dla i 1] 
Udx Ma dx T, 
l j ( 1 ‘a 
dp, —* dMa : (12) 
py dx qn To 


where 


and obtain from Eq. (10) 


6, dMaT, (* + 6,* ) 
— + = — \Jaé,} + 
dx dx Ty) Ma 
6 j 162, Oo, d6 dMaT, 


2r Idx J 6 dx = dx 7 


260, + 62,* 
| ee i a te 
Ma f py? 


If Eq. (13) is rearranged and written in the following 
way, 


d ( 6 ) 1 dMaT, 
: 6, + Ao, + 
dx 2r Ma dx Tp) 


r ee, — ' 6 
26, + 6,*% + — (260, + 62,,*) — Ma*| 6, + Ao, = 
2r > 


T 


Pr U? 
it becomes evident that for values of the free-stream 
Mach Number 


5 
26, + 5,* + — (26, + bm,*) 


Ma? = = (14) 


the term with diV/a/dx as factor vanishes. At these 
particular values of the Mach Number the effects upon 
the growth of the momentum thickness produced by 
simultaneous changes of velocity and density cancel 
each other. At higher Mach Numbers the density 
changes predominate, so the resulting effect of a pres- 
momentum thickness will 


sure gradient upon the 


be opposite to that known for low speeds. 

From Eq. (14) follows that the critical value of the 
Mach Number depends upon the parameter of the 
velocity profile and the parameter of lateral curva- 


ture 6/r. The equation has been solved graphically 
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Mach Number values obtained are shown 


They are plotted versus n for values of 6 r 


and the 
in Fig. 1. 
equal to 1.0, 0.0, and —1.0 representing a boundary 
layer at a convex, a plane, and a concave surface, re- 
spectively. It can be easily shown that for large values 
of m the solution becomes independent of n and 6/r 


and approaches a value of V3/(2 — ¥ 
1.400 is 2.236. 


In order to keep the analysis within reasonable limits 


which for 4 


n will in the following be considered as a constant. 
This will restrict, of course, the applicability to rela 
tively small variations of the free-stream Mach Num 
ber, particularly in cases of flow deceleration. Since, 
however, at higher Mach Numbers the shape of the 
velocity profile becomes less significant, due to the 
dominating role of the density changes, this assump- 
tion is not as crude as it would be for low speed flow. 
Moreover, the pertinent functions will be evaluated 
for different n-values so that variations of this param 
eter in the course of the boundary-layer development 
can be taken into account to some extent by changing 
it intervalwise. 
We now introduce the nondimensional terms 


Qi), = Oe),/5 and Aw), * (9,4 


which for a constant ” depend on x only indirectly over 
the variation of the stream Mach Number with the 
flow length MJa(x). The boundary-layer thickness, 
on the other hand, may be directly dependent on x. 


We can express this by 


62), = Qe), [| Ma(x) | 6[.Ma(x); x] (15) 
5.2)” Aw), {| Ma(x)] 6[Ma(x); x] (16 
It follows that 
6.2), dé dQ»), dMa = 
=: () + 6 (14 
dx dx dMa dx 


Introducing Eqs. (15), (16), and (17) into Eq. (13) 


we obtain 


dé | dQ, Ty (= + A, ) | 
: + 6 + VWaQ, 
' dx dx | dMa T\ Ma 


Q 


0 j dé dMa Fa r 
2Q, rs 
29r ts dx dx (tdMa 
J ( = — MaQ I 4 = is 
te Ma py WV? 


or, after rearrangement and if setting for briefness 


dQ»), en (7 et Ag, 
Oe” die * Th Ma 


dé ( 6 )) dMa 
2, +- Qe) )} + 6 x 
dx r dx 


oT al = ~  \sY) 
2r . pul F 


Ma%x) (19) 
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b) Shear Stress Relation 


In order to obtain from Eq. (20) a detailed and real- 
istic picture of the turbulent boundary-layer develop- 
ment in flows with variable Mach Number we have to 
find a relation which connects the shear stress at the 
wall to local flow conditions and suitable boundary- 
layer parameters including thereby, at least approxi- 
mately, also the effects of Mach Number and Mach 
Number gradient. Since there is, at the present time, 
not even for low speeds an appropriate theory available 
such a relation must necessarily contain empirical ele- 
ments. For lack of suitable experiments with turbu- 
lent boundary layers under pressure gradients in com- 
pressible flow we have to derive the basic relation from 
incompressible flow data and allow for compressibility 
by a separate correction. 

(1) Relation for Incompressible Flow. 
ments with accelerated and retarded two-dimensional 
flow at low subsonic speeds Wieghard,'* Ludwieg and 
Tillmann,’ and independently Ross and Robertson'® 
observed that the influence of the pressure gradient upon 
the shape of the velocity profile is restricted to the 


From experi- 


outer portion of the boundary layer, while a region 
near the wall remains relatively unaffected. 

According to reference 15 the variation of velocity 
with distance from the wall can, for this region, be well 
approximated by the power law relation 


u C (" ) - 
u* y 


where u* = V 1r/p is the shear stress velocity based on 
the wall shear stress and Cand n’ are functions of Reyn- 


(21) 


olds Number. 
From Eq. (21) follows for the shear stress coefficient 


2n’/(n’'+1) TaN (2/(n’+1)] 
aT tt kale | =) : *) (22) 
pl 2 l v 


Eq. (22) indicates that r/pU? will be raised or lowered 
in the same sense as the velocity «7 at a certain wall dis- 
tance y is changed. In cases where the profile is not 
disturbed by sudden changes of the border velocity U 
as along flat plates and also in fully developed (in- 
compressible) pipe flow, the velocity distribution as 
determined by the ‘law of the wall’’ Eq. (21) extends 
practically to the free stream and the entire profile can 
be represented by 


u/U = (y/6)'”" (23) 
In such cases in Eq. (22) may be set wu = U and y = 6 


which reduces it to 


= — [2/(n’+1)] 
To Co n/n" 2G *) (24) 


pl y 


With n’ = 7 and the empirically determined value for 
the proportionality factor Eq. (24) further specializes 
to the old Blasius formula. 
or deceleration of the outer flow the experimental data 


For cases of acceleration 
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of reference 15 indicate that Eq. (21) is still applicable 
up to a distance from the wall, which roughly equals the 
momentum thickness where «4 = w,’.. Hence in this 
general case Eq. (22) may be written 


9 


As it was stated before by Eq. (3) the profile will, 
also in the general case of nonuniform border velocity, 
be approximated by a power law with the average ex- 
ponent 1/n. x will be so determined that for a given 
6 the ratios 6;/6 are the same for the actual profile and 
the power law approximation, that is 


n _ 6, en u ; MY 
(n+1)(n+2) 6 ° J, U uj”? 


or 
1 — 39, 1 — 3a,\2 . 
n= - + —2 (26) 
20; 22; 
From Eq. (3) follows for the velocity at y = 6 
lin 

ie 0 , 

“= (3a) 
l 6 


If only moderate pressure gradients are considered, 
the difference between wu,’ and wu, will not be great and 
no excessive error will be caused by setting up,’ = 
ue;, Under this condition Eq. (25) can be written 


[2(n’ n)]/[n(n’ + 1) 
T i a 6; ; , 
T “aaidailie ws ) x 
cU* 6 


Us [2/(n" 1) 
°) (27) 
v 
or with 
a= Ca ee’ n 1)] 
n [2(n’ n)]/[n(n’ + 1)] 
8 = 
(n + 1) (m + 2) 
and 
U6 
Re; = ‘ 
v 
T , , — [2/(n’+ 
= = a(n’) B(n'; n) Re, satiated (27a) 
pU? 


Since ’ is a function of Re; Eq. (27) expresses the 
wall shear stress in terms of Reynolds Number and a 
profile parameter. Similar relations for plane incom- 
pressible flows with pressure gradients have been de- 
veloped by Ludwieg and Tillmann,' Ross and Robert- 
son,'® Granville, '’ Rotta.’ In 
formulas the ratio 6;*/@; = H; is used as profile param- 


and most of these 


eter and 6; as characteristical length for the Reynolds 


Number. In the present paper which is dealing with 
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compressible flows and curved surfaces the use of n and 
§appears more suitable. 

Values for the exponent in wall proximity ’ and the 
coeflicient of Reynolds Number a@ are plotted against 
log Res in Fig. 2. They have been derived from the 
well known von Karman-Nikuradse relation for the fric- 
tion coefficient in fully developed pipe flow*’ and experi- 
mental data from references 24 through 29 as it is shown 
in Appendix 2. 

Values of the velocity gradient factor 8 as a function 
of n with n’ as parameter are given in Fig. 3. 

In Fig. 4 values of 7/pl™” are plotted against 7 with 
Re; as parameter according to Eq. (27a). Values of 
a and n’ for a particular value of Re; have been taken 
from Fig. 2. The dotted lines are obtained from the 
interpolation formula of Ludwieg and Tillmann" which 
closely represents their own test results covering a 
range of Re; between approximately 10* and 3 X 10°. 
For a power law profile this formula reads in terms of 
dand n 
0.268 

0.123 
— = Re; 


= 0.678 [(n +2)/n 


cU? 10 


The dash-dot lines represent values given by Rotta."® 


0.268 


They are derived by basically similar, but more refined 
considerations which include application of an energy 
theorem for turbulent boundary layers. 

One notices from Fig. 4 that, though differences exist 
in some regions, the three methods yield essentially 
the same result, namely, increasing 7/pl* with in- 
creasing m-values and decreasing 7 pl° with increasing 
Re;. 

Differences between experimental results are still 
larger, particularly in the region of low n-values which 
corresponds to adverse pressure gradients. The hot 
wire tests by Schubauer and Klebanoff"’ lead to values 
for +/pU? which are about forty per cent higher than 
those of Ludwieg and Tillmann.’” Data obtained 
more recently by Sandborn” with a heat-transfer instru- 
ment similar to that of reference 15 agree within twenty 
per cent with the data of that reference. Further in- 
vestigations appear necessary to narrow the present 
range of uncertainty. t 

+ After completion of this paper the author learned that 
Clauser*! recently has conducted tests with retarded boundary 
layers in which he carefully tried to avoid the possible sources 
of error in previous investigations. He derived a relation for 
the effect of pressure gradient upon the friction coefficient by 
evaluating an auxiliary function in the general logarithmic skin 
friction law with the aid of his test data. Values computed 
from this relation and converted to our parameter system have 
been inserted as dash-double dot lines into Fig. 4. Consider- 
ing the entire range of Reynolds Number the agreement is par 
ticularly good with the results of the present analysis and adds 
weight to them if not to the analysis itself 

The difficulty to verify strictly two-dimensional boundary- 
layer flow experienced in the above tests emphasizes the im- 
portance of a theory correlating axisymmetric with two-dimen- 


sional flows. 
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(2) Compressthility Factor—The influence of com 
pressibility upon the friction coefficient has been de- 
termined with reasonable accuracy from a number of 
experiments with turbulent boundary layers in con- 
stant pressure flows up to a Mach Number of about 
five. If Reynolds Number is based upon the free- 
stream conditions LU’, uy, pr, then it turns out that at 
this Mach Number the friction coefficient is less than 
half the value for incompressible flow at the same 
Reynolds Number. Since values for the ratio of fric 
tion coefficients obtained by tests of this kind reflect 
only the influence of Mach Number upon the momen- 
tum flow in the boundary layer, they may not be 
applicable in cases where a pressure gradient largely 
affects the skin friction. It is not known how to com- 
bine accurately the effects of compressibility and pres- 
At the present state of knowledge it 
with the 


sure gradient. 


appears the best to multiply Eq. (27) 


ratio 


T, pr l 2 Cr 
(7/pU) in Cy 
flows of uniform 


tests with pres- 


as obtained from 


sure. 

The data of Coles*! which have been obtained from 
direct measurements of the friction force upon a float 
ing element and which have been corroborated re- 
cently by Chapman and Kester** appear to be the most 
reliable at the present time. They are shown in Fig. 5. 
Since these c¢,/c;; values refer to Reynolds Numbers 
based on x rather than on 6 they might not be directly 
applicabie unless the ratios 6 x are the same for com 
pressible and incompressible flows. We can check 


this in the following way 


With 


where c, denotes the mean value of cy, on a flat plate 
with constant Mach Number and constant velocity 
profile, we obtain for the thickness ratio of a com- 
pressible to an incompressible boundary layer at the 


same position x 


6 €76;6 Cr, 
— = 29 


The product cy, ¢;,, 2; 2 with c, c,, taken from reference 
2landQ, Qforn = 5,7, and 9 from reference 5 is plotted 
versus Mach Number in Fig. 6. The values are close 
to one and show little influence of the profile parameter. 
This result indicates that at least up to Mach Number 
5, the thickness of a turbulent boundary layer on a flat 
plate can with good approximation be considered as 
independent of Mach Number. There are also some 
data scattered in the literature to support this con 
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Fic. 1. 
the influence of a pressure gradient upon momentum defect in 
Velocity profile parameter 


the boundary layer vanishes. 27: 
6/r: Parameter of lateral curvature. 























Fic. 2. Velocity profile parameter as determined by wall 
friction n’ and coefficient of generalized Blasius law a in depend- 
ence of Reynolds number Res. 
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Fic. 3. Variation of velocity gradient coefficient 8 with param- 
eter of outer part of the velocity profile » for different values of 
parameter of inner part of velocity profile n’. 
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Fic. 4. Variation of shear stress coefficient 7/p U2, with veloc- 
ity profile parameter » and Reynolds Number Res for incom 
pressible flow. Comparison of results of different investigators 
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Fic. 5. Effect of Mach Number on skin-friction ratio for 
flat plates and cylinders in axial flow 
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+ Measurements of O x/2r=8 
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friction force by 
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Ux/vu = 107 





clus 
enc 
stre 


gral 


Ur) 
pap 
of E 
dric 
Re;, 
tain 
wit! 
A su 
wall 
shea 
abot 
wall 
char 
tota 
two- 
It 
the 
gene 
Eq. 
data 
Kart 
an il 
part 
cave 
W 
clud: 
corre 
b/r < 
It 
pres: 


later 


e 
the s' 
Mach 
and r 
Incotli 
crease 
ducti 
only { 
this c 
in ten 
confo: 
and J 
effect 
flow 
functi 
proac 


also t 





0° 





/ 

a 
——| 
—| 
— 
an} 
a 
— 
ne 








- 


reloc- 
com 
tors 


kert 


for 





ANALYSIS OF TURBULENT BO 


We thus can combine the values of refer- 
to obtain an approximate shear 


clusion. f 
ence 21 with Eq. (27) 
stress formula for compressible flows with pressure 
gradient as follows 

T ) 
Re, * 30 


cy 
- = a(n’) B(n'sn (Ja 
pol “i Cy 


No relation can be given for a possible effect of m upon 
There are indications, however, that such an 


Ce Cfz5 
effect will be small. 

3) Sensitivity of the Shear Stress Formula to Lateral 
Surface Curvature.—In the of this 
paper, consideration must be given to the applicability 
of Eq. (30) to external and internal flows along cylin- 


line with purpose 


drical surfaces, that is, whether for given values of 
Re;, n, and Ja the same value for 7/ pyl will be ob- 
tained irrespective of the variation of the ratio 6/7 
within certain ranges of positive and negative values. 
A suggestion of this kind is favored by the fact that the 
wall layer to which Eq. (21), the basic relation for our 
shear stress formula, applies has a thickness of only 
about ten per cent of the total boundary layer. The 
wall layer will therefore retain its two-dimensional 
character also in cases of surface curvature where the 
total boundary layer can no longer be considered as 
two-dimensional. 

It is also stated by Ludwieg and Tillmann' that 
the “law of the wall’ is valid for flow in pipes. The 
general agreement between the formula of these authors 
Eq. (28) which is based upon two-dimensional flow 
data and our Eq. (27) based primarily upon the von 
Karman-Nikuradse relation for pipe flow may serve as 
an indication for the insensitivity of the incompressible 
part of Eq. (30) to lateral surface curvature in the con- 
—I. 

With respect to the above reasoning one may con- 
clude that the same insensitivity will exist also in the 


cave range that is, for 0 é6/r S 


A 


corresponding range of convex surfaces that is, for 0 
b/r < 1. 

It is difficult to predict accurately how the com- 
pressibility factor will behave under the influence of 
lateral surface curvature. If we assume that the con- 

t Pitot tube surveys from various sources indicate that also 
the shape of the velocity profile is not noticeably changed by 
Mach Numbers within the considered range. With profile shape 
and rate of increase of the full thickness being the same as for an 
incompressible turbulent boundary layer the reduced rate of in- 
crease of the momentum thickness which equals the observed re 
duction of the friction coefficient by Mach Number can result 
only from a reduction of density. For a constant pressure layer 
this density reduction in turn can be produced only by the rise 
From this it appears, in 
‘and by Young 


in temperature due to frictional heating 
conformity with suggestions made by Monaghan 
and Janssen,* that in order to account for the ‘‘compressibility”’ 
effect upon turbulent skin-friction modification of incompressible 
flow formulas by a ratio rather than a 
function of Mach Number is the logical and more general ap 
proach; more general, because it includes, at least approximately, 
also the cases of heat transfer to and from the surface 


suitable temperature 
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layer thickness for various velocity power law profiles 


dition 6/62 = 1, which as we have seen, appears to 
hold reasonably well for flat plate flow can be extended 
also to curved surfaces then we can make a simple esti 
mation of this influence from the differences in momen- 
tum flow. 

The relation between the compressibility factor of a 
cylindrical surface and that of a flat plate may be given 


by 


Cf a a Q) 6 
Cri Cy Q; 6 
is equally valid for cylinder and plate and with 6/6, = 1 
we can write Eq. (31) 
Cyl Plate 


q) ) 
_— (32) 
Q Q 


For a profile with the constant exponent 1/ is accord- 
ing to Eq. (1.6) of Appendix | 
Cyl Plate Plate 

2, = 2, + 


(6/27)Qp, (33) 


Under this condition Eq. (32) can be written 


27) Qe, Q, 


= ¢ (34) 
2r) Qs, Q ni 


«) + (6 


8) + (6 


where the superscript “‘plate’’ has been omitted for sim- 


plicity. Eq. (34) yields 
1 + [(6/2r) (Q,/2,) ] 


«= 

1+ [(6/2r) (Q2,;/2,;)] 
The ratio 2,/Q, Mach Number for 
5, 7, and 9 in Fig. 1 of reference 1. 


is plotted versus 
n= 

In order to facilitate comparison with experimental 
results of reference 22 we correlate « to ratios of mean 


friction coefficients. 
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With 
Cy cf 
Cri a 


valid for cylinder and plate follows from Eq. (31) 


Cyl Cyl. 
re (36) 
Plate © Plate : 
Cy Cri 
Cyl. Plate 


A relation for the ratio c,;/ c; has been derived as 
Eq. (14) in reference 1 for the special case of the Blasius 
shear stress law. For incompressible flow this equa- 


tion reads: 


6 Yon; 
Cyl ] 
; 9 
Cri ZY Ui, ies 
Plate ~~ ane 1/5 (a4) 
‘ DO Sloan; 
cr 1 + 
Ir Qa; 
Cyl. Plate 


Eliminating in Eq. (36) cy, / ¢,, by Eq. (37) and e by Eq. 


(35) yields 


0 Qo, 
Cyl ] + ; 
= Dr CO 
es il. (38) 
Plate ~ 5 6 ©. 15 
Cy | 4 se2y 
or Q, 
Cyl. Plate 


This alternate expression for c; / c; based on the as- 
sumption that there is no influence of Mach Number 
upon the thickness of a turbulent boundary layer in 
constant pressure flows predicts for compressible speeds 
a higher effect of 6/r upon the mean friction coefficient 
and is in better agreement with the experimental 
findings of reference 22 than Eq. (14) of reference 
be 

Considering, however, that no quantitative agree- 
ment has been established yet and that for the majority 
of cases the effect is not excessive compared to errors 
introduced by other limitations of the analysis it is 
not felt justified to complicate the evaluation by this 
refinement. We will therefore in the following con- 
sider € as equal to one and Eq. (30) as independent of 


lateral curvature within —1 < (6/r) @ + 1. 


n’ 


a +3 i 5\ 5 
- gi x: e” 
‘ i or A: r) 


n 


3 
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<2 | ie 6 F 
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(c) Approximate Solution of the Momentum Equation 


Combination of Eqs. (20) and (30) yields 


dé 6 _adMa 6 
v2. + tk, | + 0 Cn T P2n = 
dx r dx 2r 
” I r [2/(n’ + 1)] ; 
ap cf ( ) ge FS (239 
Cy; Vr 


No integration of Eq. (39) in closed form appears feas- 
ible except for the two-dimensional case 6/r = 0 where 
it reduces to a differential equation of the Bernoulli 
type. However, if the ratio 6/r can be considered as 
a parameter and kept constant within small steps of 
x, Eq. (39) reduces to the same type of equation as for 
the two-dimensional case. For conciseness we intro- 
duce the functions 


dMa 6 o>, + (6/2r) go, dMa 
f (110; -- -° n) =< ope (40) 
: dx =r Q, + (6/r)Qs, dx 
and 
¢ fU\-* 
ap ) 
= as Ba. (41) 


The function f depends on x, strictly speaking, by 
way of the functions \/a(x), n(x) and [(6/r) (x)]. As 
mentioned in section IIIa the analysis will consider 
only such cases where ” does not change with x over 
longer intervals. If also 6/r is kept constant within 
small intervals of x function f depends upon x only over 
Ma(x) 
Within such an interval f may therefore be considered 


the function which is assumed to be known.t 
as a known function. 

The same consideration applies to the function g. 
The additional variable n’ is a function of Reynolds 
Number which in turn depends on 6, Mach Number 
and the stagnation conditions of the flow py and 7». 
For isentropic free stream assumed in this analysis 
they do not change with x. We can therefore regard 
also g as a function known within steps and Eq. (39) 
may be written 


w+)) —@ (42) 


> 
= 

> 

| 

79 
= 

mo 
' 


For a step x, 41 — X, an approximate solution of 


Eq. (42) can be established by 


3 Penis [> (?) Ja i + , 


- +e” (43) 


t In the cases where the displacement effect of the boundary layer upon the potential flow cannot be disregarded as in pipes 
and channels, the variation of Mach Number and effective cross section can be brought into agreement by an iteration procedure. 
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ANALYSIS OF TURBULENT 
where 6, 4 , and 6, denote the boundary-layer thickness 
at positions x, 4 , and x,, respectively, and (6/r) is the 
mean value of 6/r between those positions. If the 
changes in 6/r are small the error will not be appre- 
ciable if we replace (6/r) ; by (6/7), as parameter 
in functions f and g. 

Since with ” and 6/r as constants the function f de- 
pends upon x only over the Mach Number variation 
Ma(x) and since it has as a factor d\/a/dx we obtain 
°XK *Ma(xk 
fax— bd Ma 
XK YJ Maixt 


« 


where 

6 

On T 2 
2?r 

?P,, is (44) 

6 

q) alee () 
> 


® is a function of Mach Number with the parameters 


é/r and n. The integrals in the exponents Eq. (48) 


can therefore once and for all be calculated and 
tabulated. 
If we linearize the Mach Number variation a(x) 


within intervals as it has been done for the two-dimen 
sional case by Walz‘ and Tucker’ already we can evalu- 
ate and tabulate also the main integral of Eq. (43) as 
a function of Mach Number with 6/r, ”, and n’ as 
parameters. 

Tables of these functions together with practical ap 
plications will be presented in a later paper. 

For a refined treatment of the boundary-layer de- 
velopment the analysis should include a relation which 
predicts the deviation of » from ’ under the effects 
of pressure gradient, lateral curvature and Mach Num- 
ber. At the present time there is, however, neither a 
sufficient theoretical nor an experimental basis from 
which to develop such a relation. 


(IV) SUMMARY 


An approximation method is derived from the mo- 
mentum theorem for calculating the thickness of a tur- 
bulent boundary layer developed along the convex or 
the concave side of a smooth, thermally insulated 
cylindrical shell placed into an isentropic stream of 
variable Mach Number with the axis parallel to the 
flow flow is included 
in the analysis. 
total temperature across the boundary layer are con- 


Two-dimensional 
It is assumed that static pressure and 


direction. 


stant and that the velocity distribution can be repre- 
sented by a power law. The exponent of the velocity 
profile can be chosen to suit a particular type of flow 
but is then considered as invariable with flow condi- 
tions and geometry. 

It is shown that some corrective terms arising from 
the curvature of the surface relative to the boundary- 
layer thickness can be evaluated from a two-dimensional 


BOU 


ELOPMENT 7 


NDARY-LAYER DEV 
boundary layer having a profile exponent which is one- 
half that of the actual profile. A relation for the wall 
shear stress in terms of profile exponent and boundary- 
layer thickness is developed from experiments with 
flow The 
influence of compressibility is taken into account by an 


incompressible under pressure gradient. 
empirical factor obtained from constant pressure flows 
Estimations of the sensitivity of the shear stress rela 
tion to lateral surface curvature lead to the conclusion 
that effects of curvature are small that in first 
approximation the relation may be considered as uni 
versally applicable to plane and curved surfaces. 
The differential equation for the boundary-layer 
growth resulting from combination of the momentum 
equation with the shear-stress formula can be integrated 


and 


in closed form only for the special case of two-dimen 
sional flow. An approximate integration is possible 
in the general case if the ratio of boundary-layer thick 
ness to cylinder radius is considered as parameter and 
kept constant within intervals. By combining this 
procedure with linearization of the potential flow for 
the same intervals and a one-time evaluation of the 
differential equation as function of Mach Number 
practical application of the method can be reduced to a 


routine operation. 


(V) APPENDIX 1—RELATION 3ETWEEN INTEGRAL 


PARAMETERS OF Two-DIMENSIONAL AND ROTATIONALLY 
SYMMETRICAL BOUNDARY LAYERS 


The momentum thickness of a boundary layer 


generated by axial flow along a cylindrical surface can 
be defined in analogy to the two-dimensional case from 


Cyl 
pul *al(r + 6)? — r*] an | pul u) (r+ y) dy 


wall curvature (internal flow) r is con- 


Disregarding the term 6° for being 


For concave 
sidered as negative. 
small of the second order and dividing by 2zrpyl™ one 


obtains from Eq. (1.1 
A pu ( , ul )( 14 *) ™ 9 
« Pr U'\ l r , 


setting for conciseness 


and splitting the integral yields 


Cyl 


Q } PV(1 ss V)dVY t 


If the velocity profile is given by 
V= Y/Y 


the density profile is accordingly ** 
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I . | eal ° 
P = - sa, Where K == Ma? 
1+ A(1l— Y~") 2 
provided that static pressure and total temperature 


do not change with Y. 
Eliminating V and P one obtains from Eq. (1.3 


Cyl 
Qs 
sen 
0 


yee a em YY 
i+ Kki-~ ys 
y'/*(1 y'/") VdV 


| - 
gJe 1+ Kii — Y*") 


Substituting Y = Z" yields from Eq. (1.5) 
_ if Z"(1 — Z)dZ 
01+ K(1 — 2?) 
é a8 | Z"(1 — Z)dZ “ss 
(1.9) 
01+ A(1 — 2?) 


(1.4) 


r 2 


Comparison of both integrals reveals their identical 
structure except that the second integral employs the 
quantity 2m instead of n. Since the first integral repre- 
sents the momentum thickness for a flat plate (related 
to the total thickness) the second one can also be inter- 
preted as a flat plate momentum thickness derived from 


a velocity profile 


V on y! 2n 
Eq. (1.5) can therefore be written 
Cyl Plate $ Plate 
0 . 
Q,, = Q,, + Qs, ( ] .O) 
2r 


and it follows from Eqs. (1.2) and (1.3) that 
: 
f pu(U — u)y dy = = pu U"602, 
0 _ 
which is Eq. (6) of Part IIIa. 


In the same way one obtains for the displacement 
. * - 
thickness 6 = Aé from 


[ evav= — A, 
0 

1 - 1 
[ evva = oil =~ de 


(1.7) 
that 


(1.8) 


or 


5 
l * 
f pu y dy = = puU 5(5 — bx ) 
0 ~ 
which is Eq. (7) of Part IIa. 


(VI) APPENDIX 2—EVALUATION OF THE SKIN FRICTION 
LAW FOR INCOMPRESSIBLE FLOW 


In order to make practical use of the shear stress rela- 
tion Eq. (27a) 


SCIENCES—OCTOBER, 1954 


- = a(n’) B(n'; n) Re, 
é 
we have to determine the relations between a, n’, and 
Re;. Since there is no adequate theoretical foundation 
we have to employ empirical results. The experiments 
on turbulent fluid friction which cover the widest range 
of Reynolds Numbers (about 2.5°10* < Re; < 1.9° 108 
appear to be those of Nikuradse.** 
His results are condensed in the relation 


1/VX = —0.8 + 2 log (ReV Xd) (2:4 


where \ = 87/(pu*) is the friction coefficient of fully 
developed flow in pipes and Re = (iid), v is the Reynolds 
Number based on the average velocity across the pipe 
section and its diameter.t For fully developed pipe 


flow isd = 26. Re and Re; are therefore related by 


u 


5) 
Re = Res. (2:3 


where for power law profiles 


u 2n? 


U (n + 1) (2n + 1) 


Since LU’ 
assume that the velocity profile develops undisturbedly 


is constant for incompressible flow we can 


over the entire boundary layer according to Eq. (21). 


It is then » = n’ and 8 = | so that Eq. (27a) reduces to 


(r/pU*) = a(n’)Re, 2 


Te os 

In order to give a survey on experimental findings 
values of n’ have been plotted against log Re; in Fig. 7 
as they were obtained from measurements in pipes by 
Nikuradse;** on flat plates by Dryden,** Schultz- 
Grunow,”~ Townsend,”® Klebanoff and Diehl;*? and in 
two-dimensional and circular channels by Laufer.*” ** 
The scattering is considerable. The flat plate data 
are generally lower than those for channel flow which 
indicates that a favorable pressure gradient (as present 
in subsonic channel flow) tends to raise the value of n’ 
also if no flow acceleration takes place. This is also 
made evident in Fig. 20 of reference 29. A more re- 
fined treatment will therefore require different relations 
between 7’ and Re; for flat plates and channels. With 
respect to the limitations of the present paper it appears 
permissible, however, to disregard such differences and 
consider the experimental data of Fig. 7 as represented 
by the relation 


n' = 2(log Re; — 1) (2.5) 
or 


Re, = 10!" + 2/21 (2.5a) 


6 

1 The change in the constants of Eq. (2.1) by a few per cent 

which has been suggested recently by Ross* in a reanalysis of 

Nikuradse’s data is within the range of uncertainty of the present 
analysis and has no relevant effect upon the final results. 
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ANALYSIS OF TURBULENT 
With the relation between n’ and Re; fixed we will 
now determine a(n’) in such a way that the von Karman- 
Nikuradse relation Eq. (2.1) and Eq. (2.4) which has 
the form of the Blasius law become equivalent ex- 
pressions. 
From combination of Eqs. (2.2) and (2.5) follows 


log Re = [(n' + 2) 2] + log 2 + 


log [(4/U) (n’)] (2.6) 


On the other hand we can write Eq. (2.1) 


l l 
0.4 — ( tog» - 
2 Vi 


Eliminating log Re from Eqs. (2.6) and (2.4a) yields 


(2.1a) 


log Re = 


u ; l l 
-(n') + 1.802 = log -+- . (: 


l Vd 


n’ + 2 log 
Eq. (2.7) has been evaluated graphically by making use 
4 to 11 is 
found to be representable with 0.5 per cent accuracy 


of Eq. (2.3) and between values of 1 from 


by the relation 


(1/V A) = 0.911 n’ + 0.600 (2.8)T 
With 
T Al 4a 2 
= - (# ) 
pl? 8 Ll 


and a combination of Eqs. (2.4), (2.5a), and (2.8) we 


finally obtain 
[((@/U) (n")fP- 10% + Oe +! 
7 8(0.911 n’ + 0.600)? 


a (2.9) 

It can be seen from the numerical values of a and n’ 
plotted in Fig. 2 that for n’ = 
the Blasius law 


7 Eq. (2.1) specializes to 


r/pU? = 0.0227 Re, ‘'"4 
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lysis of Institute of Fluid Mechanics, Technical University of Braunschweig, ' iT } \ 
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2 5 l T 
‘bout N A TWO-DIMENSIONAL CASCADE the local loss coefficient of the | c 
— ) , 30 4-be }—» Z Cc 
behead blade is constant along the length of the blade In this case 4) o 7 
the local loss coefficient ¢),.(z) is identical with the mean loss =) i l 
ment, : é : c A 3 
5-102 coefficient Cora. of the whole cascade.t For the cascade installed ® lk Ls 1 4 
¢ ye : a . > = a 
between two parallel walls of distance / the interference between XY 0,25 o 5 
S \ Tunnel Side -Walls 
aie the blades and the walls causes a secondary flow and thus an z | 
, additional loss. We assume that this secondary flow loss is § 0,20 
c } . . , " 6 T 1 
. independent of the aspect ratio //c, if the two interference zones 4 ; 
22465 , haa " Blade Middle -Section | 
it the walls have no connection in the middle-plane of the cas % t 
: ; 4 ° f 1 1 3 5 
cade. In this case we suggest the following formula for the a 0,15 * nm or ge = , 
ndaary 7 
mean loss coefficient: ’ A.) 
1, pp S y.) 
° | 
] F e-.. 8 | 
= } sS=C¢ + . 1 (1 | | 
Ftotat = Croc(s) ¢ i 5 wa a SE en Cee} —_ , 
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| 
* These investigations have been reported briefly by H. Schlichting at the : 
IAS Meeting in New York on Jan. 26, 1954 0,05 a ti. -dim 
+t The dimensionless loss coefficient { is defined by Ah/(qi+Ap), where | | 
Ah means the loss of total head, averaged over one spacing, q; the dynamic | | 
pressure of the inflow velocity, and Ap the pressure drop across the cascade = 1 | 
rhis is nearly equal to ¢ Ak where is the dynamic pressure of the 05 1 Y 
sis ¢ i é n/ q re ge 1s 4 ress > oO 
: 2 | a ry ar 
outflow velocity 
z 
Wall Distance i 
Fic. 2. Local loss coefficient fi.< against dimensionless wall 
distance z/c, for three different aspect ratios //¢ The measure 
ments of the left-hand and right-hand side of the blade are 
plotted together d/c = 1.0; Bg = 35 
inte 
, La 45° v2 In order to check Eq 1) we carried out experiments on turbine 
fice cascades as shown in Fig. 1 with aspect ratios //c = 1, 3, 5, the 
ibre = solidity ratio being d/c = 1, with blade angles between 8g = 50 
and 35 The flow was incompressible, the Reynolds Number 
a being Res = wec/vy = 8.104 
lers d In Fig. 2 the local loss coefficient is given for different aspect 
ratios //< Plotting the local loss coefficient against the dimen 
gs at sionless wall distance s/c yields the same curve for all aspect ra 
835, : tios //c. This proves the above assumption that the secondary 
flow loss, as given by ¢'wai, is independent of the aspect ratio, 
ull thus confirming the validity of Eq. (1 At the wall the local 
loss coefficient climbs to about four times the two-dimensional 
loss coefficient, whereas at a distance from the wall of about one 
chord the local loss coefficient comes down to the two-dimen 
sional loss coefficient This is to compare with the boundary 
layer thickness at the side-walls in front of the cascade, which 
amounts toé = 0.2 
Fic. 1. Turbine cascade. Solidity ratio d/c = 1.0; blade angle In Fig. 3a the secondary flow loss coefficient {'weu, as defined 
between Bg = 50° and 35°; angle of inflow 8; = 45°. by Eq. (1), is plotted against the blade loading Aw, /waz Here 
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Fic. 3. (a) Secondary flow loss coefficient, ¢‘war, against 
blade loading Awy/war. (b) Relative contribution of the two 
dimensional losses and the wall-losses to the total loss 


Aw, is the difference of the tangential velocity components 
across the cascade and we, is the axial velocity component, see 
It is remarkable that the coefficient (’,a1) is nearly inde- 


The influence of the blade loading 


Fig. 1. 
pendent of the blade loading 
on the secondary flow loss coefficient is caused by the variation 
of the pressure distribution over the blade. But in the present 
case this variation of the pressure distribution with the biade 
loading is very small and is confined to the vicinity of the nose 
because of the large pressure drop across the cascade in our tur 
bine cascade. For this reason ¢'waiz is nearly independent of the 
blade loading in the present case 

But for a compressor cascade it is to be expected that (' yar 
depends considerably on the blade loading. The influence of the 
spacing of the blades on €'wa1. has not yet been explored. Fig 
3b shows how the total loss splits up into two-dimensional losses 
and secondary flow losses. For aspect ratios //c between 1 and 
2, the two-dimensional losses and the secondary flow losses are 
nearly equal. 

Eq. (1), 
only if, with a variation of the aspect ratio //c, the Reynolds 
Number and therefore also the blade chord are kept constant, be- 
cause Ftwo-dim. aS Well as €’wars depend on the Reynolds Number 


as confirmed by our experiments, is, of course, valid 


In practical applications it often occurs that the through-flow 
area, and therefore the length of the blades, is fixed, but the 
In this case varying the aspect ratio 
There- 


blade chord may be varied. 
of the blades means a variation of the Reynolds Number. 
fore it may happen that the increase of the total loss coefficient 
ftotat, Which is caused by an increase of the blade chord c accord- 
ing to Eq. (1) is compensated by the decrease of ftwo-dim, and 
€'watt, Which results from the increased Reynolds Number (com- 


pare reference 4). In such a case, when applying Eq. (1), the 
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influence of the Reynolds Number on ¢to-aim, ANd ©’ wats has to b 
taken into account. 
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A Mach 3.106 Two-Dimensional Adjustable 


Nozzle for Low-Density Flow* 


Lawrence Talbotf 

Assistant Professor, Mechanical Engineering, University of 
California, Berkeley 

April 19, 1954 


- A RECENT ISSUE OF THIS JOURNAL, G. H. Backer! reported on 
the design and performance of an adjustable two-dimensional 
nozzle for use in the small No. 2 Wind Tunnel of the Low Pres 
sures Project of the University of California, Berkeley. This 
note is a brief report on subsequent tests made with a larger 
nozzle of the same type. 

Low Pressures Project interest in two-dimensional nozzles arose 
out of difficulties encountered in obtaining good test section flows 
with early axi-symmetric designs, due to the focusing of disturb- 
ances on the nozzle axis. Backer’s nozzle represented the first 


attempt at producing a two-dimensional nozzle for use in low- 
density wind tunnels. 


allow detailed measurements of the flow field (a probe large 


Although this nozzle was too small to 


enough to avoid serious viscous and thermal transpiration effects 
and time constant problems tended to obliterate small variations 
in the flow), its performance appeared promising and a larger 
version for the No. 3 Wind Tunnel® was constructed 

The performance of the No. 3 Wind Tunnel nozzle was in most 
Usable test section flows were obtained 


respects satisfactory. 
this variability in 


over a considerable range of Mach Numbers; 
Mach Number makes the nozzle a valuable tool for wind tunnel 
research. The principal undesirable feature in nozzle perform- 
ance was the existence of secondary motions in the flow which, 
because of the thick boundary layers produced at low densities, 


* This work was sponsored by the National Advisory Committee for 


Aeronautics 
+t The author would like to express his gratitude to Mr 
his supervision of the shop work necessary in the fabrication of the nozzle 


G. J. Maslach for 


and for his considerable assistance in the performance of the nozzle tests 


TABLE | 
p p 
(meas- calcu- 
Po; po’, ured lated) 
Run G, ZR, microns microns microns microns 
No. Ib./hr in Hg Hg Hg M, Hg 
] 12 0.663 $,184 1,248 86.5 >. 1] 96.7 
S 12 1.690 1,642 823.2 101.6 2.50 96.1 
3 12 0.159 17,180 2.012 87.5 4.20 86.9 
} 12 0.071 38,830 2,287 73.7 5.06 67.6 
5 12 0.021 116,500 3,016 60.8 6.20 60.3 
6 12 0.014 161,500 3,350 64.9 6.53 60.5 
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air must be throttled. For the selected design conditions, the 


stagnation chamber state was calculated to be po = 3.658 mm 
Hg, with 7) = room temperature.) 
The overall dimensions of the nozzle at design setting were 
L = 11.05 in. (length from throat to exit) 
2hx = 0.663 in. (throat height) 
2h, = 4.116 in. (exit height, including boundary-layer correc- 
tions) 


b = 3.899 in. (exit width, including boundary-layer corrections 


In evaluating the performance of the nozzle, the Mach Numbers 
within the uniform core of the flow were determined from the 
impact-stagnation pressure ratio, Po’ / Po, assuming isentropic flow 
up to the detached shock wave ahead of the probe. (This assump 


tion is supported by recent tests with an axisymmetric noz- 
zle.4) Mach Numbers calculated in this fashion did not agree 
with those determined from the ratio p/p) (also assuming isen- 


tropic flow). The lack of agreement is believed to be due to the 
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presence of secondary flow transverse pressure gradients, causing 
the wall pressure measured at the nozzle exit to be differen 
from the static pressure within the uniform core. Consequently 
where jet static pressures were required in reducing the data, they 
were calculated from the isentropic relations using the Mac} 
Number determined from po’/po. Table 1 gives the conditions 
at which the nozzle was tested, and some of the more important 
data obtained at each of the test conditions 

Mach 


Similar distributions were obtained 


Number ot 


Fig. 1 shows the centerline variation in 
tained at design conditions. 
at the other nozzle block settings. A satisfactory explanation for 
the decrease in Mach Number within the test rhombus has not 
been found. The same kind of variation, though not nearly so 
pronounced, has been observed in axisymmetric nozzles,‘ so it 
appears to be an effect not unique to two-dimensional nozzles 
0.4 in. at 


exit) it is difficult to say whether the increase in Mach Number 


Because of the thick nozzle wall boundary layers (6* 


downstream of the exit plane was due to an underexpanded flow 
condition or to incomplete cancellation of expansion waves within 
the nozzle. 

Mach 


settings are shown in Figs. 2, 3, 4, 


obtained at different nozzk 
All exhibit an in 


These con 


Number ‘‘contour maps” 
5, and 6 
denting of the contours toward the corner regions 
tour shapes are evidence that secondary motions were present 
possibly of the form sketched in Fig. 7. A silk thread tuft was 
used to explore the stream, and in the sidewall boundary layer 
the maximum deviation from axial direction was estimated visu 
ally to be about 10°, but it was not possible to carry these ex 
plorations far enough to definitely establish the nature of the 
secondary flow. Secondary flows also may have been present in 


Backer’s nozzle, for as he has remarked, the smallest usabk 
probe still occupied a considerable portion of the stream and prob 
ably obscured the fine details of the flow. 

The usable test region of the flow became smaller as the Mach 
Number was increased. At the design setting the uniform flow 
region (deviations in Mach Number less than +1°%) was about 
along the 


2 in. high by 1'/: in. wide, and extended about 2 in 


axis. At the highest centerline Mach Number, 1/; = 6.53, the 


nozzle was completely filled with boundary layer. There was 
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Stall Propagation in Cascades of Airfoils 


Alan H. Stenning 

nstructor in Mechanical Engineering, 
Technology, Cambridge, Mass 

April 19, 1954 


Massachusetts Institute of 


[ HAS BEEN SHOWN by Sears! and Emmons? that stall propaga 
tion along a cascade of airfoils is possible if a time delay exists 
between changes in angle of attack and changes in cascade pres 


sure coefficient. Experimental results*® indicate that nonviscous 
inertia effects may be primarily responsible for this time delay 
Using a simple model, it is shown below that these effects are of 
the right order to explain stall propagation phenomena found in 
ixlal compressors 

Consider a cascade of airfoils lving along the y axis in the x-y 
series of channels in 


The fluid enters the 


plane. The cascade is represented by a 
parallel, with variable-area outlets (Fig. 1) 
iscade at (1) with entry angle §;, is turned to angle 82 at (3) ina 
short distance, and leaves the channel at (2 The assumption 


is made that A2/A; is a function of 8; 





indirect evidence of air condensation at this highest Mach Num- 


appears 
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jer As A =a = F(cot B 
oe 
= ~e Ao/A corresponds to the real-cascade parameter cos 
SW Ww : ; eens pins 
cas B2V (2/p) (pu — pz) obtainable as a function of 8, from test, 
~ = a where Po is total pressure entering the cascade and pf» is stream 
pressure leaving 
SOLUTION IN THE FIELD BEFORE THE CASCADI 


Perturbations are considered from a steady flow with angle 3 
ind velocity components C,, Cy 


In the unsteady flow, 


From continuity, u, + = 0: @ + © ( Solutions to this 


equation are obtained in the form 


ny nry Nw 
o = : 3 ( {) cos + sin : 1) 
= ) ) 
Satisfving the boundary condition, 
o o (Pat x 
From Bernoulli's equation, 
b, + H constant 
where 
, F “24 F 
H = + + 
») ») 
2 f 2 f 
For small perturbations from the steady flow, 
bp 
oO - €O + ) (2 
p 
Eqs. (1) and (2) are valid in the region — \ 0 
ENTRY TO THE CASCADE 
The distance (1)-(3) is assumed small so that inertia effects 


between (1) and (3) may be neglected. 
Thus, 
C"* p ( p 
» T » 
2 p 2 f 
bp ) 
6c; + 6 + 3) 
p 
At (1), Eq may be rewritten 
6p 
Q + Cy bt t ) } 


MOMENTUM EFFECTS IN THE CASCADE 


Considering the flow between (3) and (2) as 1-dimensional, the 


momentum equation in the /-direction is 


On 0 1 Op 
we | 0 
ot O p Oo 
Integrating with respect to / from to (2 
rey - tf od ( f 
L } } 0 
ol 2 2 p p 


For small perturbations from a steady flow, 
Od« _ : f 
or i ; 
rhis 


downstream 


with the assumption that p. does not change issumption 


to be valid if there is no interference. 


Eliminating (c3;6c; + 6p;/p) between (4) and (5 
O6« 
> +- (dg; 4 CadC “} 6 
, 
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From continuity, 
1 ( == A C2 
1,6¢ = 1od¢ i 2 626 | 
Alse, from continuity, « = ¢; cos Bz. Substitution for ce, cs yields 
the following equation for @ at point (1) 
Lo C.D ri : 0 1 
ad: + + — i = (0 (7 
COS Be a’ cos? 8 a?’ cos B 1. 
But, 
641 ba l da 
= = d(cot 8 
1 a a d(cot B,) 
a’ a’ ds — cot Bidy 
= d(u/v) = 
a a Cy 
Thus, 
Lébzt Crbz cot Bia’ cot?B,a’ 
ots 4 a -_ ~ -Crdy = 0 
cos By a’ cos~ Bo a a? cos? By 
(8) 


This equation must be satisfied by ¢ at the entrance to the cas- 


cade, this is, at x = 0. Substituting the solution for ¢ into (8) 
r — 4 nny nny ‘ 
and separating coefficients of cos »”’ sin two simultaneous 
) } 


differential equations are obtained in a, and b 


Assuming solutions of the form 
Gd, =A nernt 


b, = Bed 


the characteristic equation 


1954 


Lnx nw C ( cot p;- 
cos $82 + An T : — 
a b b a? cos Bo a 


nw cot? Bia’ 
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D a’ cos Be 


is obtained. 
Roots with positive damping are obtained when a’ a/cot 
8,). Stall propagation is possible when a’ = (a/cot 8 At this 


point, the characteristic equation becomes 


nr cot py ’ 
b a’ cos Be 
A?+ = 0 
Lnr 
: + COS Pe 
D 


representing waves of frequency 


since ¢, = ¢3; cos B 


The velocity of wave propagation of the mth harmonic along 


ss w 2b 
the cascade LU’, = 
2n n 
cot 9, -¢ 
= a? 
l 
Luz 
+ cos Bo 
D 
l cot 2; 
( Lur 
a ~ t+ cos g 
D 
For a cascade with 8; = 45°, 6 25°, a = 0.9 (a common value 


for the first harmonic is 0.3 for 6/L 
1.35 for b/L 


range of experimentally obtained propagation velocities 


at stall), the ratio u/< 


and tends to very large. These values cover the 


INTERFERENCE EFFECT OF A BLADE Row UPSTREAM OF THE 


CASCADE 


Upstream interference effects are very easily determined. If 


a row of closely spaced blades of large chord is situated distance 
p from the cascade (Fig. 2), the perturbation ¢, will be forced to 
zeroatx = —p. 


A solution for @ may be found in the form 


unx un > 7 
x+2p 
b h 


L< T ¢ J 


Substitution of this solution into Eq. (8) as before gives, for 


the case of undamped waves, when a’ = (a/cot ;), 


mr\ cot Pi 2nxp 
h a? 1+ ¢ b a 
Lur Z 
a eae t+ COS Bo FF ig 


Interference increases the velocity of stall propagation for 





of order unity but does not change the velocity of propagation for 
large b/L. 

For wave lengths of the same order as the blade chord, the 
boundary layer time delay should be approximately the same as 
the nonviscous inertia delay of the free stream, and, in conse- 
quence, the theory should overestimate stall propagation veloci- 


ties. However, when the wave length is many blade chords, in- 
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effects will dominate the flow, and the theory should repre- 


ert! 
sent a good approximation to reality. 
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Pressure Distribution at Mach Number One 
over a Flat Plate at Angle of Attack 


Robert Wesley Truitt 

Professor of Aeronautical Engineering, Virginia Polytechnic Institute, 
Blacksburg, Va. 

Aay 5, 1954 


I A RECENT PAPER, Guderley determined the pressure distri- 
bution over a flat plate in two-dimensional flow with a free- 
stream Mach Number one.! It is interesting to note that a very 
simple approximate solution of this same problem was found 
previously (see Fig. 5 of reference 2) 

Using the notation of reference 1, the results of reference 2 are 
compared with those of Guderley in Fig. 1. The pressure distri- 
butions on the pressure side are seen to be quite similar in shape 
The numerical values of reference 2 are slightly larger than 
Guderley’s values but the order of magnitude of the pressure is 
the same. The center of pressure at the pressure side is at 38 
per cent from the results of reference 2; this compares favorably 


with Guderley’s value of 35 per cent 
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On the upper (or suction) side of the flat plate, the general 
character of the pressure distribution curves for the two results is 
In reference 2, it was assumed that 
1 Meyer 


expansion which deflects a sonic parallel flow by an angle equal 


seen in Fig. 1 to be dissimilar 
the pressure on the suction side could be obtained by 


to the angle of attack. In terms of Guderley’s pressure param 


eter, the pressure on the suction side, according to reference 


supersonic-type 


pres 


that 


gives /((3/2)6 = |, ie., a constant or 


sure distribution. As noted by Guderley,! this is the value 
would be obtained at the suction side of the plate if the flow ex- 
tended to infinity beyond the trailing edge or, in other words, 


this is the value of the asymptot of Guderley’s curve n/[(3/2)6 
versus x on the suction side. Guderley’s pressure distribution on 
the suction surface appears to be more reasonable than that of 
reference 2. The results of Guderley show the center of pressure 
of the suction side to lie at 46 per cent while reference 2 shows it 
Guderley found the center of 


to be, of course, at 50) per cent 


pressure of the entire wing to be at 42 per cent compared to 44.5 
per cent found in reference 2 

Guderley’s value of the total lift is slightly larger than that 
found in reference 2, due primarily to the peak pressures at the 
upper side of the leading edge. It is interesting to note, however, 
that there is less than 7 per cent difference between the predicted 


total lifts 
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On Slender-Body Theory and Apparent Mass 


Alvin H. Sacks 

Aeronautical Research Scientist, NACA, Ames Aeronautica 
Laboratory, Moffett Field, Calif 

May 12, 1954 


tbe USEFULNESS of the apparent mass concept for applica- 
tions to slender-body theory (as employed by Munk! for bod- 
ies of revolution and by Jones? for flat pointed wings) has been the 
subject of some speculation ever since the appearance of Ward's 
analysis® of a general class of slender bodies. The purpose of 
this note is to establish the relationship between Ward's result 
and the apparent mass integral and thus to clarify the place of 
apparent mass in slender-body theory 

For the sake of simplicity, the discussion will first be limited 
to the lift on a slender body in steady flow Ward's result for 


this lift is (in the notation of Fig. 1 


mad U’ pf ¢ dy ] 
x l 


where ¢ is the perturbation velocity potential which satisfies the 


L= 


two-dimensional Laplace equation near the body in planes x 
and the contour integral is taken round the body cross sec 


const 

tion With axes fixed in the body (see Fig. 1), the line of centroids 
of the cross sections (the ‘“‘camber line’’) is given by y-(x) and 
s(x). If we assume at first that y.(x) = 8 = O (where 8 side 


slip angle), we can represent ¢ as the sum of two parts, one corre 
sponding to a rigid-body translation of the cross section as a 
whole in the s direction, the other corresponding to variations in 
This 


latter portion of ¢ is required in order to satisfy the boundary 


the x direction of the shape and size of the cross section 


condition at the body surface * Thus, 
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where 
W. = velocity of centroid of cross section in s direction 
. dz, = 
= —Ula— (see Fig. 1) 
dx 
¢; = potential for unit velocity of cross section in s direction 
¢’ = potential corresponding to variations of cross section 


with x 


With these substitutions, Eq. (1) becomes 


BS dz. 
es: fp ~~ 
dx / x l ey 
pi B gy’ dy (2 


x l 


and it is noted that the quantity [a — (dz./dx)| x 


jis simply the 
angle that the line of centroids makes with the flight velocity in 
the xz plane at the base of the body (i.e., the angle of attack of 
the base section). 

At this point we recognize p { ¢g; dy as the quantity defined 
in incompressible flow theory as the additional apparent mass of 
the cross section in the z direction for pure translation in the 
direction. An alternative form of this apparent mass is the ex- 
pression — of vsldes ‘On) ds which represents the kinetic energy 
associated with a unit velocity of the cross section in the s direc- 
tion. It is essential to note, however, that it is oniy for rigid- 
body motions of the cross section (as represented for exaimnle by 
¢;), that the two integrals are identical, since only for such motions 
does 0g/On = —(dy/ds) at the boundary. Thus p f gy’ dy 
—p f g'(O¢g’/On) ds even if ¢’ is given proper dimensions by 
dividing by a velocity. 

Now, denoting the apparent mass in the s direction for pure 
translation in the z direction by m;;, we can rewrite Eq. (2) as 


. dz- ’ 
L = U,? mal a ) — pU,y p ¢'dy (3) 
dx <=] seal 7 


Without 
namely, 


and we recall that y.(x) and 8 were assumed to be zero. 
this assumption, Eq. (3) would contain another term; 
—l 2} m32[8 + (dy-/dx)] a 1 Where #32 is the apparent mass in 
the s direction for pure translation in the y direction. Thus, two 


general conclusions can be drawn: 


SCIENCES—OCTOBER, 1954 
(1) In all cases a differentiation with respect to @ yields 
OL 
Oa 
so that the lift-curve slope is, in all cases, given directly by the 


apparent mass ™;; of the base cross section 


(2) In cases where f gy’ dy = 0), the lift itself is giver 
x ‘ 


by the apparent masses of the base cross section and the angles 


that this cross section makes with the flight direction 
Let us consider first a cambered body with no wings whose 


cross sections near the base are all circular. In this case 


so that gy’ dy = 0 and the apparent mass gives the lift 


x=l 
as well as the lift-curve slope. On the other hand, consider a 
flat plate having deflected partial-span flaps. Here (when th 
plate is aligned with the flight direction), the entire potential 


‘ 


at the trailing edge is given by ¢’,=) whose integral 


¢’ dy does not vanish Hence, in this case, the lift-curve slope is 


given by the apparent mass but the lift is not. This is also true 
for a body of revolution having flat-plate wings at incidenc 
to the body. 

The above discussion can be extended to the unsteady case by 
starting with more recent results obtained by Miles‘ for unsteady 
flow rather than with Ward’s result [Eq. (1)] for steady flow 
It can then be shown that (1) in all cases, all the stability deriva- 
tives (except drag) for rigid slender bodies can be obtained from 
the apparent masses (or, more generally, the “‘inertia coeffi 
cients’’), and (2) if certain other integrals involving ¢g’ vanish, 
then all the forces and moments (except drag) are likewise ob 
tainable from the inertia coefficients of the cross sections 
Further details can be found in reference 5 
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On the Nonlinear Thickness Effect for a 
Supersonic Oscillating Airfoil 


John W. Miles 

Assuciate Professor, Department of Engineering, 
California, Los Angeles 

May 10, 1954 


University of 


— had occasion to extend to more general motions Van 
Dyke’s calculation! of the nonlinear effect of thickness on 
the pressure distribution over an airfoil executing a pitching 0os- 
cillation in a supersonic flow, we present the results here 

In Van Dyke’s notation, let the displacements of the upper 
and lower surfaces of the airfoil be specified by 


y*(x, t) = g(x) + RI {h(x)e™}, y 2 0 


where 2g(x) is the thickness distribution, h(x) the amplitude 
(not necessarily real) of the oscillatory motion, w the circular fre 
quency, ¢ the time, and x a dimensionless, chordwise coordinate 
having the values 0 and 1 at the leading and trailing edges, re 
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spectively. It is assumed that the airfoil is closed, viz. 
(0) = g(1) = O (Z) 


Then, the pressure jump across the airfoil, positive up, is given by 


9 7 wtl . 
2p l 2R/} f(x e’ ; (0 


\| 


Ply=y- —P 


where 


B = (M? — 1)’ (6 
V/ being the free stream Mach Number, y the usual specific heat 
ratio, and & the usual reduced frequency, viz 


airfoil chord - w “ 


free stream velocity 

We remark that if the airfoil thickness is 0(6) and its steady 
angle of attack is a then the terms neglected in the approximation 
1) are O(67h, ah, h Terms of O( ah, h?), although entering the 
calculation of the pressure on either surface of the airfoil, make 
no contribution to the pressure jump in virtue of symmetry con- 


siderations. 
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Nonlinear Thickne Ee ffe 


A Theory of Shock-Wave Boundary-Layer 


Interaction* 


Theodore R. Goodman 

Aerodynamic Research Department, Cornell Aeronautical 
Laboratory, Inc., Buffalo, N.Y 

May 21, 1954 


HE problem to be considered is that of a linearized shock 

wave (a Mach wave) impinging on a laminar boundary layer 
along a flat plate. The early investigations of Howarth,' 
Tsien-Finston,? and Lighthill,* while each one improved the 
mathematical model of its predecessor, were not in accord with 
experiment as to the prediction of the distance which disturb- 
ances, caused by the interaction, were propagated upstream via 
the subsonic portion of the boundary layer. Each one of these 
investigations presupposed an inviscid fluid, and it became clear 
that it is this very assumption which inhibited the solution and 
caused the errors in the calculation of the distance of upstream 
influence. 

More recent investigations which include the effects of vis 
cosity were far more successful in predicting the results of ex 
periment. For example, Kuo‘ considered the boundary layer 
to be divided into two parts—the inner layer near the wall being 


* This is a precis of the author’s Ph.D. thesis submitted to the Graduate 
School, Cornell University, September, 1953 

Copies of the complete thesis may be obtained by writing to Elma 1 
Evans, Librarian, Cornell Aeronautical Laboratory, Inc., Buffalo 21, N.Y., 


and requesting Report No. CAL-61. 
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viscous and incompressible, and the outer layer being viscous and 


compressible. The conclusions are confirmed favorably by 
experiment. 

In an attempt to modify and correct his earlier theory, Light 
hill himself has proposed a more sophisticated mathematical 
model. The present study is also an attempt to improve the 
early work of Lighthill, but viscosity is introduced into the model 
in a different fashion leading to different results 

In both of these studies, the boundary layer is divided into two 
lavyers—the sublayer near the wall where viscosity acts, and a ro 
tational but nonviscous layer beyond. Lighthill assumes that 
in the sublayer the viscous forces are of the same order of mag 
nitude as the inertia forces, but, because the velocity is low, com 
pressibility is negligible. He is thus led to use the Orr-Sommer 
field equation for small disturbances in an incompressible viscid 
fluid in the sublayer. Certain assumptions then lead him to ‘‘d« 


rive’ the height of the sublayer 


In the present paper the assumptions are quite different. In 
the sublayer the speeds are assumed so small that the inertia forces 
are negligible with respect to the viscous forces. Hence, the 
equations of motion become Stokes’ equations for creeping flow 
Two other physical assumptions are made, viz., that the sublayer 
is small compared to the complete boundary layer, and that it can 
actually be expressed as certain fraction of the boundary-layer 
thickness 

The success of the theory depends to a considerable degree on 
estimating the thickness of the Stokes’ flow region. If one were 
to solve the undisturbed boundary-layer problem by using Stokes 
equations and compare the results with the Blasius solution it 
would be found that they agree in a small neighborhood near th« 
wall. In order to fix this distance more precisely it is chosen 
to be the momentum thickness of the boundary layer, which 
is about one-seventh the boundary-layer thickness 

At the wall the no-slip condition is applied; at the interfac« 
between the sublayer and the rotational layer continuity of the 
pressure and velocity components is imposed. Similar condi 
tions are applied at the interface between the rotational layer 
and the potential flow beyond 

Only the case of an impinging wave on a laminar layer is con 
sidered and fundamentally two phenomenon are studied in this 
connection, viz., the distance of upstream influence, and the criti 
cal shock strength necessary to cause separation 

The ratio of local to free stream static pressure is found to d¢ 

A p 2yeM? 
cay exponentially along the wall: =l1l-— 
r 8 

where x is distance along the wall measured negative upstream 


from the point where the shock first impinges the boundary 


layer. According to the present theory the quantity wy is found 
to be 
34a 1 aed 
w ! = R 
B T; 
where 
v = is distance from leading edge to point of impinge 
ment 
8 = VY M?—-1 
V = free stream Mach Number 
T,,/T; = ratio of wall temperature to free-stream temperature 
R = UU, x)/m is the free-stream Reynolds Number 
U7 = free-stream velocity 
v = free-stream kinematic viscosity 
Y = ratio of specific heats 
€ = angle of the wedge which induces the shock 


This may be compared with the same quantity as computed 
by Lighthill> (x wo ) 


48x T \+4 
cl = — ke 
p’* Ti 
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A comparison between theory and experiment for free 


Fic. 1. 
stream Mach Number (17,) 2.05. 


It is interesting to compare these values with those obtained 


by Kuo,‘ by Kuo and Ritter,’ and by Lees’ (x wo r) 


A“ wR ° 2 ‘ (Kuo) 


1/ 


a i _ (Kuo and Ritter) 


Krk a RB (Lees) 


Of the five, Lees’ result and the result of Kuo and Ritter give the 
slowest variation with Mach Number and Reynolds Number, 
while the present result gives the most rapid variation. 

As Lighthill has stated,’ his result predicts a distance of up- 
stream influence of about ten boundary-layer thickness for typical 
values of Rand M. The present result will predict about forty 
or fifty boundary-layer thicknesses 

In Fig. 1 the theoretical pressures along the wall are compared 
with two experimental curves taken from reference 8. Both 
curves are for W= 2.05 but the value of R and e is different in 


The agreement for «€ = 1° is quite good, whereas 


each case 

when « = 3° the agreement is less favorable. It is quite possible 
that the shock in the latter case is strong enough to cause the 
boundary layer to separate. To determine the possibility of this 
occurrence the phenomenon of separation was studied, and it was 
found that the slope of the velocity profile vanished somewhere 


along the wall if « > e-, where 


- (z y" 
é,» = 0.00788 ae 
T» 


This result was derived by applying a formula valid far away 
from the point of impingement at the point of impingement. 
Hence, it can only be said that if € > €-, no separation occurs, and 
if ¢ < €-, separation may occur. 

For WV = 2.05 (the case given in Fig. 1) €.- = 0.4° which, being 


below both the test values, discloses nothing. 
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Effect of Cable Model Mounting on Wind- 
Tunnel Static Pressures at High Speeds 
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SUMMARY 


It is shown that a typical cable mount for a wind-tunnel model can result 
in an effective increase of the Mach Number of 0.02 and a decrease in tunnel 
static pressure of 2 per cent at a clear jet Mach Number of 0.90, with even 


larger effects occurring at higher Mach Numbers heoretical and experi 


mental results are compared, and charts are presented to facilitate cor 


rections for any cable installation 


DISCUSSION 


ie ORDER TO MEASURE the static pressure in the wake of a simu- 
lated jet engine or behind or near the base of any body, it is 
sometimes necessary to support the model by cables to avoid a 
large effect from the usual sting setup. When a cable setup is 
used, a special application of the blocking corrections may be 
applied to allow for the effect of the cables, as follows: 
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SET-UP FOR MEASURING EFFECT OF CABLES ON 





TUNNEL STATIC PRESSURE 





A body in a wind tunnel has an effective increase of dynamic 
pressure due to its displaced volume and a second increase due to 
its wake. The volume effect diminishes rapidly with distance 
downstream, and may be neglected for a cable mount. The 
wake effect becomes constant a few diameters downstream, and 
must be allowed for. 

The wake effect for a body is 


AV 1S. 
«= = Cp (1 
V. 4C 
where 
C = tunnel area, sq. ft. 
S = model or cable reference area, sq.ft 
V. = clear jet velocity 
Cp = model or cable drag coefficient 


For a typical swept cable installation, special attention must be 
paid to the treatment of the cable reference area and the drag 
coefficient. First of all, since the wake blocking correction is in- 
herently a drag correction, it is reasonable to account for the drag 
reduction due to sweep by reducing the cable reference area (the 
total cable length times its diameter) by the (cosine)? of the sweep 


angle A. The cable drag coefficient may be found from Fig. 1 


which is an adaptation of Fig. 11:6 of reference 1; the change 
being merely to use .\/, cos A as a parameter rather than .l/ 
This action makes an allowance in line with the cosine assump 


tion above. It is of interest (see Fig. 11:6 of reference 1) that 
above VV = 0.5 the drag coefficient of a circular cylinder is domi 
nated by Mach rather than Reynolds Number effects 

Since Mach Number effects are so strong, it seems reasonable 
to assume that the drag coefficient of a cable is the same as that of 
a circular cylinder 

A final and important change to Eq 
body the wake blocking is double that at 


case of swept cables, we then have 


1 is that downstream of a 
the body For the 


1 S cos? A 
«= , ( - 
2 ( 


The velocity increment at high subsonic Mach Number becomes 


and the corrected Mach Number is then? 
U= M,[1 + (1 + 0.2.11, 6 } 


when y = 1.4 
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COMPARISON OF THEORY AND EXPERIMENT FOR CABLE 
BLOCKING CORRECTION 








The corrected static pressure is then 


L + 02M? f 7 
> = (¢ 
: si 1 + 0.21 


An experimental check of the above assumption was made by 


installing a 2'/» in. static pipe in an 8'/» by 12 ft. nearsonic wind 
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EFFECT OF SEVERAL CABLE INSTALLATIONS ON 
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EFFECT OF SEVERAL CABLE INSTALLATIONS ON TUNNEL 
MACH NUMBER 





tunnel and reading the static pressure distribution with and with 
out a cable mount (Fig. 2 The correction outlined above was 
then applied and ratios of static pressure cables-in to static pres 
sure cables-out were then computed for the theoretical and ex 
perimental cases. The eight cables were 75 in. long, '/s in. diam 
eter, swept an average of 29°15’ to the tunnel centerline 

rhe agreement between theory and experiment is shown in Fig 
4, and is in general, good. It is noted that the cable effect be 
comes quite appreciable at high Mach Numbers, and represents 
about 0.02 in Mach Number or 2 per cent in g at 1/7, = 0.90 

Corrections for the general case are presented in Figs. 4, 5, and 
6. It should be noted that these corrections are for the cables 
only, and customary wake and solid blocking corrections must 
be applied for the model in addition 

The best procedure to follow is to compute the cable effect 
according to Eqs. 2 and 8, and add it to that for wake and solid 
blocking for the model. Eqs. 4+ and 5 could then be used to com 
pute the corrected Mach Number and static pressure for the 
model, or conversely to compute the blocked conditions to meet 


desired test conditions 
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Approximate Structural Analysis by the 
Method of Collocation 

Bertram Klein and Hugh L. Cox 

trength Engineers, Douglas Aircraft Company, Inc 
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May 24, 1954 


anta Monica, 


A’ IMPORTANT SIMPLE METHOD of approximate structural 
inalysis, which appears to have been overlooked by many 
urcraft structural engineers, is the method of collocation 
[his method consists essentially of satisfying a given differential 
equation, or set of equations, at a finite number of points. Engi 
neers are able to combine practical judgment with the mathe 
matics of the method in such a manner that a reasonable solution 
toa complicated problem may be obtained with a relatively small 
mount of computational work 
It is assumed that the differential equations for a structure 
subjected to given loads are known. The governing differential 
equation, or equations, can be obtained by various methods for 
numerous types of structures.° A deflection function, or set of 
functions, which satisfies the boundary conditions of the structure 
is selected. (For purposes of discussion, deflection functions are 
however, stress functions or combined functions also 


This deflection function, which contains gen 


considered ; 
may be treated 
ralized coefticients, is substituted into the governing differential 
equation for the structure For stress distribution problems, 
the unknown generalized coetticients are determined from the re 
sulting equation by substituting coordinates for as many points 
is there are generalized coeflicients. For instability and free 
vibration problems a characteristic equation is determined 

Often the engineer desires only a rough approximate answer 
to a problem that is difficult to solve by more formal methods; 
under such conditions, a reasonable deflection function that con 
tains only one unknown generalized coefficient may be used 

An example will illustrate the speed and simplicity of the 
method. Consider a square clamped plate that is axially com 
pressed along all edges by uniformly distributed edge thrusts 
Coordinate axes, x and y, are taken along two adjacent edges of 


the plate. It is desired to find the lowest critical value of the 
The governing differential equation for the plate 


_f Ow O*u 
DAlw = N + (1 
Ox? oy? 


deflection positive downward, D = plate stiffness pet 
unit of length, and .V The fol 
lowing deflection function, which satisfies the boundary condi 


buckling load 


is‘ 


where u 

= axial load per unit of length 

tions and contains only one generalized coeflicient, is assumed: 
w=aX) é 


= length and width of plate» 
The 


mav be written immediately bv use of 


where a = generalized coetticient, a 
Y = x(x — a) sin (ax)/a, and V = y(y — a) sin (ry)/a 
desired derivatives of w 
Leibnitz’s Rule? as 


O72 d*X | mx 2a wx 
= a hi =aj2sin + 2x 1) cos 
Ox dx* ’ ( a a 
T wx rs 
x(x — a)sin } 3 
a* a 


FORUM 719 


Oru d?) 
= a .X } 
Oy dy 
0% eX ¢ 
a ) 
Ox*Oy" ax ay 
2) dix n wr tx 
a Yy=a 12 sin x 
Ox dx 
1 TT T 
2x a COS t \ \ 1 sin ) hy 
O*u d‘*) 
a JX 7 
Ov dy 
Eqs. (3) through (7) are evaluated at the point a y a/2 


Substitution of the results into (eq. | 


gives 


from which 


r*D 8/1 + 8S + Tt } n*D 
Ver - : = 5.05 


Phe correct answer® is 5.3 (a?D/a? Thus, the per cent error in 


taking one generalized coefficient in the above solution is 4.7 per 


cent If enough generalized coetlicients are taken in the deflec 
tion function, greater accuracy, of course, can be obtained 

Phe above example illustrates several important features of the 
collocation method If the governing differential equation for 
the structure is not an integral differential equation, then no 
integrations involving the assumed deflection functions are neces 
sary in order to obtain a solution of the problem This fact may 
be important if the chosen deflection functions are of a type that 
ire difficult or impossible to integrate in closed form. In energy 


methods the functions must be differentiated and integrated 


and quadratic forms are present even in linear problem: rhere 
fore, if the functions involved are not orthogonal the integrations 
may become excessive and laborious rhe collocation method, 
however, is not affected by the use of nonorthogonal functions, 
and linear terms appear throughout in linear problems 

Certain nonlinear differential equations may be solved rapidly 
by the collocation method especially when certain knowledge 
regarding the physical representation of the equations is avail 
able 


The method may merit when used in conjunction with 


have 


experimental work. Often experiments reveal the types of de 
flection curves that are to be expected from certain loading con 
ditions When the general shapes of the deflection curves are 
obtained, it may be possible to express these curves mathemati 
cally by deflection functions so that solutions may be obtained for 


structures similar to those being tested 
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On Buckling of a Clamped Rhombic Thin Plate 
in Shear 


M. Hasegawa 
Shimane University, Matsue, Japan 
July 9, 1954 


| By US CONSIDER A CLAMPED RHOMBIC thin plate submitted to 
the action of shearing forces S per unit length of the edges 
uniformly distributed along the edges. This problem becomes 
important for the swept wings. The Ritz energy method will 
be applied to this problem 

Let the oblique coordinates axes (x, y) be taken in the middle 
plane of the plate of uniform small thickness in such a way that 
the origin coincides with the center of the plate and that the axes 
are parallel to the sides, and let the length, the thickness, the 
density, Young's modulus, and Poisson's ratio of the plate be 
denoted by 2a, h, p, E and », respectively. Let the angle be- 
tween coordinates axes be denoted by a. If w be the transverse 
displacement of a point on the middle plane, the total energy U 
(See Fig. 1.) 


of the plate is expressed by the following: 


D ea °%a j Pe) 
SU = — 4 (Aw)? 4 cos @ 
GT ae a Ox0' 
o*w \2} 
4 cos? a ; 
axdy/) § 


2AD sin a@ ie “ Ow Ow 
i dxdy 
a? J ~aee ~@ OF OF 


Eh Sa? a 0? 


- : A; = a : (2) 
12(1 — vy?) D Ox? Oy? 


dxdy + 


where 
D = 


As the edges are clamped, the boundary condition will be satis 


fied, if we use for the following expression: 


y?)? Javo + anxy + aoo(x? + y?) + 


- £*)* (a* 
d22x*y? + dai(x8y + xy) } (3) 


w = (a? 


The constants dyo, di, . must be determined so as to make 


Uaminimum. From this requirement, we have 

ou ou ou 

= 0; 0; : - =); 
O01) O09 Od22 


ou 

Ooo 
oU 
0031 


= () (4) 


Eqs. (4) are homogeneous linear in the a’s, and give in general 
doo = Qn = . = O unless A takes particular values. This means 
that displacement of all the points on the plate does not take 
place. 

In order to equilibrate in buckled state, the values of A must 
be determined in such a way that the solutions exist with some of 
the a’s different from 0. This requires that \ must be equal to 
some root of the determinantal equation derived from the co- 
efficients of ayo, au, in Eqs. (4). The smallest root A» of the 
determinantal equation corresponds to the critical value of the 
shearing stress S. The values of \» have been calculated for a = 
55°, 60°, 70°, 75°, and 90° and are given in Table 1. 

The true value of A» calculated by the author! for square plate 
is 36.13 and B. Budiansky and R. W. Connor? found that 35.73 S 
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Pd 
Va 


TABLE 1 





io 90° 


1816 36.298 


70° 


60° 
114.301 


69.9599 57 


5f 
153.513 


A» S 36.50. By comparing these values, the values in the table 


seem to be of good accuracy. 
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A 


14.4nb 
QO =- ~ ; 
(rw + 27 + 2w +1) cos o Vdr cos ¢ 
{ 
(w™-1)(1—~7) .. (w—-1)yr+(1-—~7r) é 
Vo + 9 Yy, +r¥o] + 


6 3 


TRANSCRIPTION error from my original notes was made 
The equation for ‘‘Q,’’ should read: 


j 3b 


Lv . 
0.105ra + 0.362r + 0.095 + 0.105} 


This equation differs from the equation printed in the Journal 
’ marked by the vertical arrow and by 


by the addition of the ‘‘r’ 
The 


changes in the constants marked by the slanted arrows 
same corrections apply to Eqs. (13), (18), (19), and (20) 
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